shuoj 组装零件 状压dp 位运算得到状态的子集

零件组装
发布时间: 2017年7月9日 20:20 最后更新: 2017年7月10日 21:12 时间限制: 1000ms 内存限制: 128M

描述
现有n个零件，小Y花费了很多时间来收集它们，现在他想把零件拼在一起，拼完就可以召唤神龙了。已知零件之间存在相邻的关系，拥有相邻关系的零件在最终的组装结果中就是相邻的，并且组装过程中每次只能通过相邻关系来组合零件。小Y每次可以选择两个零件（也可以是两个零件块，或一个零件与一个零件块）拼起来，成为一个零件块，但要求拼接时必须在两个零件块（或零件）之间存在相邻的零件。除此之外这些零件两两之间有类似于磁力的排斥关系，当将两个零件或者零件块拼接在一起的时候，会受到两边的零件间的排斥力，排斥力的大小=两边零件的相互排斥对数*单侧零件个数的最大值（拼接完成的零件组合体中的零件之间排斥不计）。现在已知零件间的相邻关系和排斥关系，小Y自然想知道如何拼接不费力，因此需要求出将这些零件组装起来的最优方案，使得所有步骤的排斥力之和最小。

输入
第一行有一个整数T表示数据组数。（T<=20）
接着有T组数据，每组数据第一行是整数n表示零件个数。
接着依此有两个n∗n的矩阵，都只由0和1构成。(2<=n<=14)
其中第一个矩阵表示零件两两之间的相邻关系，第i行第j列为1表示第i个零件与第j个零件相邻，
第二个矩阵表示零件两两之间的排斥关系，第i行第j列为1表示第i个零件与第j个零件排斥。
数据保证矩阵根据对角线对称，并保证通过零件的相邻关系可以最终拼接完成。

输出
每组输入一个整数表示拼接过程的最小排斥力之和。

样例输入1 复制
1
4
0 0 1 1
0 0 1 0
1 1 0 0
1 0 0 0
0 1 0 1
1 0 1 1
0 1 0 0
1 1 0 0
样例输出1
6
提示
样例的组装方法是现将零件2和零件3拼起来，受到排斥力1*1

然后将零件1和零件4拼起来，受到排斥力1*1

最后将零件块{2,3}和{1,4}拼起来，受到排斥力2*2.

这题题目表达就很不好，说的是两边的排斥对数，其实是左边和右边的相互排斥对数，每一边内的不算，然后相连接还有两个条件都要满足，第一个条件是，结合的两部分必须有相邻的，而且必须有可以互斥的，有点抽象，所以要快速求出每两部分的是否可以结合，而且贡献是多少。
状压dp，预处理。
预处理出每个状态的1的数目，和每个状态的结合总数和每个状态的相邻总数。

然后转移方程就是 把一个状态分成两个子集状态，这个地方运用到一个很巧妙的位运算
for(i 状态）
j=((i-1)&i);j;j=((j-1)&i) 这个位运算就算把所有的 1的的子集状态 求出。
如果它的状态时 1010 那么j依次是 1000 0010 0000
得到的是所有的1的子集，这样就可以分成两部分

#include <bits/stdc++.h>
using namespace std;

int dp[(1<<14)+100];
int cnt1[(1<<14)+100],cnt2[(1<<14)+100];
int lin[20][20],lin2[20][20];
int num[(1<<14)+100];
int n,END;
const int inf = 0x3f3f3f3f;
void init()
{
    for(int i=1;i<END;i++)
    {
        num[i]=1+num[i-(i&(-i))];
    }
    for(int st=0;st<END;st++)
    {
        cnt1[st]=cnt2[st]=0;
        for(int i=0;i<n;i++)
        {
            for(int j=i+1;j<n;j++)
            {
                if((st&(1<<i))&&(st&(1<<j)))
                {
                cnt1[st]+=lin[i][j];
                cnt2[st]+=lin2[i][j];
                }
            }
        }
    }
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        END=(1<<n);
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                int x;
                scanf("%d",&x);
                lin[i][j]=x;
            }
        }
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {

                int x;
                scanf("%d",&x);
                lin2[i][j]=x;
            }
        }
        init();
        memset(dp,inf,sizeof(dp));
        for(int i=0;i<n;i++)
        {
            dp[(1<<i)]=0;
        }
        dp[0]=0;

        for(int i=0;i<END;i++)
        {
            for(int j=((i-1)&i);j;j=((j-1)&i))
            {
                int sx=j,sy=i-j;
                int num1=cnt1[i]-cnt1[sx]-cnt1[sy],num2=cnt2[i]-cnt2[sx]-cnt2[sy];
               if(num1) dp[i]=min(dp[i],dp[sx]+dp[sy]+num2*max(num[sx],num[sy]));
            }
        }
        printf("%d\n",dp[END-1] );
    }
}