XJOI网上同步训练DAY5 T1

 思路:考虑得出，最终的集合一定是gcd=1的集合，那么我们枚举n个数中哪个数必须选，然后把它质因数分解，由于质数不会超过9个，可以状态压缩，去得出状态为0的dp值就是答案。

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<iostream>
int n,l[305],c[305],p[3005],val[30005],f[305][1025];
int read(){
    int t=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while ('0'<=ch&&ch<='9'){t=t*10+ch-'0';ch=getchar();}
    return t*f;
}
int main(){
    n=read();
    int ans=1<<30;
    for (int i=1;i<=n;i++)
     l[i]=read();
    for (int i=1;i<=n;i++)
     c[i]=read();
    for (int i=1;i<=n;i++){
        int tot=0,t=l[i];
        for (int j=2;j*j<=t;j++){
            if (t%j==0){
                tot++;p[tot]=j;
                while (t%j==0) t/=j;
            }
        }
        if (t!=1) p[++tot]=t;
        for (int j=1;j<=n;j++)
         if (j!=i){
            val[j]=0;
            for (int k=1;k<=tot;k++)
             if (l[j]%p[k]==0)
              val[j]|=1<<(k-1);
         }
        int sum;
        val[i]=sum=(1<<tot)-1;
        for (int j=1;j<=n;j++)
         for (int k=0;k<=sum;k++)
          f[j][k]=1<<30;
        f[i][sum]=c[i];
        for (int j=i;j<=n-1;j++)
         for (int k=0;k<=sum;k++)
          if (f[j][k]!=1<<30){
            f[j+1][k]=std::min(f[j+1][k],f[j][k]);
            f[j+1][k&val[j+1]]=std::min(f[j+1][k&val[j+1]],f[j][k]+c[j+1]);
          }
        if (f[n][0]!=1<<30){
            ans=std::min(ans,f[n][0]);
        }    
    }
    if (ans!=1<<30) printf("%d
",ans);
    else 
    printf("-1
");
    return 0;
}