Hdu 3718（KM算法）

Hdu 3718

 

（1）题意：

给出一个长度为n的字符序列，字符序列是由‘A’~‘Z’组成的，

先给出一个标准的字符序列，表示正确答案的顺序，其中不同的字符数量为k。

接下来给出m个字符序列，每个字符序列长度为n，

每个字符序列中的不同字符数量不大于k，其中的字符与正确序列的字符可能不一样，

请问如果将询问序列中的字符与正确序列中的字符匹配，才能使正确率更高。

输出每个询问字符序列的正确率。

（正确率 = 匹配正确的字符数量/字符序列长度）。

 

（2）思路：

将每个序列的不同字符与正确序列的字符匹配，

求出最多可以使多少个字符成功匹配。

其实统计每个询问字符所对应的不同的正确字符的位置与数量就行了，

就相当于在询问字符与它配对的字符之间建立了一条边，边的权值是匹配成功的数量。

然后求出这张图的最佳完备匹配就好了。

 

（3）代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn = 30;
const int INF = 1e9+10;
int pre[maxn],vx[maxn],vy[maxn],lx[maxn],ly[maxn],cha,n,m,k,t1,val[10005];
int mp[maxn][maxn],v1[maxn],v2[maxn];
int MIN(int x,int y){
	return x<y?x:y; 
}
int MAX(int x,int y){
	return x>y?x:y;
}
bool dfs(int u){
	vx[u] = 1;
	for(int v=1;v<=k;v++)
	if(vy[v]==0&&lx[u]+ly[v]==mp[u][v]){
		vy[v] = 1;
		if(pre[v]==-1||dfs(pre[v])){
			pre[v] = u;
			return true;
		}
	}
	else if(vy[v]==0){
		cha = MIN(cha,lx[u]+ly[v]-mp[u][v]);
	}
	return false;
}
double km(){
	memset(ly,0,sizeof(ly));
	for(int i=1;i<=t1;i++){
		lx[i] = mp[i][1];
		for(int j=2;j<=k;j++) lx[i] = MAX(lx[i],mp[i][j]);
	}
	memset(pre,-1,sizeof(pre));
	for(int i=1;i<=t1;i++){
		while(true){
			memset(vx,0,sizeof(vx));
			memset(vy,0,sizeof(vy));
			cha = INF;
			if(dfs(i)==true) break;
			for(int j=1;j<=t1;j++)
			if(vx[j]==1) lx[j] -= cha;
			for(int j=1;j<=k;j++)
			if(vy[j]==1) ly[j] += cha;
		}
	}
	int ans = 0;
	for(int i=1;i<=k;i++) 
	if(pre[i]!=-1) ans += mp[pre[i]][i];
	double res = (ans*1.0)/(n*1.0);
	return res;
}
int main(void){
	int T;scanf("%d",&T);
	while(T--){
		scanf("%d%d%d",&n,&k,&m);
		memset(v1,0,sizeof(v1));
		t1 = 0;
		for(int i=0;i<n;i++){
			char ss[5];scanf("%s",ss);
			int x = ss[0]-'A';
			if(v1[x]==0) v1[x] = ++t1;
			val[i] = v1[x];
		}
		for(int i=0;i<m;i++){
			memset(v1,0,sizeof(v1));
			memset(mp,0,sizeof(mp));
			t1 = 0;
			for(int j=0;j<n;j++){
				char ss[5];scanf("%s",ss);
				int x = ss[0]-'A';
				if(v1[x]==0) v1[x] = ++t1;
				mp[v1[x]][val[j]]++;
			}
			printf("%.4lf\n",km());
		}
	}
	return 0;
}