首页 > 技术文章 > scau 8616 汽车拉力比赛

liuweimingcprogram 2016-01-21 13:27 原文

      上次我们过了二分图的最佳匹配,现在我们看一道题目,经典的二分图的最佳匹配题目

8616 汽车拉力比赛

时间限制:500MS  内存限制:1000K
提交次数:71 通过次数:24

题型: 编程题   语言: G++;GCC

 

Description

 SCAU车队要去参加汽车拉力比赛啦。拉力比赛中,每辆汽车需要一个驾驶员和一个导航员。SCAU车队一共有有N个驾驶员和N个导航员,每个 驾驶员对每个导航员有一个默契值,现在需要求得某种N对N的配对,使得这N个配对的默契值之和最大。 

输入格式

 多Case,每一个Case第一行是一个整数N,范围是(1<=N<=16)。接下来有N行,每行有N个数,第i行第j个数表示第i个驾驶员和第j个导航员的默契值。 测试数据会有多组,当N=0时表示数据结束,这行不用处理。

输出格式

 输出配对后默契值之和的最大值。

输入样例

1
5
2
10 20
25 5
0

输出样例

5
45
至于如何做,我就不说了,看看我以前的博客吧,思路应该很清晰的
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#define inf 99999
//保存一个我们认为的无穷大
int n;
int match[50];//表示这个[驾驶员]搭配那个驾驶员
int e[50][50];//用来保存一张图,表示i-->j的默契值
int vx[50];//用来表示 驾驶员 来找过partner
int vy[50];//用来表示 导航员 被邀请过搭配
int fx[50];//一开始是用来表示每一个驾驶员与1--n个导航员的max,但是随着程序执行,会变化哦
int fy[50];//一开始是0的。我们称这两个东西为 标干
int dfs (int u)//我要来找u的partner
{
    vx[u]=1;//u已经来找过partner啦,标记一下先,待会的km算法需要
    int i;
    for (i=1;i<=n;i++)//筛选n个 导航员  可以知道,驾驶员比较“主动”
    {
        if (vy[i]==0&&fx[u]+fy[i]==e[u][i])//这个导航员没有被邀请过..导航员为u
        //驾驶员和导航员的【标杠】值相加,有这个默契值的话,这样选出来的是最优的
        {
            //vy[i]==1;//这个导航员被邀请过了,但是,接不接受?还是另外一回事
            //不是等于等于
            //因为他还可能是被别人搭配了嘛
            vy[i]=1;
            if (match[i]==0||dfs(match[i]))//如果他没被人搭配了或者,它可以找其他人搭配
            {
                match[i]=u;//i搭配u  改变条件的
                //意思就是导航员i和驾驶员u搭配
                return 1;//搭配成功
            }
        }
    }
    return 0;//我找不到啊,后面,就会执行km
}
void do_km()//km算法
//对每一个在vx[]数组中的驾驶员,我们对他进行操作
//操作对象是每一个不在vy[]数组中的导航员,在了的导航员是不能动的,是当前的最优解
{
    int i,j;
    int d;//表示默契值减小得最多,就是,搭配最大
    d=inf;//认为他是无限大
    for (i=1;i<=n;i++)//遍历每一个驾驶员
    {
        if (vx[i]==1)//表示他是遍历过的
        {
            for (j=1;j<=n;j++)//对他进行遍历导航员
            {
                if (!vy[j])//没有vy过的导航员才操作
                {
                   if (d>(fx[i]+fy[j]-e[i][j]))//对于后面没vy的导航员,选出最小差值
                   {
                       d=fx[i]+fy[j]-e[i][j];
                       //printf ("%d\n",d);
                   }
                }
            }
        }
    }
    //对于每一个vx[]的 fx[] 减去d
    //对于每一个vy[]的 fy[] 加上d
    //这样,总值会减小哦,就是,我们有机会加入新的导航员,留下空位
    for (i=1;i<=n;i++)
    {
        if (vx[i]==1)
        {
            fx[i] -= d;
            vx[i]=0;//请0
        }
        if (vy[i]==1)//是访问过的哦
        {
            fy[i] += d;
            vy[i]=0;//情0
        }
    }
    return ;
}
void work ()
{
    int i,j;
    memset (e,0,sizeof(e));
    memset (vx,0,sizeof(vx));
    memset (vy,0,sizeof(vy));
    memset (fx,0,sizeof(fx));
    memset (fy,0,sizeof(fy));
    memset (match,0,sizeof(match));

    for (i=1;i<=n;i++)
    {
        for (j=1;j<=n;j++)
        {
            scanf ("%d",&e[i][j]);//不是无向图哦,有向的
            //就是e[i][j]!=e[j][i]  有可能相等
        }
    }
    //km算法的一部分,先初始化fx,fy
    for (i=1;i<=n;i++)//遍历每一个驾驶员
    {
        fy[i]=0;
        fx[i]= -inf;//无穷小
        for (j=1;j<=n;j++)//遍历每一个导航员
        {
             if (fx[i]<e[i][j])//默契值
             {
                 fx[i]=e[i][j];
             }
        }
    }
    for (i=1;i<=n;i++)//遍历每一个驾驶员
    {
        memset (vx,0,sizeof(vx));//每一个驾驶员都是新的开始
        memset (vy,0,sizeof(vy));//要清零,
        while (!dfs(i))//如果他找不到搭配,就实现km算法
        {
            do_km();//km完后,还是会对这个想插入的节点进行dfs的,因为他还没搭配成功嘛
        }
    }
    //进行计数
    //现在match[i]记录的是[导航员i]和[驾驶员match[i]]搭配
    int ans=0;
    for (i=1;i<=n;i++)//遍历导航员
    {
        ans += e[match[i]][i];//输入的第一个是驾驶员,第二个是导航员
    }
    printf ("%d\n",ans);
    return ;
}
int main ()
{
    while (scanf ("%d",&n)!=EOF && n)
    {
        work ();
    }
    return 0;
}
View Code

 

 

推荐阅读