「国庆训练」Kingdom of Obsession（HDU-5943）

题意

给定$s,n$，把$s+1,s+2,...,s+n$这$n$个数填到$1,2,...,n$里，要求$x$只能填到$x$的因子的位置（即题目中$x\%y=0$那么x才能放在y位置的要求），问是否能够把所有数都填上去。

分析

建模之后就是条二分图的板子题。
注意到一个有趣的事实：如果$[s+1,s+n]$中有两个质数，那他们肯定完蛋。因为它们一定都会和1连边然后翻车。
这意味着什么呢？在$10^9$的范围里面，质数最大的间隔是282（我抄网上的结论的，正确不保证，反正不多），那么我们就可以认定如果n比这个数更大一定是No。
然后就可以快乐建图了，把取模为0的建一条有向边，看看最后能否找到完美匹配即可。
实现上，s应当比n大。为什么？如果$s<n$,那么$[s+1,n]$直接放在自己的位置上即可，因为不可能放在比他们更大的地方了（集合的交）：那爽死了，我们只需要考虑$[1,s]$与剩下$[n+1,s+n]$的匹配即可，这个东西的实质就是n、s交换。有没有其他数$k$放在这些位置上面（如一个数$x$）然后这些数字没法放到其他数字的反例呢？没有哦，因为如果$k$能放在$[s+1,n]$的对应位置，那么作为这个数本身的$x$一定能放在$k$上。所以得到结论：可以直接将$s$和$n$调换来缩小$n$。

代码

#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <stack>
#define MP make_pair
#define PB push_back
#define fi first
#define se second
#define ZERO(x) memset((x), 0, sizeof(x))
#define ALL(x) (x).begin(),(x).end()
#define rep(i, a, b) for (int i = (a); i <= (b); ++i)
#define per(i, a, b) for (int i = (a); i >= (b); --i)
#define QUICKIO                  \
    ios::sync_with_stdio(false); \
    cin.tie(0);                  \
    cout.tie(0);
#define MS(x,y) memset(x,y,sizeof(x))
using namespace std;
typedef long long ll;

const int MAXN=2018;// actually, max distance: 282
bool mat[MAXN][MAXN];
bool is_prime(int x)
{
	if(x<=1) return false;
	else for(ll k=2;k*k<=ll(x);++k)
		 {
			 if(x%k==0) return false;
		 }
	return true;
}

int linker[MAXN];
bool vis[MAXN];

int n;
bool dfs(int x)
{
	rep(i,1,n) if(mat[x][i] && !vis[i])
	{
		vis[i]=true;
                if(linker[i]==-1 || dfs(linker[i]))
		{
			linker[i]=x;
			return true;
		}
	}
	return false;
}

int hungary()
{
	MS(linker,-1);
	int ans=0;
	rep(i,1,n)
	{
		ZERO(vis); // 注意这一行，应当是每一次清空
	 	if(dfs(i)) ans++;
	}
	return ans;
}

int main()
{
	int T; scanf("%d",&T);
	rep(kase, 1, T)
	{
		int primecnt=0;
		int s;
		ZERO(mat);
		scanf("%d%d", &n, &s);
		if(s<n) swap(s,n); // s should be bigger
		if(n>1000)
		{
			printf("Case #%d: No\n", kase);
			continue;
		}
		rep(i,1,n)
		{
			if(is_prime(s+i)) primecnt++;
			if(primecnt>1) break;
		}
		if(primecnt<=1)
		{
			rep(i,1,n)
			{
				int t=s+i;
				rep(j,1,n) if(t%j==0)
					mat[i][j]=true;
			}
			if(hungary()!=n)
				printf("Case #%d: No\n", kase);
			else printf("Case #%d: Yes\n", kase);
		}
		else printf("Case #%d: No\n", kase);
	}
	return 0;
}