牛客CCPC.2017秦皇岛站-重现赛，H（二分图匹配）

一开始确实想到了用匈牙利算法去做，但是错在了两个地方：
1.无向图应建立双边；
2.想到了用匈牙利算法求出所有的元素不重复的prime set后，设当前总集合数为Set，k-Set>0的情况想得过复杂了而导致错误。
另外，这题本人觉得有一点玄学的地方就是，我原本是用前向星建图的，但是超时了，后来看别人的代码改用了vector建临接矩阵才过的。。。（感觉前向星应该是比vector快挺多的才对，之前有遇到过dijkstra用vector建图不过而用前向星稳稳地过了的情况。。。STL按道理不是应该比较慢的吗？）

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<vector>
#include<set>
using namespace std;
const int maxn=3e3+5;
const int maxm=2e6+5;
int a[maxn],match[maxn],n;
bool vis[maxn],num[maxm];

void getprime()
{
    memset(num,1,sizeof(num));
    num[0]=num[1]=0;
    for(int i=2;i<maxm;++i)
    {
        if(num[i])
        {
            for(int j=i+i;j<maxm;j+=i)
                num[j]=0;
        }
    }
}

bool find(int x)		//匈牙利算法模板
{
    vis[x]=true;		//注意，因为本题的图是无向图，所以vis[x]也要设为1（可以自己思考一下），有向图就可以不用
    int len=vt[x].size();
    for(int i=0; i<len ;++i)
    {
        int v=vt[x][i];
        if(!vis[v])
        {
            vis[v]=true;
            if(match[v]==0||find(match[v]))
            {
                match[v]=x; match[x]=v;
                return true;
            }
        }
    }
    return false;
}
 
int main()
{
    int T;
    getprime();
    scanf("%d",&T);
    while(T--)
    {
        int k;
        scanf("%d%d",&n,&k);
        memset(match,-1,sizeof(match));
		cnt=2;
        for(int i=1;i<=n;++i)
        {
            scanf("%d",&a[i]);vt[i].clear();
            for(int j=1;j<i;++j)
            	if(num[a[i]+a[j]])
            		vt[i].push_back(j),vt[j].push_back(i),match[i]=match[j]=0;	//这里非常巧妙和充分的利用了match数组,
            		//match[i]=0说明a[i]可以跟某个数形成prime set,那所以的match[i]=0就组成了所以prime set的广义并集合
        }
        int ans1=0,ans2=0;
        for(int i=1;i<=n;++i)
        {
            if(match[i]==0)
           	{
           		for(int j=1;j<=n;++j)    vis[j]=0;
           		 if(find(i))    ++ans1;
			}
        }
        for(int i=1;i<=n;++i)	if(match[i]==0)	++ans2;
        if(ans1>k)	printf("%d\n",2*k);
        else printf("%d\n",2*ans1+min(k-ans1,ans2));
    }
    return 0;
}