poj 3686 The Windy's

http://poj.org/problem?id=3686

 #include <cstdio>
 #include <cstring>
 #include <algorithm>
 #define maxn 30000
 using namespace std;

 const int inf=<<;
 int weight[][maxn],lx[maxn],ly[maxn],match[maxn],a[][];
 bool sx[maxn],sy[maxn];

 int n,m;

 bool path(int u)
 {
     sx[u]=true;
     for(int v=; v<=n*m; v++)
     {
         if(!sy[v]&&lx[u]+ly[v]==weight[u][v])
         {
             sy[v]=true;
             if(!match[v]||path(match[v]))
             {
                 match[v]=u;
                 return true;
             }
         }
     }
     return false;
 }

 int bestmatch()
 {
     memset(lx,-inf,sizeof(lx));
     memset(ly,,sizeof(ly));
     memset(match,,sizeof(match));
     for(int i=; i<=n; i++)
     {
         for(int j=; j<=n*m; j++)
         {
             lx[i]=max(lx[i],weight[i][j]);
         }
     }
     for(int u=; u<=n; u++)
     {
         while()
         {
             memset(sx,false,sizeof(sx));
             memset(sy,false,sizeof(sy));
             if(path(u)) break;
             int dx=inf;
             for(int i=; i<=n; i++)
             {
                 if(sx[i])
                 {
                     for(int j=; j<=n*m; j++)
                     {
                         if(!sy[j]) dx=min(lx[i]+ly[j]-weight[i][j],dx);
                     }
                 }
             }
             for(int i=; i<=n; i++)
             {
                 if(sx[i]) lx[i]-=dx;
             }
             for(int j=; j<=n*m; j++)
             {
                 if(sy[j]) ly[j]+=dx;
             }
         }

     }
     int sum=;
     for(int i=; i<=n*m; i++)
             sum+=weight[match[i]][i];
     return -sum;
 }
 int main()
 {
     int T;
     scanf("%d",&T);
     while(T--)
     {
         scanf("%d%d",&n,&m);
         for(int i=; i<=n; i++)
         {
             for(int j=; j<=m; j++)
             {
                 scanf("%d",&a[i][j]);
             }
         }
         for(int i=; i<=n; i++)
         {
             for(int j=; j<=n; j++)
             {
                 for(int k=; k<=m; k++)
                 {
                     weight[i][(k-)*n+j]=-a[i][k]*j;
                 }
             }
         }
         printf("%.6f\n",1.0*bestmatch()/n);
     }
     return ;
 }