素数伴侣--最大二分匹配

#include<bits/stdc++.h>
using namespace std;
#define N 100
int edge

,cx
,cy
;//edge记录两点的关系，如果两点相连，则edge【i】【j】为1
int visited
;//判断该店是否被访问过
int nx,ny,res;
const int M=60000+100;
bool prim[M];
void initprim()
{
memset(prim,true,sizeof(prim));
prim[0]=prim[1]=false;
for(int i=2; i*i<=M; i++)
for(int j=i*i; j<M; j+=i)
{
if(prim[j])prim[j]=false;
}
}
int path(int u)
{
int v;
for(v=0; v<ny; v++)
{
if(edge[u][v]&&!visited[v])
{
visited[v]=1;
if(cy[v]==-1||path(cy[v]))////如果y集合中的v元素没有匹配或者是v已经匹配，但是从cy[v]中能够找到一条增广路
{
cx[u]=v;
cy[v]=u;
return 1;
}
}
}
return 0;
}
int main()
{
int n;
initprim();
while(scanf("%d",&n)!=EOF)
{
int i,j,a[100]= {0},a1[100]= {0},a2[100]= {0};
nx=0,ny=0,res=0;
memset(cx,0xff,sizeof(cx));
memset(cy,0xff,sizeof(cy));//初始值为-1表示两个集合中都没有匹配的元素！
memset(edge,0,sizeof(edge));
for(i=0; i<n; i++)
{
scanf("%d",&a[i]);
if(a[i]%2==1)
a1[nx++]=a[i];
else
a2[ny++]=a[i];
}
for(i=0; i<nx; i++) //初始化edge数组，如果两个数之和为素数，则edge[i][j]置1
{
for(j=0; j<ny; j++)
{
if(prim[a1[i]+a2[j]])edge[i][j]=1;
}
}
for(i=0; i<nx; i++)
{
if(cx[i]==-1)
{
memset(visited,0,sizeof(visited));
res+=path(i);
}
}
printf("%d\n",res);
}
}