【Codeforces 364D】Ghd

最新推荐文章于 2020-12-05 21:22:14 发布

_tan45

最新推荐文章于 2020-12-05 21:22:14 发布

阅读量199

点赞数

分类专栏： Codeforces 文章标签：随机化 Codeforces

本文链接：https://blog.csdn.net/C20181503csy/article/details/102772302

版权

Codeforces 专栏收录该内容

3 篇文章 0 订阅

订阅专栏

给你 $n$ 个数，求一个尽量大的数，使得数列中有超过 $\frac{n}{2}$ 个数能被该数整除。
数 $a_i \leq 10^{12}$

发现每一个数出现在能被整除的集合中的概率均 $≥\frac{1}{2}$ ，于是考虑随机。每一次猜测一个数在这个集合中，然后计算与其它数的gcd。最后查找一个最大的出现次数超过 $\frac{n}{2}$ 的gcd即可。（若 $r e s [i] ∣ r e s [j]$ 则 $s u m [i] + = s u m [j]$ ）
我们分析一下这样做的正确率。每一次选对的概率 $≥\frac{1}{2}$ ，所以说一次选错的概率 $\leq \frac{1}{2}$ ，即随机 $x$ 次的成功率为 $\Large 1-\frac{1}{2^x}$ 。当 $x$ 取约 $14$ 时，其成功率可以接受。

        #include<bits/stdc++.h>
        #define reg register
        #define ll long long
        using namespace std;
        const int mn = 1000005;
        ll a[mn], res[mn];
        int sum[mn];
        inline ll gcd(ll a, ll b) {return b ? gcd(b, a % b) : a;}
        int main()
        {
            srand(time(NULL));
            int n;
            scanf("%d", &n);
            for(int i = 1; i <= n; i++)
                scanf("%I64d", &a[i]);
            ll ans = 1;
            for(int T = 0; T < 13; T++)
            {
                int pos = ((rand() << 15) + rand()) % n + 1;
                if(a[pos] <= ans) continue;
                ll maxs = 1;
                for(reg int i = 1; i <= n; i++)
                    res[i] = gcd(a[i], a[pos]), maxs = max(res[i], maxs);
                if(maxs <= ans) continue;
                sort(res + 1, res + 1 + n); int cnt = 0;
                for(reg int i = 1; i <= n; i++)
                {
                    if(res[i] != res[i - 1])
                        res[++cnt] = res[i], sum[cnt] = 0;
                    ++sum[cnt];
                }
                for(reg int i = cnt; i; i--)
                {
                    if(ans >= res[i]) break;
                    int tmp = 0;
                    for(int j = i; j <= cnt; j++)
                        if(res[j] % res[i] == 0) tmp += sum[j];
                    if((tmp << 1) >= n) {ans = res[i]; break;}
                }
            }
            printf("%I64d\n", ans);
        }