• <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

            <2025年8月>
            272829303112
            3456789
            10111213141516
            17181920212223
            24252627282930
            31123456

            導航

            統計

            公告

            統計系統

            留言簿(1)

            隨筆分類(227)

            文章分類(2)

            OJ

            最新隨筆

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            久久99国产亚洲高清观看首页| 久久久久久久人妻无码中文字幕爆| 欧美亚洲国产精品久久蜜芽| 国产精品99久久久久久www| 狠狠人妻久久久久久综合| 久久亚洲高清综合| 欧美黑人激情性久久| 国产视频久久| 久久久久亚洲精品天堂| 日本精品久久久久久久久免费| 少妇久久久久久被弄高潮| 91久久精品视频| 久久无码人妻一区二区三区| 久久国产精品波多野结衣AV| 久久偷看各类wc女厕嘘嘘| 久久久久久午夜精品| 97久久精品人人做人人爽| 麻豆亚洲AV永久无码精品久久| 国産精品久久久久久久| 国产99精品久久| 久久久久亚洲国产| 模特私拍国产精品久久| 嫩草影院久久国产精品| 久久久久国产精品熟女影院 | 日本欧美久久久久免费播放网| 精品久久久久久无码国产| 久久er国产精品免费观看2| 亚洲va中文字幕无码久久| 欧美精品九九99久久在观看| 久久久久九九精品影院| 国产激情久久久久影院| 久久久综合九色合综国产| 99久久超碰中文字幕伊人| 亚洲va久久久噜噜噜久久狠狠| 成人午夜精品无码区久久| 亚洲精品无码久久一线| 亚洲日韩中文无码久久| 亚洲成色www久久网站夜月 | 大香网伊人久久综合网2020| 国产精品久久久久久福利69堂| 国产精品久久波多野结衣|