正題

luogu
CF364D

---

題目大意

給你一個(gè)大小為n的集合，選擇一個(gè)至少為一半的子集，另其gcd最大

---

解題思路

由于數(shù)字個(gè)數(shù)很多，考慮隨機(jī)

隨機(jī)選10個(gè)數(shù)，對(duì)于每個(gè)數(shù)，先處理出約數(shù)，然后求出所有數(shù)和當(dāng)前數(shù)的gcd，這些gcd的約數(shù)權(quán)值全部加1（權(quán)值意義為gcd為這個(gè)點(diǎn)，集合最多選多少個(gè)點(diǎn)），最后找到最大滿足權(quán)值大于一半數(shù)即可

求出10個(gè)數(shù)的答案之后求個(gè)max即可，因?yàn)榧洗笥谝话?#xff0c;所以正確率為 $1?\frac{1}{2^{10}}$

時(shí)間復(fù)雜度為 $O(n+d^2)$，隨機(jī)的數(shù)再多就會(huì)TLE了

---

code

#include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #define ll long long #define N 1001000 using namespace std; int n,w,now,v[N]; ll x,ans,a[N],s[N]; void get(ll x) {w=0;for(ll i=1;i*i<=x;++i)if(x%i==0){w++,s[w]=i,v[w]=0;if(x/i!=i)w++,s[w]=x/i,v[w]=0;}sort(s+1,s+1+w);return; } ll gcd(ll x,ll y) {return (y?gcd(y,x%y):x); } int main() {scanf("%d",&n);for(int i=1;i<=n;++i)scanf("%lld",&a[i]);for(int g=1;g<=10;++g){x=a[rand()*rand()%n+1];get(x);for(int i=1;i<=n;++i){now=lower_bound(s+1,s+1+w,gcd(x,a[i]))-s;v[now]++;}for(int i=1;i<=w;++i)for(int j=i+1;j<=w;++j)if(s[j]%s[i]==0)v[i]+=v[j];for(int i=w;i>0;--i)if(v[i]>=(n+1)/2)ans=max(ans,s[i]);}printf("%lld",ans);return 0; }

總結(jié)

以上是生活随笔為你收集整理的【随机】Ghd（CF364D）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。