这个题目刚开始看没思路,听别人说是费用流,打死也没看出来,后来发现自己傻了,确实很简单,这个题是说n个人,m道题,然后每答完一道题,任意两个人答题的数目的差不能大于1,所以比如说4个人,10个题,那么一定是,前四个题目每人一个,5到8题没人一个,后两个随机,但也必须每人最多做一个,这样的话其实就是跑了m/n+1次费用流,因为要以费用跑最短路,所以要把费用变成其相反数来做,还有就是注意是浮点数,不能硬套模板,需要把很多地方改成浮点值,真的好烦啊,还有就是边的数量要开大一点,我开始自认为开的挺大,结果居然一直WA。。。。。
#include <iostream>
#include <string.h>
#include <stdio.h>
using namespace std;
const int mm=66666;
const int mn=5555;
const int oo=1e9;
int src,dest,node,edge;
int ver[mm],flow[mm],_next[mm];
double cost[mm],dis[mn];
int head[mn],p[mn],q[mn];
bool vis[mn]={0};
void prepare(int _node,int _src,int _dest)
{
node=_node,src=_src,dest=_dest;
for(int i=0;i<node;++i)head[i]=-1;
edge=0;
}
void addedge(int u,int v,int f,double c)
{
ver[edge]=v,flow[edge]=f,cost[edge]=c,_next[edge]=head[u],head[u]=edge++;
ver[edge]=u,flow[edge]=0,cost[edge]=-c,_next[edge]=head[v],head[v]=edge++;
}
bool spfa()
{
int i,u,v,l,r=0;
double tmp;
for(i=0;i<node;++i)dis[i]=oo;
dis[q[r++]=src]=0;
p[src]=p[dest]=-1;
for(l=0;l!=r;(++l==mn)?l=0:l)
for(i=head[u=q[l]],vis[u]=0;i>=0;i=_next[i])
if(flow[i]&&dis[v=ver[i]]>(tmp=dis[u]+cost[i]))
{
dis[v]=tmp;
p[v]=i^1;
if(vis[v])continue;
vis[q[r++]=v]=1;
if(r==mn)r=0;
}
return p[dest]>-1;
}
double spfaflow()
{
int delta,i;
double ret=0;
while(spfa())
{
for(i=p[dest],delta=oo;i>=0;i=p[ver[i]])
if(flow[i^1]<delta)delta=flow[i^1];
for(i=p[dest];i>=0;i=p[ver[i]])
flow[i]+=delta,flow[i^1]-=delta;
ret+=delta*dis[dest];
}
return ret;
}//以上是黄大神的费用流模板
int n,m,T;
double a[15][1005];
int main()
{
cin>>T;
int t=1;
while(T--)
{
double sum=0.0;
cin>>n>>m;
memset(a,0,sizeof(a));
for(int i=1; i<=n; i++)
for(int j=1; j<=m; j++)
{
cin>>a[i][j];
}
for(int k=0; k<=m/n; k++)
{
prepare(2*n+2,0,2*n+1);
for(int i=1; i<=n; i++)
{
addedge(src,i,1,0);
addedge(i+n,dest,1,0);
for(int j=1; j<=n; j++)
addedge(i,j+n,1,-a[i][j+k*n]);
}
sum+=(-spfaflow());
}
cout<<"Case #"<<t++<<": ";
printf("%.5lf\n",sum);
}
return 0;
}