hihocoder1158 质数相关

题意：

两个数$a$和 $b$ ($a<b$)被称为质数相关，是指 $a × p = b$，这里$p$是一个质数。一个集合$S$被称为质数相关，是指$S$中存在两个质数相关的数，否则称$S$为质数无关。如$\{2, 8, 17\}$质数无关，但$\{2, 8, 16\}$, $\{3, 6\}$质数相关。现在给定一个集合$S$，问$S$的所有质数无关子集中，最大的子集的大小。

数据范围：$|S| \le 1000$

思路：

我们把质数相关想成两个点之间有一条无向边，那么问题就是转化为求点的子集，子集内所有点之间都没有边相连。

很裸的求二分图的最大独立集的问题，一个定理：二分图的最大匹配数 = 最小边覆盖点集个数，最大独立集 = $N$ - 最大匹配数。（浅显的证明）

但是，上面集合构成的图一定是二分图吗？当然是二分图，最直接的想法，因为$a \cdot p_1$和 $a \cdot p_2$（$p_1，p_2为质数$）不是质数相关的！

代码：

#include <bits/stdc++.h>
using namespace std;

const int maxn=500000+5;

bool vis[maxn];
int prime[maxn],num;

void init()
{
    static int i,j;
    memset(vis,0,sizeof(vis));
    num=0;
    vis[1]=vis[0]=1;
    for(i=2;i<maxn;i++)
    {
        if(!vis[i])prime[num++]=i;
        for(j=0;j<num && i*prime[j]<maxn;j++){
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)break;
        }
    }
}

int n;
int f[1005];
bool p[1005][1005],p1[1005];
int ans[2];

void dfs(int u,bool pd)
{
    //printf("dfs %d\n",u);
    p1[u]=true;
    ans[pd]++;
    for(int v=1;v<=n;v++)
    {
        if(p1[v])continue;
        if(p[u][v])dfs(v,!pd);
    }
}

int main()
{
    init();
    static int T,Cas;
    scanf("%d",&T);
    for(Cas=1;Cas<=T;Cas++){
        scanf("%d",&n);
        static int i,j;
        for(i=1;i<=n;i++){
            p1[i]=false;
            scanf("%d",&f[i]);
        }
        sort(f+1,f+n+1);
        for(i=1;i<=n;i++){
            p[i][i]=0;
            for(j=i+1;j<=n;j++)
                if(f[j]%f[i]==0 && !vis[f[j]/f[i]])p[i][j]=p[j][i]=1;
                else p[i][j]=p[j][i]=0;
        }
        int Ans=0;
        for(i=1;i<=n;i++){
            if(p1[i])continue;
            ans[0]=ans[1]=0;
            dfs(i,0);
            Ans+=max(ans[0],ans[1]);
        }
        printf("Case #%d: %d\n",Cas,Ans);
    }
    return 0;
}