题解思路:因为题目说有一个数除以另外一个数如果等于一个素数的话就只能一个加入集合,那么这两个数矛盾并且可以构成一个二分图,就是求二分图的最大独立子集,然后如何快速的求这些数相互矛盾呢?
把这些数拆分成看有多少个不同的素数并且除以素数,然后看下有没有在出现过,如果有就两个位置上加入一条无向边即可。最后就是跑一遍HK就可以了。
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
#include<vector>
#include<queue>
using namespace std;
const int mx = 500000;
const int inf = 0x3f3f3f3f;
int vis[mx+5];
int a[mx/10];
int p[mx];
int num[mx/10];
int n;
int dx[mx/10],dy[mx/10];
int dis;
int x[mx/10],y[mx/10];
bool mark[mx/10];
struct node{
int next,v;
}E[mx];
int head[mx/10];
int tot;
void init(){
memset(vis,0,sizeof(vis));
for(int i = 2; i <= mx; i++){
if(!vis[i]) p[++p[0]] = i;
for(int j = 1; j <= p[0]&&p[j]*i <= mx; j++){
vis[p[j]*i] = 1;
if(i%p[j]==0)
break;
}
}
}
void add(int u,int v){
tot++;
E[tot].v = v;
E[tot].next = head[u];
head[u] = tot;
}
bool search(){
memset(dx,0,sizeof(dx));
memset(dy,0,sizeof(dy));
dis = inf;
queue<int>q;
for(int u = 1; u <= n; u++)
if(!x[u]){
dx[u] = 1;
q.push(u);
}
while(!q.empty()){
int u = q.front();
if(dx[u]>dis) break;
q.pop();
for(int i = head[u]; i; i = E[i].next){
int v = E[i].v;
if(!dy[v]){
dy[v] = dx[u]+1;
if(!y[v]) dis = dy[v];
else{
dx[y[v]] = dy[v]+1;
q.push(y[v]);
}
}
}
}
return dis!=inf;
}
bool find(int u){
for(int i = head[u]; i; i = E[i].next){
int v = E[i].v;
if(!mark[v]&&dy[v] == dx[u]+1){
mark[v] = 1;
if(y[v]&&dis == dy[v]) continue;
if(!y[v]||find(y[v])){
x[u] = v;
y[v] = u;
return true;
}
}
}
return false;
}
int MaxMatch(){
int ans = 0;
memset(x,0,sizeof(x));
memset(y,0,sizeof(y));
while(search()){
memset(mark,0,sizeof(mark));
for(int u = 1; u <= n; u++)
if(!x[u]&&find(u))
ans++;
}
// cout<<ans<<endl;
return ans/2;
}
int main(){
int t;
init();
scanf("%d",&t);
for(int casei = 1; casei <= t; casei++){
scanf("%d",&n);
for(int i = 1; i <= n; i++)
scanf("%d",&a[i]);
tot = 0;
memset(head,0,sizeof(head));
memset(vis,0,sizeof(vis));
for(int i = 1; i <= n; i++)
vis[a[i]] = i; //表示这个数的位置
//cout<<n<<endl;
for(int i = 1; i <= n; i++){
int sum = 0;
int x = a[i];
for(int j = 1; j <= p[0]&&p[j]*p[j]<=x; j++){
if(x%p[j]==0)
num[++sum] = p[j];
while(x%p[j]==0)
x/=p[j];
}
if(x>1)
num[++sum] = x;
for(int j = 1; j <= sum; j++){
x = a[i]/num[j];
if(vis[x]){
add(vis[x],i);
add(i,vis[x]);
}
}
}
// cout<<n<<endl;
printf("Case %d: %d\n",casei,n-MaxMatch());
}
return 0;
}