题意:
给你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);
}
感觉自己要蠢死了。。。。