2019 icpc ecfinal H.King 思维模拟随机数

binggui2

于 2022-01-09 16:05:25 发布

阅读量536

点赞数

CC 4.0 BY-SA版权

分类专栏：思维文章标签：算法 c语言

本文链接：https://blog.youkuaiyun.com/binggui2/article/details/122394832

思维专栏收录该内容

4 篇文章

订阅专栏

本文介绍了一道思维题的编程解决方案，通过使用C++实现，并针对随机数生成进行了优化以提高程序的正确率。该方案利用了快速幂运算来解决特定问题，并通过多次尝试寻找最优解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

一个思维题，切入点是n/2

写代码后交疯狂wa18，然后发现rand函数的上限是32767，太小了，换了一个随机数就过了。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int N = 200020;

int n, mod;

std::mt19937 rnd(time(0));
int a[N];

ll qmi(ll a, int k) {
    ll res = 1;
    while(k) {
        if(k & 1) res = (res * a) % mod;
        a = (a * a) % mod;
        k >>= 1;
    }
    return res;
}

int main() {
    int t;
    scanf("%d", &t);
    while(t --) {
        scanf("%d%d", &n,&mod);
        for(int i = 1;i <= n;i ++) scanf("%d", &a[i]);

        int ans = 0;
        for(int Case = 1;Case <= 500;Case ++) {
            ll x = rnd() % n + 1;

            for(ll l = 1;l <=2 && l + x <= n;l ++) {
                int num = 2;
                ll q = (a[l+x] * qmi(a[x], mod-2)) % mod;
                ll last = (a[l+x] * q) % mod;
                for(int i = l + x + 1;i <= n;i ++) {
                    if(a[i] == last) {
                        num ++;
                        last = (last * q) % mod;
                    }
                }

                q = qmi(q, mod-2);
                last = (a[x] * q) % mod;

                for(int i = x-1;i >= 1;i --) {
                    if(a[i] == last) {
                        num ++;
                        last = (last * q) % mod;
                    }
                }
                ans = max(ans, num);
            }
        }
        if(ans >= (n+1)/2) printf("%d\n", ans);
        else printf("%d\n", -1);
    }
    return 0;
}