poj3686

题解：

KM算法

把每一个点拆成n个

然后改变编圈

代码：

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=;
char s[],q[N][];
int cnt,a[N][N],z,e,cas,b[N*N];
int visr[N],T,x,y,exl[N],exr[N],visl[N],match[N],slack[N],n,m;
int in(char s[])
{
    for (int i=;i<=cnt;i++)
     if (q[i][]==s[])return i;
    q[++cnt][]=s[];
    return cnt;
}
int dfs(int x)
{
    visl[x]=;
    for (int i=;i<=m;i++)
     if (!visr[i])
      {
          int k=exl[x]+exr[i]-a[x][i];
          if (k==)
           {
               visr[i]=;
               if (!match[i]||dfs(match[i]))
                {
                    match[i]=x;
                    return ;
                }
           }
          else slack[i]=min(slack[i],k);
      }
    return ;
}
int km()
{
    memset(exl,,sizeof exl);
    memset(exr,,sizeof exr);
    for (int i=;i<=m;i++)
     for (int j=;j<=m;j++)exl[i]=max(exl[i],a[i][j]);
    memset(match,,sizeof match);
    for (int i=;i<=m;i++)
     {
        memset(slack,0x3f,sizeof slack);
        while ()
         {
              memset(visl,,sizeof visl);
              memset(visr,,sizeof visr);
              if (dfs(i))break;
              int d=1e9;
              for (int j=;j<=m;j++)
               if (!visr[j])d=min(d,slack[j]);
              for (int j=;j<=m;j++)
             if (visl[j])exl[j]-=d;
            for (int j=;j<=m;j++)
             if (visr[j])exr[j]+=d;
             else slack[j]-=d;
         }
     }
    int ans=,flag=;
    for (int i=;i<=m;i++)ans+=a[match[i]][i];
    return ans;
}
int main()
{
    scanf("%d",&T);
    while (T--)
     {
         scanf("%d%d%d",&n,&m,&e);
         cnt=;
         memset(q,,sizeof q);
         for (int i=;i<=n;i++)
          {
              scanf("%s",s);
              b[i]=in(s);
          }
         while (e--)
         {
             cnt=;
             memset(q,,sizeof q);
             memset(a,,sizeof a);
              for (int i=;i<=n;i++)
               {
                   scanf("%s",s);
                   a[b[i]][in(s)]++;
               }
              printf("%.4lf\n",(double)km()/n);
         }
     }
}