hdu5045：带权二分图匹配

题目大意 ：

n个人 做m道题，其中 每连续的n道必须由不同的人做

已知第i人做出第j题的概率为pij，求最大期望

思路：
考虑每连续的n道题 都要n个人来做，显然想到了带权的二分图匹配

然后就是套模板了

代码：

#include <iostream>
#include <stdio.h>
#include<string.h>
#include<algorithm>
#include<string>
#include<ctype.h>
#include<math.h>
#include <cstring>
using namespace std;
#define MAXN 15
int n,m;
int m1[15];
int m2[15];
double value[1010][15];
bool xiangdeng(double a,double b)
{
    if(fabs(a-b)<0.00000001)
        return 1;
    return 0;
}

double fuck(int m,int n,double Gragh_gailv[][MAXN])
{
    int s[MAXN],t[MAXN];
    double l1[MAXN],l2[MAXN];
    int p,q,index,j,k;
    double value_return=0;
    for(index=0;index<m;index++)
    {
        l1[index]=-10000000;

        for(j=0;j<n;j++)
            l1[index]=Gragh_gailv[index][j]>l1[index]?Gragh_gailv[index][j]:l1[index];
        if(xiangdeng(l1[index],-10000000))
            return -1;
    }

    for(index=0;index<n;index++)
        l2[index]=0;
    memset(m1,-1,sizeof(m1));

    memset(m2,-1,sizeof(m2));


    for(index=0;index<m;index++)
    {
        memset(t,-1,sizeof(t));
        p=0;q=0;
        for(s[0]=index;p<=q&&m1[index]<0;p++)
        {
            for(k=s[p],j=0;j<n&&m1[index]<0;j++)
            {
                if(xiangdeng(l1[k]+l2[j],Gragh_gailv[k][j])&&t[j]<0)
                {
                    s[++q]=m2[j];
                    t[j]=k;
                    if(s[q]<0)
                    {
                        for(p=j;p>=0;j=p)
                        {
                            m2[j]=k=t[j];
                            p=m1[k];
                            m1[k]=j;
                        }
                    }
                }
            }
        }

        if(m1[index]<0)
        {
            index--;
            double pp=10000000;
            for(k=0;k<=q;k++)
            {
                for(j=0;j<n;j++)
                {
                    if(t[j]<0&&l1[s[k]]+l2[j]-Gragh_gailv[s[k]][j]<pp)
                       pp=l1[s[k]]+l2[j]-Gragh_gailv[s[k]][j];
                }
            }
            for(j=0;j<n;j++)
               l2[j]+=t[j]<0?0:pp;
            for(k=0;k<=q;k++)
               l1[s[k]]-=pp;
        }
    }
    for(index=0;index<m;index++)
        value_return+=Gragh_gailv[index][m1[index]];
    return value_return;
}
int main()
{
    #ifndef ONLINE_JUDGE
        freopen("shu.txt","r",stdin);
    #endif
    int t;
    scanf("%d",&t);
    int cas=0;
    while(t--)
    {
        double ans=0;
        cas++;
        scanf("%d%d",&n,&m);
        memset(value,0,sizeof(value));
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<m;j++)
                scanf("%lf",&value[j][i]);
        }
        for(int i=0;i<m;i+=n)
        {
            int e=i+n-1>=m?m-1:i+n-1;
            ans+=fuck(e-i+1,n,value+i);
        }
        printf("Case #%d: %.5f
",cas,ans);
    }


    return 0;
}