暑训日记Day3 T3《苹果》题解 (随机化算法)

T3 苹果

【问题描述】

$Y$有 $n$ 个苹果，他给每个苹果定义了一个味觉值，但是$Y$获得舒适度的机制非常特别，他在吃完一些苹果后只能获得等同于它们味觉值最大公约数的舒适度，现在$Y$需要吃掉恰好一半的苹果来填饱肚子，他想知道他能获得的最大舒适度是多少。

【输入格式】

第一行一个正偶整数$n$。

第二行 $n$ 个正整数 $a_i$ 表示每个苹果的味觉值。

【输出格式】

一行一个整数表示最大舒适度。

【输入样例】

6
1 2 3 4 5 6

【输出样例】

2

【数据范围】

对于50%的数据，$n\le 1000，a[i]\le 1010$

对于100%的数据，$n\le 2\ast 10^5，a[i]\le 1012$

数据有一定梯度

【题解】

这道题是最逆天的一道，正解竟然是随机化。。。

随机出一个数并枚举其因数，贪心出“这些因数中较大的，且仍然为不小于$n/2$个元素的因数，的那个因数”，来更新最终的答案。

由于所选的被枚举因数的元素，只有一般的概率为最终所选的元素，所以我们多随机几次，就有很大的概率得到正确答案。

#include<bits/stdc++.h>
#define REG register
#define LL long long 
using namespace std;
const int maxn = 200005;
LL d[maxn],a[maxn],f[maxn],s,T = 10;
int n;
inline LL Read(){
  char c = getchar();
  LL x = 0,f = 1;
  while(c < '0' || c > '9'){
    if(c == '-') f = -1;
    c = getchar();
  }
  while(c >= '0' && c <= '9'){
    x = x * 10 + c - '0';
    c = getchar();
  }
  return x * f;
}
int main(){
  freopen("apple.in","r",stdin);
  freopen("apple.out","w",stdout);
  scanf("%d",&n);
  for(int i = 1;i <= n;++i) a[i] = Read();
  srand(time(NULL));  
  while(T--){
	LL g = a[(rand()*32768 + rand()) % n + 1];
	//注意，要保证2*10^5全部被取到 
	LL t = 0;  	
	for(REG int i = 1;i * i <= g;i++){		
	  if(g % i == 0){
		d[++t] = i;  	
	    if(i * i != g)
		  d[++t] = g / i;
	  }
	}
	sort(d + 1,d + t + 1);//贪心 
	for(REG int i = 1;i <= t;i++) f[i] = 0;
	for(REG int i = n;i;i--)
	  f[lower_bound(d + 1,d + t + 1,
	                __gcd(a[i],g)) - d]++;
	for(REG int i = 1;i <= t;i++){
	  for(REG int j = i + 1;j <= t;j++)
	    if(d[j] % d[i] == 0)
		  f[i] += f[j];//扫一遍，防止遗漏 
	  if(f[i] * 2 >= n)
	    s = max(s,d[i]);
	}
  }
  printf("%lld",s);
  return 0;
}