算法提高 金陵十三钗 状态亚索DP

问题描述

　　在电影《金陵十三钗》中有十二个秦淮河的女人要自我牺牲代替十二个女学生去赴日本人的死亡宴会。为了不让日本人发现，自然需要一番乔装打扮。但由于天生材质的原因，每个人和每个人之间的相似度是不同的。由于我们这是编程题，因此情况就变成了金陵n钗。给出n个女人和n个学生的相似度矩阵，求她们之间的匹配所能获得的最大相似度。
　　所谓相似度矩阵是一个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列。

输出格式

　　仅一行，一个整数，表示可获得的最大相似度。

样例输入

4
97 91 68 14
8 33 27 92
36 32 98 53
73 7 17 82

样例输出

354

数据规模和约定

　　对于70%的数据，n<=10
　　对于100%的数据，n<=13

样例说明

　　最大相似度为91+92+93+73=354

对于70%的数据,直接用全排列枚举出每个女生对应的人然后求解,取最大值,时间复杂度是O(n*n!),对于n<=10，是可以接受的

对于100%的数据，可以采用状态压缩的办法进行动态规划：

dp[i][j]表示现在进行到i个女人，被占用学生的状态为j的最大制。

将状态j表示成二进制，第k位为1 表示第k个学生已经被占用了。

比如j=3,二进制是00000011,表示第1和第2个学生已经被占用。

对于初始状态：

dp[0][j]=0    

转移方程是：

dp[i][j]=max(dp[i-1][k]+like[i][t]） 

其中t是j中其中1个1的位置，k是j把第t位变成0的数，把所有t找一边取个最大值再+like[i][t],就是dp[i][j] 

答案就是dp[n][2^n-1]

时间复杂度应该是O(n*2^n)

考虑到对于所有i，只和i-1项有关，所以可以减少一维空间消耗

对于n<=13是可以秒出的。

最后再贴一份刚写的新鲜的代码:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53

#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #define MAXN 15 using  namespace  std; const  int  maxzt=(1<<13);    //最大的状态数, int  dp[maxzt]; int  like[MAXN][MAXN],n; int  numberOfOne( int  num){  //num 二进制中1的个数      int  cnt=0;      while (num){          cnt+=(num&1);          num>>=1;      }      return  cnt; } int  lowbit( int  x){   //num 二进制中只保留最后一个1 如: num=20 二进制10100 返回 二进制100，也就是4      return  x&(-x); } int  posOfOne( int  num){  //num 二进制中最后一个1的位置  如: num=18 二进制10010 返回 2      int  pos=0;      while (num){          pos++;          if (num&1)              return  pos;          num>>=1;      }      return  pos; } void  work( int  x){      int  maxstatus=1<<n;      for ( int  i=0;i<maxstatus;i++){          int  nowstatus=i,t=numberOfOne(nowstatus);          if (t!=x)  continue ;    //第x位女人有x个1，不是就继续找下一个数          while (t--){              int  pos=lowbit(nowstatus);              dp[i]=max(dp[i],dp[i-pos]+like[x][posOfOne(pos)]);              nowstatus-=pos;          }      } } int  main(){      scanf ( "%d" ,&n);      for ( int  i=1;i<=n;i++)          for ( int  j=1;j<=n;j++)              scanf ( "%d" ,&like[i][j]);      memset (dp,0, sizeof (dp));      for ( int  i=1;i<=n;i++)          work(i);      printf ( "%d\n" ,dp[(1<<n)-1]);      return  0; }