牛客国庆集训派对Day2 - E 数据排序状压dp

最新推荐文章于 2020-10-13 22:55:35 发布

原创最新推荐文章于 2020-10-13 22:55:35 发布 · 277 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#牛客国庆集训派对Day2 #E 数据排序 #状压dp

动态规划同时被 2 个专栏收录

15 篇文章

订阅专栏

状压dp

4 篇文章

订阅专栏

本文介绍了一种使用状态压缩动态规划（状压DP）的方法来解决一组照片的评分冲突问题，通过计算每张照片与其他照片之间的相对美观度，找到使评分冲突值最小化的方案。

题目链接

题意：有n张照片，有2*n(n-1)组数据，每组数据<x,y>表示照片x比照片y好看，现要给每张照片赋予一个评分 $c_{i}$ ，记g(x,y)为数据中<x,y>出现的次数，则有冲突值计算公式f(c)= $\sum g(j,i)\left ( c_{i}>c_{j} \right )+\sum \left | g(i,j)-g(j,i) \right |\left ( c_{i}=c_{j} \right )$ ,求f(c)的最小值。

思路：由于n最大只有15，所以很容易想到是状压dp，0表示当前位置照片未考虑，1表示已经考虑在内。那么对于当前状态i，扩展到之后的状态，我们可以考虑枚举扩展的方式(也就是枚举状态i^((1<<n)-1)的子状态)。假设这些状态内的照片评分均相等且大于i状态内任一已考虑的照片的评分。设f[i]表示状态i评分均相等产生的冲突值，d[i][j]表示第i张照片评分大于j状态内任一已考虑的照片的评分时产生的冲突值，dp[i]存答案，那么转移方程为dp[i|k]=dp[i]+f[k]+d[][i](k状态里每个已考虑照片的和)。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<string>
#include<map>
#include<vector>
#include<queue>
#include<deque>
using namespace std;
#define ll long long
#define PI acos(-1)
#define INF 0x3f3f3f3f
#define NUM 33010
#define debug true
#define lowbit(x) ((-x)&x)
#define ffor(i,d,u) for(int i=(d);i<=(u);++i)
#define _ffor(i,u,d) for(int i=(u);i>=(d);--i)
#define mst(array,Num,Kind,Count) memset(array,Num,sizeof(Kind)*(Count))
const int p = 1e9+7;
int n;
int g[15][15]={},dp[NUM],f[NUM]={},d[15][NUM]={},sum[NUM];
map < int , int > mm;
vector < int > pos;
vector < int > num;
template <typename T>
inline void read(T &x){
    char ch = getchar();x = 0;
    for (; ch < '0' || ch > '9'; ch = getchar());
    for (; ch >='0' && ch <= '9'; ch = getchar()) x = x * 10 + ch - '0';
}
template <typename T>
inline void write(T x)
{
    int len=0;char c[21];
    if(x<0)putchar('-'),x*=(-1);
    do{++len;c[len]=(x%10)+'0';}while(x/=10);
    _ffor(i,len,1)putchar(c[i]);
}

inline void AC()
{
    char c[35];
    read(n);
    int x,y,k,u,len;
    ffor(i,1,(n*(n-1))<<1)read(x),read(y),++g[x-1][y-1];
    u=1<<n;--u;
    mst(dp,INF,int,u+1);
    ffor(i,0,n-1)mm[1<<i]=i;
    ffor(i,1,u)
    {
        f[i]=f[i-(k=lowbit(i))],x=k=mm[k];
        while(i>>(++x))f[i]+=(((i>>x)&1)?fabs(g[x][k]-g[k][x]):0);
    }
    ffor(i,0,n-1)
    {
        k=u^(1<<i);
        ffor(j,1,u)
        {
            if((j>>i)&1)continue ;
            x=lowbit(j),d[i][j]=d[i][j-x]+g[mm[x]][i];
        }
    }
    dp[0]=0;
    ffor(i,0,u)
    {
        k=u^i,len=0;
        num.clear();
        for(int j=k;j;j=(j-1)&k)num.push_back(j);
        pos.clear();
        while(k)
        {
            if(k&1)pos.push_back(len);
            k>>=1,++len;
        }
        x=(1<<pos.size())-1,sum[0]=0;
        len=num.size();
        ffor(j,1,x)
        {
            k=num[len-j];
            y=lowbit(j);
            sum[j]=sum[j-y]+d[pos[mm[y]]][i];
            dp[i|k]=min(dp[i|k],dp[i]+f[k]+sum[j]);
        }
    }
    write(dp[u]),putchar('\n');
}
int main()
{
    AC();
    return 0;
}