HDU 3718 Similarity

HDU 3718 Similarity
km最优匹配（百度）
建图的时候对应位置的边权值加一（最后建出来的图是[A-Z]->[A-Z],而不是字符串位置对应字符串位置），km匹配完之后把匹配成功的边的权值加起来除以n就是答案。

代码：

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
using namespace std;
const int nl = 26, nr = 26;
const int N = 1e4+10,inf = 9e7;
char m_str[N],mt_str[N],str[5];
int mp[nl][nr],slack[nr],ll[nl],lr[nr],visl[nl],visr[nr],match[nr],cnt=0;

bool dfs(int k)
{
    visl[k]=1;								
    for(int i=0;i<nr;i++)
    {
        if(visr[i])continue;
        if(ll[k]+lr[i]==mp[k][i])
        {
            visr[i]=1;
            if(match[i]==-1||dfs(match[i]))
            {
                match[i]=k;
                return true;
            }
        }
        else slack[i] = min(slack[i],ll[k]+lr[i]-mp[k][i]);
    }
    return false;
}


int km()
{
    memset(ll,0,sizeof(ll));
    memset(match,-1,sizeof(match));
    for(int i=0;i<nl;i++)
        for(int j=0;j<nr;j++)
        {
            lr[j]=0;
            ll[i]=max(mp[i][j],ll[i]);
        }
    //cout << cnt++ << endl;
    for(int i=0;i<nl;i++)
    {
        while(1)
        {
            memset(slack,0x3f3f3f,sizeof(slack));
            memset(visl,0,sizeof(visl));
            memset(visr,0,sizeof(visr));
            if(dfs(i))break;
            int d = inf;
            for(int j=0;j<nr;j++)
                if(!visr[j])d = min(d,slack[j]);
            for(int j=0;j<nl;j++)
                if(visl[j])ll[j]-=d;
            for(int j=0;j<nr;j++)
                if(visr[j])lr[j]+=d;
        }
    }
    int ans =0;
    for(int i=0;i<nr;i++)
        ans+=mp[match[i]][i];
    return ans;
}



int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        cnt=0;
        int n,k,m;
        scanf("%d%d%d",&n,&k,&m);
        //cout << n << k << m << endl;
        for(int i=0;i<n;i++)
        {
            scanf("%s",str);
            m_str[i]=str[0];
        }
        m_str[n]='\0';
        //cout << m_str << endl;
        for(int i=0;i<m;i++)
        {
            memset(mp,0,sizeof(mp));
            for(int j=0;j<n;j++)
            {
                scanf("%s",str);
                mt_str[j]=str[0];
                mp[mt_str[j]-'A'][m_str[j]-'A']++;
            }
            mt_str[n]='\0';
            //cout << mt_str << endl;
            //int ans=0;
            int ans = km();
            printf("%.4f\n",(ans*1.0)/n);
        }
    }
}