• <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精品国产麻豆蜜芽| 久久ZYZ资源站无码中文动漫| 国産精品久久久久久久| 久久综合九色综合欧美狠狠| 香蕉久久av一区二区三区| 伊人热热久久原色播放www| 欧美777精品久久久久网| 精品久久一区二区三区| 香蕉久久影院| 久久99热这里只有精品国产| 国产69精品久久久久观看软件| 久久精品视频91| 久久久国产亚洲精品| 久久成人国产精品免费软件| 久久久久久精品免费看SSS| 男女久久久国产一区二区三区 | 国产精品99久久久久久人| 97精品伊人久久久大香线蕉| 久久精品一区二区三区中文字幕 | 热久久国产欧美一区二区精品| 久久久久婷婷| 亚洲中文字幕久久精品无码喷水| 欧洲人妻丰满av无码久久不卡| 色8久久人人97超碰香蕉987| 国产真实乱对白精彩久久| 欧美丰满熟妇BBB久久久| 久久久午夜精品| 99久久99久久精品国产| 伊人久久大香线蕉综合网站| 69国产成人综合久久精品| AV无码久久久久不卡蜜桃| 久久天天躁狠狠躁夜夜不卡| 国产99久久久国产精品小说| 欧美精品一区二区久久 | 26uuu久久五月天|