算法提高 金陵十三钗
时间限制:1.0s 内存限制:256.0MB
金陵十三钗
本题难度:难
本题占分比例:5%
本题难度:难
本题占分比例:5%
问题描述
在电影《金陵十三钗》中有十二个秦淮河的女人要自我牺牲代替十二个女学生去赴日本人的死亡宴会。为了不让日本人发现,自然需要一番乔装打扮。但由于天生材质的原因,每个人和每个人之间的相似度是不同的。由于我们这是编程题,因此情况就变成了金陵n钗。给出n个女人和n个学生的相似度矩阵,求她们之间的匹配所能获得的最大相似度。
所谓相似度矩阵是一个n*n的二维数组like[i][j]。其中i,j分别为女人的编号和学生的编号,皆从0到n-1编号。like[i][j]是一个0到100的整数值,表示第i个女人和第j个学生的相似度,值越大相似度越大,比如0表示完全不相似,100表示百分之百一样。每个女人都需要找一个自己代替的女学生。
最终要使两边一一配对,形成一个匹配。请编程找到一种匹配方案,使各对女人和女学生之间的相似度之和最大。
所谓相似度矩阵是一个n*n的二维数组like[i][j]。其中i,j分别为女人的编号和学生的编号,皆从0到n-1编号。like[i][j]是一个0到100的整数值,表示第i个女人和第j个学生的相似度,值越大相似度越大,比如0表示完全不相似,100表示百分之百一样。每个女人都需要找一个自己代替的女学生。
最终要使两边一一配对,形成一个匹配。请编程找到一种匹配方案,使各对女人和女学生之间的相似度之和最大。
输入格式
第一行一个正整数n表示有n个秦淮河女人和n个女学生
接下来n行给出相似度,每行n个0到100的整数,依次对应二维矩阵的n行n列。
接下来n行给出相似度,每行n个0到100的整数,依次对应二维矩阵的n行n列。
输出格式
仅一行,一个整数,表示可获得的最大相似度。
样例输入
4
97 91 68 14
8 33 27 92
36 32 98 53
73 7 17 82
97 91 68 14
8 33 27 92
36 32 98 53
73 7 17 82
样例输出
354
数据规模和约定
对于70%的数据,n<=10
对于100%的数据,n<=13
对于100%的数据,n<=13
样例说明
最大相似度为91+92+93+73=354
题目分析:本题如果直接暴力,复杂度是n!无法全部AC,为了优化时间,采用状压dp,我们用dp[i][j]存储第i个人被替换时状态为j,这里的j是二进制的表示形式,第x个位置为1的话代表第x个人被替换了,所以状态转移方程为dp[i][j] = max(dp[i][j],dp[i - 1][j - p] + like[i][t]),p代表只有第x个人的状态,t代表第x个人的位置,求第x个人的状态用lowbit.为了节省空间,用滚动数组代替二维数组更新状态
代码:
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<string>
#include<queue>
using namespace std;
int like[15][15],dp[1<<13];
int n;
int lowbit(int x)
{
return x & (-x);
}
int numofone(int x)
{
int cnt = 0;
while(x)
{
if(x & 1)cnt++;
x >>= 1;
}
return cnt;
}
int locofone(int x)
{
int cnt = 0;
while(x)
{
x >>= 1;
cnt++;
}
return cnt;
}
void work(int x)
{
int maxstatus = 1<<n;
for(int i = 0; i < maxstatus; i++)
{
if(numofone(i) != x)
continue;
int t = numofone(i);
int now = i;
while(t--)
{
int pos = lowbit(now);
int k = locofone(pos);
dp[i] = max(dp[i],dp[i - pos] + like[x][k]);
now -= pos;
}
}
}
int main()
{
cin>>n;
for(int i = 1; i <= n; i++)
for(int j = 1; j <= n; j++)
cin>>like[i][j];
for(int i = 1; i <= n; i++)
work(i);
/*for(int i = 0; i < (1 << n); i++)
cout<<dp[i]<<" ";
*/
cout<<dp[(1 << n) - 1]<<endl;
return 0;
}