codeforces1305F 随机

题目传送门

题意：

给你  个数  ，每次操作可以使一个数增加  或 减少  。操作后需要保证这个数是正数。

问你最少操作多少次使这些数的  大于  。

数据范围：  。

题解：

我们让每个数都变成  的倍数，上界是  。

如果有超过一半的数变化超过  ，那么就超过上界了。

所以最多有一半的数的变化不能超过  。

我们在序列中随机一个数  ，然后枚举  的质因子作为整个序列的  。

因为每次随机失败的概率是  ，那么  次随机失败的概率大概是  。 

感受：

常数写的有点大。

看题解之前能想到必须得知道gcd是什么才能做，但接下来推不下去了。

其实这种题也只能是随机了，但是不看题解确实想不到。

代码：

#include<bits/stdc++.h>
using namespace std ;
typedef long long ll ;
const int maxn = 2e5 + 5 ;
const int maxm = 1e6 + 5 ;
const int inf = 0x3f3f3f3f ;
mt19937  rnd(chrono::high_resolution_clock::now().time_since_epoch().count()) ;
int n , m ;
ll a[maxn] ;
int cnt = 0 ;
bool vis[maxm] ;
int prime[maxm] ;
void get_prime()
{
   memset(vis , 0 , sizeof(vis)) ;
   vis[1] = 1 ;
   for(int i = 2 ; i <= maxm - 5 ; i ++)
   {
     if(!vis[i]) 
     prime[++ cnt] = i ;
     for(int j = 1 ; j <= cnt && i * prime[j] <= maxm - 5 ; j ++)
     {
       vis[i * prime[j]] = 1 ;
       if(i % prime[j] == 0) break ;
     }
   }
}
ll work(ll x)
{
	ll num = 0 ;
	for(int i = 1 ; i <= n ; i ++)
	{
		ll y = a[i] % x ;
		if(a[i] < x)  num += x - y ;
		else  num += min(y , x - y) ;
	}
	return num ;
}
ll cal(ll x)
{
	ll ans = 1e18 ;
	if(x <= 1)  return ans ;
	for(int i = 1 ; i <= cnt && prime[i] <= x ; i ++)
	{
		if(x % prime[i] != 0)  continue ;  
	    ans = min(ans , work(ll(prime[i]))) ;
	    while(x % prime[i] == 0)  x /= prime[i] ;
	}
	if(x > 1)  ans = min(ans , work(x)) ; 
	return ans ;
}
ll solve(ll x)
{
   ll ans = 1e18 ;
   ans = min(ans , cal(a[x] - 1)) ;
   ans = min(ans , cal(a[x])) ;
   ans = min(ans , cal(a[x] + 1)) ;
   return ans ;	
}
int main()
{
	ll ans = 0 ;
	scanf("%d" , &n) ;
	for(int i = 1 ; i <= n ; i ++)  scanf("%lld" , &a[i]) ;
	get_prime() ;
	for(int i = 1 ; i <= n ; i ++)
	  if(a[i] % 2 == 1)  ans ++ ;
	for(int i = 1 ; i <= 30 ; i ++)
	{
		int x = rnd() % n + 1 ;
		ans = min(ans , solve(x)) ;
	}
	printf("%lld\n" , ans) ;
	return 0 ;
}