关闭

HDU 3718 Similarity(KM算法)

36人阅读 评论(0) 收藏 举报
分类:
把对应位置的字母边练起来,权值加一,最后求二分图最大权值即可。
#include<iostream>
#include<cstring>
using namespace std;
const int maxn=30+5;
int w[maxn][maxn],n,m;
int lx[maxn],ly[maxn];
int Left[maxn];
bool s[maxn],t[maxn];
char ss[35][20005];
char tt[20005];
bool match(int i)
{
    s[i]=true;
    for(int j=1;j<=m;j++) if(lx[i]+ly[j]==w[i][j]&&!t[j]){
        t[j]=true;
        if(Left[j]==-1||match(Left[j])){
            Left[j]=i;
            return true;
        }
    }
    return false;
}
void update()
{
    int a=(1<<30);
    for(int i=1;i<=n;i++) if(s[i])
        for(int j=1;j<=m;j++) if(!t[j])
            a=min(a,lx[i]+ly[j]-w[i][j]);
    for(int i=1;i<=n;i++) if(s[i]) lx[i]-=a;
    for(int i=1;i<=m;i++) if(t[i]) ly[i]+=a;
}
void KM()
{
    memset(Left,-1,sizeof(Left));
    memset(lx,0,sizeof(lx));
    memset(ly,0,sizeof(ly));
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++)
            lx[i]=max(lx[i],w[i][j]);
    }
    for(int i=1;i<=n;i++){
        for(;;)
        {
            memset(s,0,sizeof(s));
            memset(t,0,sizeof(t));
            if(match(i)) break;else update();
        }
    }
}
int r,k,c;
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d%d",&r,&k,&c);
        n=m=26;
        getchar();gets(tt);
        for(int i=0;i<c;i++) gets(ss[i]);
        for(int i=0;i<c;i++)
        {
            memset(w,0,sizeof(w));
            for(int j=0;j<=2*(r-1);j+=2) {
                w[ss[i][j]-'A'+1][tt[j]-'A'+1]++;
            }
            KM();
            int ans=0;
            for(int i=1;i<=26;i++) ans+=w[Left[i]][i];
            printf("%.4f\n",(double)ans/r);
        }
    }
    return 0;
}

0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:50391次
    • 积分:6983
    • 等级:
    • 排名:第3257名
    • 原创:690篇
    • 转载:1篇
    • 译文:0篇
    • 评论:4条
    文章分类
    最新评论