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);
}
}
}