• <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天天拍| 久久久精品国产免大香伊| 久久天天躁狠狠躁夜夜躁2O2O| 九九热久久免费视频| 久久精品国产亚洲77777| 久久婷婷五月综合成人D啪| 久久亚洲AV成人无码国产| 久久精品亚洲乱码伦伦中文| 久久精品卫校国产小美女| 99久久国产综合精品成人影院| 久久综合偷偷噜噜噜色| 午夜不卡888久久| 久久久精品人妻一区二区三区四| 亚洲欧洲精品成人久久奇米网| www.久久99| 久久亚洲精品中文字幕| 亚洲综合精品香蕉久久网| 日本亚洲色大成网站WWW久久| 狠狠色丁香久久综合五月| 久久国产精品无码一区二区三区| 久久精品中文字幕一区 | 麻豆亚洲AV永久无码精品久久 | 亚洲综合精品香蕉久久网| 香港aa三级久久三级老师2021国产三级精品三级在 | 久久婷婷五月综合成人D啪| 久久亚洲中文字幕精品一区四 | 久久久久久国产精品无码下载| 久久精品国产半推半就| 久久国产精品久久| 国产精品美女久久久久| 嫩草伊人久久精品少妇AV| 日韩AV无码久久一区二区| 麻豆一区二区99久久久久| 国产成人久久精品激情| www久久久天天com| 国产精品伊人久久伊人电影| 久久99精品久久久久久野外| 久久久噜噜噜久久| 久久亚洲国产成人影院|