LightOJ 1356 Prime Independence（质因数分解+最大独立集+Hopcroft-Carp）

  1 #include<iostream>
  2 #include<algorithm>
  3 #include<cstring>
  4 #include<cstdio>
  5 #include<vector>
  6 #include<queue>
  7 #include<cmath>
  8 #include<map>
  9 #include<stack>
 10 using namespace std;
 11 
 12 const int maxn=500000+5;
 13 const int INF=0x3f3f3f3f;
 14 
 15 int n;
 16 int cnt=0;
 17 int primes[maxn];
 18 int vis[maxn];
 19 
 20 vector<int>g[40005];
 21 int um[maxn],vm[maxn];
 22 int dx[maxn],dy[maxn],dis;
 23 
 24 void get_primes()
 25 {
 26     int m=sqrt(maxn+0.5);
 27     for(int i=2;i<=m;i++)
 28     {
 29         if(!vis[i])
 30         {
 31             for(int j=i*i;j<=maxn;j+=i)
 32                 vis[j]=1;
 33         }
 34     }
 35     for(int i=2;i<=maxn;i++)
 36         if(!vis[i])   primes[cnt++]=i;
 37 }
 38 
 39 bool searchP()
 40 {
 41     queue<int>q;
 42     dis=INF;
 43     memset(dx,-1,sizeof(dx));
 44     memset(dy,-1,sizeof(dy));
 45     for(int i=1;i<=n;i++)
 46         if(um[i]==-1)
 47         {
 48             q.push(i);
 49             dx[i]=0;
 50         }
 51     while(!q.empty())
 52     {
 53         int u=q.front();q.pop();
 54         if(dx[u]>dis)  break;
 55         for(int i=0;i<g[u].size();i++)
 56         {
 57             int v = g[u][i];
 58             if(dy[v]==-1)
 59             {
 60                 dy[v]=dx[u]+1;
 61                 if(vm[v]==-1)  dis=dy[v];
 62                 else
 63                 {
 64                     dx[vm[v]]=dy[v]+1;
 65                     q.push(vm[v]);
 66                 }
 67             }
 68         }
 69     }
 70     return dis!=INF;
 71 }
 72 
 73 bool dfs(int u)
 74 {
 75     for(int i=0;i<g[u].size();i++)
 76     {
 77         int v = g[u][i];
 78         if(!vis[v]&&dy[v]==dx[u]+1)
 79         {
 80             vis[v]=1;
 81             if(vm[v]!=-1&&dy[v]==dis) continue;
 82             if(vm[v]==-1||dfs(vm[v]))
 83             {
 84                 vm[v]=u;um[u]=v;
 85                 return 1;
 86             }
 87         }
 88     }
 89     return 0;
 90 }
 91 
 92 int maxMatch()
 93 {
 94     int res=0;
 95     memset(um,-1,sizeof(um));
 96     memset(vm,-1,sizeof(vm));
 97     while(searchP())
 98     {
 99         memset(vis,0,sizeof(vis));
100         for(int i=1;i<=n;i++)
101           if(um[i]==-1&&dfs(i))  res++;
102     }
103     return res;
104 }
105 
106 int a[maxn],b[maxn],p[maxn],num[maxn];
107 
108 int main()
109 {
110     //freopen("D:\\input.txt","r",stdin);
111     int T;
112     int kase=0;
113     get_primes();
114     scanf("%d",&T);
115     while(T--)
116     {
117         scanf("%d",&n);
118         for(int i=0;i<=n;i++)   g[i].clear();
119         memset(b,-1,sizeof(b));
120         for(int i=1;i<=n;i++)
121             scanf("%d",&a[i]);
122         sort(a+1,a+1+n);
123         for(int i=1;i<=n;i++)
124             b[a[i]]=i;
125 
126         memset(num,0,sizeof(num));
127         for(int i=1;i<=n;i++)
128         {
129             int x=a[i];
130             int sum1=0;
131             int sum2=0;
132             for(int j=0;primes[j]*primes[j]<=x;j++)
133             {
134                 if(x%primes[j]==0)
135                 {
136                     p[sum1++]=primes[j];
137                     while(x%primes[j]==0)   {x/=primes[j];sum2++;}
138                 }
139             }
140             if(x>1)  {p[sum1++]=x;sum2++;}
141             num[i]=sum2;
142             for(int j=0;j<sum1;j++)
143             {
144                 int t=b[a[i]/p[j]];
145                 if(t<i&&t>0)
146                 {
147                     if((sum2&1)==(num[t]&1))  continue;
148                     if(sum2&1)  g[i].push_back(t);
149                     else  g[t].push_back(i);
150                 }
151             }
152         }
153 
154         printf("Case %d: %d\n",++kase,n-maxMatch());
155     }
156 
157     return 0;
158 }