Codeforces Round #259 (Div. 2) D. Little Pony and Harmony Chest （状压dp || 位运算）

题意：
给你N个数（a[1],a[2],….a[n]），让你给出N个两两互质的数，使得sum(abs(a[i]-b[i]))最小。
题解：
这道题最最重要的是要想到，求一个数和一个序列其他数是否互质可以用质因子是否重合来做，维护当前序列的质因子集合，然后要加入一个新的元素时，只需要检查该元素质因子是否与集合中质因子重合，就代表是否互质。
然后检查质因子重合，还可以用位运算来加快速度，因为最多只会用到17个质因子，一共(1<<17)。最后再加上状压dp来算出最优解。

代码

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define fir first
#define se second
#define pb(a) push_back(a)
const int maxn = 1<<20;
const int inf = 0x3f3f3f3f;
const double eps = 0.000000001;
typedef pair<LL,LL> pll;
typedef pair<int,int> pii;
typedef pair<double,double> pdd;

int a[150],pp[60],p[20] = {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59};
int dp[150][1<<17],pa[150][1<<17],pr[150][1<<17];

void dfs(int len,int pos) {
    if(len == 0) return ;
    dfs(len-1,pr[len][pos]);
    printf("%d ",pa[len][pos]);
}

int main() {
    int n;
    cin>>n;
    for(int i = 1;i <= n;i++) cin>>a[i];
    for(int i  = 1;i <= 60;i++) {   
        pp[i] = 0;
        for(int j = 0;j < 17;j++) {
            if(i%p[j] == 0) pp[i] |= (1<<j);
        }
    }

    memset(dp,0x3f,sizeof(dp));
    dp[0][0] = 0;
    for(int i = 1;i <= n;i++) {
        for(int j = 1;j <= 60;j++) {
            for(int k = 0;k < (1<<17);k++) {
                if(k&pp[j] || dp[i-1][k] == inf) continue;
                int m = dp[i-1][k] + abs(a[i]-j);
                if(m < dp[i][k|pp[j]]) {
                    //if(i == 1) cout<<j<<" "<<m<<"  ffff\n";
                    dp[i][k|pp[j]] = m;
                    pa[i][k|pp[j]] = j;
                    pr[i][k|pp[j]] = k;
                }
            }
        }
    }
    //cout<<"fff"<<endl;
    int ans = inf,pos;
    for(int i = 0;i < (1<<17);i++) {
        if(ans > dp[n][i]) {
            pos = i;
            ans = dp[n][i];
        }
    }
    dfs(n,pos);
}

感觉自己要蠢死了。。。。