Poj(3686),最小权匹配,多重匹配,KM
2016-08-15 10:21
411 查看
题目链接
| Total Submissions: 4939 | Accepted: 2080 |
The manager wants to minimize the average of the finishing time of the N orders. Can you help him?
Input
The first line of input is the number of test case. The first line of each test case contains two integers, N and M (1 ≤ N,M ≤ 50).
The next N lines each contain M integers, describing the matrix Zij (1 ≤ Zij ≤ 100,000) There is a blank line before each test case.
Output
For each test case output the answer on a single line. The result should be rounded to six decimal places.
Sample Input
3
3 4
100 100 100 1
99 99 99 1
98 98 98 1
3 4
1 100 100 100
99 1 99 99
98 98 1 98
3 4
1 100 100 100
1 99 99 99
98 1 98 98
Sample Output
2.000000
1.000000
1.333333
POJ Founder Monthly Contest – 2008.08.31, windy7926778
2, 然后是拆点,分块,每一个机器都对应了nx个玩具。
3, 这里是最难的地方了。我也是看了结题报告才知道怎么处理这里,就是说可能一个机器加工多个玩具,比如说k个,耗时为a1,(a1+a2),(a1+a2+a3),......(a1+a2+......+ak);
那么,对于一个玩具(例如a1),他所贡献的时间其实为a1乘k,我们可以发现规律,对于倒数第k的玩具,他贡献的耗时为k*a1,那么拆点的时候,就应该是a1,a1乘2,......
当做模板好好学习吧!!!
The Windy's
| Time Limit: 5000MS | Memory Limit: 65536K || Total Submissions: 4939 | Accepted: 2080 |
Description
The Windy's is a world famous toy factory that owns M top-class workshop to make toys. This year the manager receives N orders for toys. The manager knows that every order will take different amount of hours in different workshops. More precisely, the i-th order will take Zij hours if the toys are making in the j-th workshop. Moreover, each order's work must be wholly completed in the same workshop. And a workshop can not switch to another order until it has finished the previous one. The switch does not cost any time.The manager wants to minimize the average of the finishing time of the N orders. Can you help him?
Input
The first line of input is the number of test case. The first line of each test case contains two integers, N and M (1 ≤ N,M ≤ 50).
The next N lines each contain M integers, describing the matrix Zij (1 ≤ Zij ≤ 100,000) There is a blank line before each test case.
Output
For each test case output the answer on a single line. The result should be rounded to six decimal places.
Sample Input
3
3 4
100 100 100 1
99 99 99 1
98 98 98 1
3 4
1 100 100 100
99 1 99 99
98 98 1 98
3 4
1 100 100 100
1 99 99 99
98 1 98 98
Sample Output
2.000000
1.000000
1.333333
Source
POJ Founder Monthly Contest – 2008.08.31, windy7926778
题意:
有N个工件要在M个机器上加工,有一个N*M的矩阵描述其加工时间。同一时间内每个机器只能加工一个工件,问加工完所有工件后,使得平均加工时间最小。分析:
1, 首先,这是一个最小权匹配,但是我没写过,然后看了一下网上的结题报告,又问了一下阳哥,发现,最小权匹配可以转换为最大权匹配,权值改为负数就可以了。2, 然后是拆点,分块,每一个机器都对应了nx个玩具。
3, 这里是最难的地方了。我也是看了结题报告才知道怎么处理这里,就是说可能一个机器加工多个玩具,比如说k个,耗时为a1,(a1+a2),(a1+a2+a3),......(a1+a2+......+ak);
那么,对于一个玩具(例如a1),他所贡献的时间其实为a1乘k,我们可以发现规律,对于倒数第k的玩具,他贡献的耗时为k*a1,那么拆点的时候,就应该是a1,a1乘2,......
当做模板好好学习吧!!!
Code(C++):
#include <iostream> #include <cstring> #include <cstdio> #include <algorithm> #include <cmath> using namespace std; #define MAXN 55 #define inf 0x7ffffff int w[MAXN][2555]; int lx[MAXN],ly[2555]; int linky[2555]; int visx[MAXN],visy[2555]; int slack[2555]; int nx,ny; int n,m; bool find(int x) { visx[x]=1; for(int y=1;y<=ny;y++) { if(visy[y]) continue; int t=lx[x]+ly[y]-w[x][y]; if(t==0) { visy[y]=1; if(linky[y] ==-1 || find(linky[y])) { linky[y]=x; return true; } } else if(slack[y] > t) slack[y]=t; } return false; } int KM() { memset(linky,-1,sizeof(linky)); memset(ly,0,sizeof(ly)); for(int i=1;i<=nx;i++) { lx[i]=-inf; for(int j=1;j<=ny;j++) if(w[i][j] >lx[i]) lx[i]=w[i][j]; } for(int x=1;x<=nx;x++) { for(int i=1;i<=ny;i++) slack[i]=inf; while(1) { memset(visx,0,sizeof(visx)); memset(visy,0,sizeof(visy)); if(find(x)) break; int d=inf; for(int i=1;i<=ny;i++) if(!visy[i] && d>slack[i]) d=slack[i]; for(int i=1;i<=nx;i++) if(visx[i]) lx[i]-=d; for(int i=1;i<=ny;i++) if(visy[i]) ly[i]+=d; else slack[i]-=d; } } int ans=0; for(int i=1;i<=ny;i++) if(linky[i] >-1) ans+=w[linky[i]][i]; return -ans; } void init() { scanf("%d%d",&n,&m); nx=n; ny=n*m; int cost; for(int i=1;i<=n;i++) { int cnt=1; for(int j=1;j<=m;j++) { scanf("%d",&cost); for(int k=1;k<=n;k++) { w[i][cnt++]=-cost*k; } } } } int main() { int t; scanf("%d",&t); while(t--) { init(); double ans=1.0*KM()/n; printf("%.6f\n",ans); } return 0; }
相关文章推荐
- POJ 3686 —— 最小费用流||最小费用匹配&KM
- poj 3686 The Windy's 二分图最小权和匹配KM
- poj 2195(KM求最小权匹配)
- POJ 3686 最小权匹配
- 【POJ 2400】 Supervisor, Supervisee(KM求最小权匹配)
- Poj 3686 The Windy's【KM匹配】
- [ACM] POJ 3686 The Windy's (二分图最小权匹配,KM算法,特殊建图)
- POJ 3686 The Windy's (二分图最小权匹配 拆点 构图)
- POJ 2400 KM最小权匹配+输出所有匹配方案
- poj 3686(拆点+最小权匹配)
- [ACM] POJ 3686 The Windy's (二分图最小权匹配,KM算法,特殊建图)
- poj 3686 The Windy's 二分匹配 KM算法求最小权匹配
- POJ 3686 The Windy's 最小权值匹配
- 二分图完备匹配(最小费用 || KM) poj 2195 GoingHome
- POJ-3686 The Windy's 最小权匹配
- poj 2195(KM求最小权匹配)
- POJ 2400 KM最小权匹配+输出所有配对方案
- 【POJ 2400】 Supervisor, Supervisee(KM求最小权匹配)
- poj 3686(拆点+最小权匹配)
- POJ 3686 The Windy's(拆点、最小权匹配)