Kuroni and the Punishment 随机排序+质因子分解

你可以每次在给定的数组里进行这样的操作：选定一个元素令它加一或者减一，不能小于1。
求令数组里的最大公因数大于一的最少操作次数。

解法

首先一定可以通过+1，-1操作让所有的元素变为偶数，这样gcd=2，最坏情况n个奇数，这样要n次操作，每个数+1。
也就是说答案是不能超过n的。
如果我要对m个数都进行两次及以上操作，那么最好的方案中：这些数的数量不能超过n/2，因为2*n/2=n，所以m<=n/2。也就是说，至少有n/2的数，我可以通过对这些数每个数进行至多一次操作，这个数x可能会变成x-1,x,x+1。

现在假设有n个数，p个数需要修改大于等于两次，q个数需要修改最多一次，p+q=n，p<=n/2, q>=n/2, 现在我随机选择一个数x，对这个数的x-1,x,x+1进行质因子分解，然后把这些因子去跟n个数比对，看是不是有某个因子可以成为所有数的因子，求修改次数的最小值就可以了，但是这个数x因为是仅仅尝试了x-1,x,x+1,所以他应该是q个数当中的数，所以我们随机选是有概率选错的，这个概率不会超过1/2，当选多个数进行这样的尝试的时候，这个概率会变得更小，当选m个数的时候，会缩小至$\frac{1}{2^m}$，取m为100，在程序里就几乎不可能出错了。

所以随机打乱给定的数组，然后取前一百个数，尝试枚举质因子取求次数的最小值即可。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstdlib>
#include <ctime>
#include <vector>
using namespace std;
typedef long long ll;
const int N = 1e6+10;

ll p[N], prime[N], m, v[N];
void primes(ll n){
	for(ll i = 2; i <= n; i++){
		if(!p[i]) p[i] = i, prime[++m] = i;
		for(ll j = 1; j <= m; j++){
			if(prime[j] > p[i] || prime[j] > n/i) break;
			p[i*prime[j]] = prime[j];
		}
	}
//	for(int i = 1; i <= 100; i++)
//		cout << prime[i] << endl;
}

vector<ll> vec;
void insert(ll x){
	if(!x) return;
	ll t = x;
	for(int i = 1; i <= m && prime[i]<= t; i++){
		if(t%prime[i] == 0){
			if(v[i] == 0) {
				v[i] = 1;
				vec.push_back(prime[i]);
			}
			v[i] = 1;
			while(t%prime[i]==0) t /= prime[i];
		}
		if(t == 1) break;
	}
	if(t > 1) vec.push_back(t);
}
ll n, a[N];
int main(){
	srand(time(NULL));
	primes(1e6);
	scanf("%lld", &n);
	for(int i = 1; i <= n; i++){
		scanf("%lld", &a[i]);
	}
	random_shuffle(a+1, a+n+1);
	for(int i = 1; i <= 100 && i <= n; i++){
		insert(a[i]);
		insert(a[i]-1ll);
		insert(a[i]+1ll);
	}
	ll ans = 1e9;
	for(int i = 0; i < vec.size(); i++){
		ll sum = 0;
		ll y = vec[i];
		for(int j = 1; j <= n; j++){
			sum += (a[j]<y)? (y-a[j]):min(a[j]%y, y-a[j]%y);
		}
		ans = min(ans, sum);
	}
	printf("%lld\n", ans);
	return 0;
}