• <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 閱讀(255) 評論(0)  編輯 收藏 引用 所屬分類: graph

            <2011年3月>
            272812345
            6789101112
            13141516171819
            20212223242526
            272829303112
            3456789

            導航

            統計

            公告

            統計系統

            留言簿(1)

            隨筆分類(227)

            文章分類(2)

            OJ

            最新隨筆

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            午夜精品久久久久久久久| 久久成人国产精品二三区| 亚洲欧洲久久久精品| 精品久久亚洲中文无码| 久久成人国产精品| 精品久久久久久国产三级| 一本久久a久久精品综合香蕉 | 亚洲国产小视频精品久久久三级| 性做久久久久久久久浪潮| 国产亚洲欧美成人久久片| 国产99久久久国产精品小说| 激情伊人五月天久久综合| 久久综合亚洲鲁鲁五月天| 精品久久久久久国产91| 欧美成人免费观看久久| 久久综合九色综合精品| 伊人久久久AV老熟妇色| 亚洲午夜久久久| 久久99精品久久久久久噜噜| 99久久人妻无码精品系列蜜桃| 日日狠狠久久偷偷色综合0| 久久电影网一区| 精品久久人妻av中文字幕| 久久精品国产99国产精品亚洲| 国内精品免费久久影院| 国产精品久久精品| 久久精品国产免费| 国产成人综合久久综合| 久久精品水蜜桃av综合天堂| 久久妇女高潮几次MBA| 国产精品亚洲综合久久| 亚洲国产精品成人AV无码久久综合影院 | 久久香蕉综合色一综合色88| 少妇高潮惨叫久久久久久| 国产美女亚洲精品久久久综合| 久久综合日本熟妇| 久久久久久国产精品无码下载| 久久精品国产亚洲精品| 青青久久精品国产免费看| 最新久久免费视频| 久久久久久久波多野结衣高潮|