• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>

            pku2400 Supervisor, Supervisee 二分圖最優匹配的全部解

            題意:
            給出職工對經理的評分和經理對職工的評分。求出所有的匹配方案,使得平均滿意度最高

            解法:
            KM匹配出最優解后,根據標號遍歷得所有解。注意,調整頂標以后合法匹配邊仍然是l[i]+r[j]=val[i][j]

            代碼:

            Source Code

            Problem: 2400User: yzhw
            Memory: 408KTime: 0MS
            Language: G++Result: Accepted
            • Source Code
            • # include <cstdio>
              # include <cstring>
              # include <cstdlib>
              using namespace std;
              # define max(a,b) ((a)>(b)?(a):(b))
              # define min(a,b) ((a)<(b)?(a):(b))
              # define N 15
              int map[N][N],match[N],ln[N],rn[N],n,c;
              bool l[N],r[N];
              bool dfs(int pos)
              {
              	//if(l[pos]) return false;
              	l[pos]=true;
              	for(int i=0;i<n;i++)
              		if(!r[i]&&ln[pos]+rn[i]==map[pos][i])
              		{
              			r[i]=true;
              			if(match[i]==-1||dfs(match[i]))
              			{
              				match[i]=pos;
              				return true;
              			}
              
              		}
              		return false;
              }
              void adjust()
              {
              	int minnum=0xfffffff;
              	for(int i=0;i<n;i++)
              		if(l[i])
              			for(int j=0;j<n;j++)
              				if(!r[j])
              					minnum=min(minnum,ln[i]+rn[j]-map[i][j]);
              	for(int i=0;i<n;i++)
              	{
              		if(l[i]) ln[i]-=minnum;
              		if(r[i]) rn[i]+=minnum;
              	}
              }
              void search(int now)
              {
              	if(now==n)
              	{
              		int ret[N];
              		for(int i=0;i<n;i++)
              			ret[match[i]]=i;
              		printf("Best Pairing %d\n",++c);
              		for(int i=0;i<n;i++)
              			printf("Supervisor %d with Employee %d\n",i+1,ret[i]+1);
              		//system("pause");
              	}
              	for(int i=0;i<n;i++)
              		if(match[i]==-1&&ln[now]+rn[i]==map[now][i])
              		{
              			match[i]=now;
              			search(now+1);
              			match[i]=-1;
              		}
              }
              int main()
              {
              	int test;
              	scanf("%d",&test);
              	for(int tt=1;tt<=test;tt++)
              	{
              		scanf("%d",&n);
              		memset(map,0,sizeof(map));
              		for(int i=0;i<n;i++)
              			for(int j=0;j<n;j++)
              			{
              				int t;
              				scanf("%d",&t);
              				map[t-1][i]-=j;
              			}
              		for(int i=0;i<n;i++)
              			for(int j=0;j<n;j++)
              			{
              				int t;
              				scanf("%d",&t);
              				map[i][t-1]-=j;
              			}
              		
              		memset(match,-1,sizeof(match));
              		memset(rn,0,sizeof(rn));
              		for(int i=0;i<n;i++)
              		{
              			int maxnum=-0xfffffff;
              			for(int j=0;j<n;j++)
              				maxnum=max(maxnum,map[i][j]);
              			ln[i]=maxnum;
              		}
              		double ans=0;
              		for(int i=0;i<n;i++)
              		{
              			memset(l,false,sizeof(l));
              			memset(r,false,sizeof(r));
              			while(!dfs(i))
              			{
              				adjust();
              				memset(l,false,sizeof(l));
              				memset(r,false,sizeof(r));
              			}
              		}
              		for(int i=0;i<n;i++)
              			ans-=map[match[i]][i];
              		printf("Data Set %d, Best average difference: %.6f\n",tt,ans/n/2.0);
              		memset(match,-1,sizeof(match));
              		c=0;
              		search(0);
              		printf("\n");
              	}
              	return 0;
              }

            posted on 2011-03-08 01:39 yzhw 閱讀(265) 評論(0)  編輯 收藏 引用 所屬分類: graph

            <2011年1月>
            2627282930311
            2345678
            9101112131415
            16171819202122
            23242526272829
            303112345

            導航

            統計

            公告

            統計系統

            留言簿(1)

            隨筆分類(227)

            文章分類(2)

            OJ

            最新隨筆

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            少妇内射兰兰久久| 久久精品国产亚洲精品| 久久精品无码一区二区无码| 久久国产精品99精品国产| 美女写真久久影院| 久久丫忘忧草产品| 99久久久久| 午夜天堂精品久久久久| 国内精品久久久久久中文字幕| 日韩久久无码免费毛片软件| 亚洲精品乱码久久久久久按摩| 久久精品国产精品亚洲| 久久久久亚洲av无码专区导航| 曰曰摸天天摸人人看久久久| 久久综合色老色| 久久精品中文字幕一区| 精品人妻久久久久久888| 少妇人妻综合久久中文字幕| 国产精品热久久无码av| 精品国产福利久久久| 亚洲精品白浆高清久久久久久 | 久久久久亚洲AV成人网人人网站| 日韩精品久久久久久免费| 久久青青草原精品国产软件| 久久91精品国产91久久户| 久久一日本道色综合久久| 亚洲午夜精品久久久久久app| 久久青青草原精品国产软件| 精品久久久久久99人妻| 国产午夜电影久久| 99热热久久这里只有精品68| 精品久久久久久综合日本| 久久亚洲精品人成综合网| 日日躁夜夜躁狠狠久久AV| 日韩人妻无码一区二区三区久久| 久久精品国产亚洲av麻豆蜜芽| 欧美亚洲另类久久综合婷婷| 久久久久综合中文字幕| 中文字幕精品久久| 色婷婷综合久久久久中文一区二区 | 91久久婷婷国产综合精品青草|