题目描述
状压dp
显然我们只要能选出一些数,使它们能够组合成1即合法。
根据裴蜀定理,显然是最大公约数为1时合法。
容易想到状态压缩dp保存每个质因数是否出现,但质因数太多了,怎么办?
可以枚举一张卡片必须被我们选,那么就只剩下最多9个质因数了。
于是做n次状压dp即可。
#include<cstdio>
#include<algorithm>
#include<cmath>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
const int maxn=300+10;
int f[maxn][2000],a[maxn],c[maxn],v[maxn],pri[maxn];
int i,j,k,l,t,n,m,tot,ans;
int main(){
scanf("%d",&n);
fo(i,1,n) scanf("%d",&a[i]);
fo(i,1,n) scanf("%d",&c[i]);
ans=1<<30;
fo(i,1,n){
k=a[i];
tot=0;
fo(j,2,floor(sqrt(k)))
if (k%j==0){
pri[++tot]=j;
while (k%j==0) k/=j;
}
if (k>1) pri[++tot]=k;
fo(j,1,n){
v[j]=0;
fo(k,1,tot)
if (a[j]%pri[k]==0) v[j]+=(1<<(k-1));
}
f[0][(1<<tot)-1]=c[i];
fo(j,0,(1<<tot)-2) f[0][j]=1<<30;
fo(j,0,n-1){
fo(k,0,(1<<tot)-1) f[j+1][k]=f[j][k];
fo(k,0,(1<<tot)-1)
if (f[j][k]!=1<<30)
f[j+1][k&v[j+1]]=min(f[j+1][k&v[j+1]],f[j][k]+c[j+1]);
}
ans=min(ans,f[n][0]);
}
if (ans==1<<30) printf("-1\n");else printf("%d\n",ans);
}