Description
有一个无限长的序列,一些位置被翻转了,你每次可以翻转连续奇质数个位置,问最少多少次操作可以将所有位置复原。
Solution
(这是LSTete的题解)
定义序列 S S 表示元件的方向, 0下1上, 特别地, . 不妨考虑差分, 定义序列 C C , 当时, Ci=0 C i = 0 , 否则 Ci=1 C i = 1 , 不难证明序列 C C 中 1 1 的个数一定为偶数个.
这样, 一次操作 [ l, r ] [ l , r ] 只会改变 Cl C l 和 Cr+1 C r + 1 的值, 那么现在对一段区间操作变成了对两点操作, 问题变成了如何两两配对使得操作数最小.
考虑任意两个数 i, j i , j 之间需要的操作数.
①当 |i−j| | i − j | 为奇质数时, 需要一次操作.
②当 |i−j| | i − j | 为偶数时, 需要两次操作. (哥德巴赫猜想)
③当 |i−j| | i − j | 为奇非质数时, 需要三次操作. (分解为偶数减奇质数之差)
希望尽量都是①情况, 其次尽量剩余②情况, 实在不行补一次③情况.
问题变得简单了, 把序列 C C 中出现的位置按奇偶变成二分图, 差为奇质数则连边, 跑最大匹配. 对于残量图, 同组(同奇偶)则需要两次操作(情况②), 若还剩余一对异组则需要三次操作(情况③).
Code
/************************************************
* Au: Hany01
* Date: Jul 16th, 2018
* Prob: ARC080 F
* Email: hany01@foxmail.com
* Inst: Yali High School
************************************************/
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef pair<int, int> PII;
#define File(a) freopen(a".in", "r", stdin), freopen(a".out", "w", stdout)
#define rep(i, j) for (register int i = 0, i##_end_ = (j); i < i##_end_; ++ i)
#define For(i, j, k) for (register int i = (j), i##_end_ = (k); i <= i##_end_; ++ i)
#define Fordown(i, j, k) for (register int i = (j), i##_end_ = (k); i >= i##_end_; -- i)
#define Set(a, b) memset(a, b, sizeof(a))
#define Cpy(a, b) memcpy(a, b, sizeof(a))
#define x first
#define y second
#define pb(a) push_back(a)
#define mp(a, b) make_pair(a, b)
#define ALL(a) (a).begin(), (a).end()
#define SZ(a) ((int)(a).size())
#define INF (0x3f3f3f3f)
#define INF1 (2139062143)
#define debug(...) fprintf(stderr, __VA_ARGS__)
#define y1 wozenmezhemecaia
template <typename T> inline bool chkmax(T &a, T b) { return a < b ? a = b, 1 : 0; }
template <typename T> inline bool chkmin(T &a, T b) { return b < a ? a = b, 1 : 0; }
inline int read() {
static int _, __; static char c_;
for (_ = 0, __ = 1, c_ = getchar(); c_ < '0' || c_ > '9'; c_ = getchar()) if (c_ == '-') __ = -1;
for ( ; c_ >= '0' && c_ <= '9'; c_ = getchar()) _ = (_ << 1) + (_ << 3) + (c_ ^ 48);
return _ * __;
}
const int maxn = 105, maxm = 1e7 + 9, M = 1e7 + 5, maxe = maxn * maxn;
int n, mk[maxm], vis[maxn << 1], beg[maxn << 1], v[maxe], nex[maxe], pr[maxm], np[maxm], prs, tot1, tot2, S[maxn << 1], T[maxn << 1], e, cnt, mat[maxn];
inline void init(int N)
{
np[1] = 1;
For(i, 2, N) {
if (!np[i]) pr[++ prs] = i;
for (register int j = 1; j <= prs && pr[j] * i <= N; ++ j) {
np[pr[j] * i] = 1;
if (!(i % pr[j])) break;
}
}
}
inline void add(int uu, int vv) { v[++ e] = vv, nex[e] = beg[uu], beg[uu] = e; }
int dfs(int u) {
for (register int i = beg[u]; i; i = nex[i])
if (!vis[v[i]]) {
vis[v[i]] = 1;
if (!mat[v[i]] || dfs(mat[v[i]])) { mat[v[i]] = u; return 1; }
}
return 0;
}
int main()
{
#ifdef hany01
File("arc080f");
#endif
n = read(), init(M);
For(i, 1, n) mk[read()] = 1;
For(i, 1, M) if (mk[i] != mk[i - 1])
(i & 1) ? (S[++ tot1] = i) : (T[++ tot2] = i);
For(i, 1, tot1) For(j, 1, tot2)
if (!np[abs(S[i] - T[j])]) add(i, j);
For(i, 1, tot1) {
Set(vis, 0);
if (dfs(i)) ++ cnt;
}
printf("%d\n", cnt + ((tot1 - cnt) >> 1) * 2 + ((tot2 - cnt) >> 1) * 2 + (((tot1 - cnt) & 1)) * 3);
return 0;
}
//槛菊愁烟兰泣露。罗幕轻寒,燕子双飞去。
// -- 晏殊《蝶恋花·槛菊愁烟兰泣露》