hdu5045 Contest（费用流）

这个题目刚开始看没思路，听别人说是费用流，打死也没看出来，后来发现自己傻了，确实很简单，这个题是说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;
}