【HDU 5943】素数间隔+二分匹配

题意：有n个人，每个人的标号以此是s+1到s+n，要求将所有人重新排序之后满足每个人的位置y能够保证被他的标号整除，就是数组重排之后满足每一位的a[i]%i==0

题解：看到数据样本其实发现非常大，一般数据结构都开不下。但是初步分析我们发现如果该区间内存在两个以上的素数那么肯定无法匹配，所以事先判断一遍，避免开大数组的情况（不避开会RE），而一个只含一个素数的区间最大1000+可以存下。

另外，新区间和1-n可能有重复的部分要删去，因为这部分一定满足所以不需要匹配，所以最后代码中判断是两端区间分别为[1,min(n,s)]和[max(n,s),n+s]。

碎碎念：这题修修改改了好几题，RE,TE,WA尝了个遍，怕是失心疯了.....

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<cmath>
#include<algorithm>
#include<set>
#define maxn 1150
using namespace std;
typedef long long ll;
ll n,s;
ll used[maxn];
int vis[maxn];
vector<ll>G[maxn];

bool match(ll u){
	ll v;
	for(ll j=0,sz=G[u].size();j<sz;j++){
		v=G[u][j];
		if(vis[v]) continue;
		vis[v]=1;
		if(!used[v]||match(used[v])){
			used[v]=u;
			return true;
		}
	}
	
	return false;
}

int maxMatch(){
	for(int i=1;i<=n+n;i++) used[i]=0;
	int ans=0;
	for(ll i=1;i<=n;i++){
		memset(vis,0,sizeof(vis));
		if(match(i)) ans++;
	}
	return ans;
}

bool isPrime(ll a){
	if(a==0||a==1)return false;
    if(a==2)return true;
    if(a%2==0)return false;
    for(ll i=3;i<=sqrt(a);i+=2)
        if(a%i==0)return false;
    return true;
}

void init(){
	for(int i=0;i<=n+n;i++) G[i].clear();
}

int main(){
	int T,cas=1;
	scanf("%d",&T);
	while(T--){
		scanf("%lld%lld",&n,&s);
		if(n>s) swap(n,s);
		int flag=2;
		for(ll i=1;i<=n;i++){
			ll tmp=i+s;
			if(isPrime(tmp)){
				flag--;
				if(flag==0) break;
			}
		}
		if(flag){
			init();
			for(ll i=1;i<=n;i++){
				ll tmp=i+s;
				for(ll j=1;j<=n;j++){
					if(tmp%j==0){
						G[i].push_back(j+n);
						G[j+n].push_back(i);
					}
				}
			}
			int tmp=maxMatch();
	//		printf("here tmp=%d\n",tmp);
			if(tmp!=n) flag=0;
		}
		
		printf("Case #%d: ",cas++);
		if(flag) printf("Yes\n");
		else printf("No\n");
	
	}	
	return 0;
}