2021“MINIEYE杯”中国大学生算法设计超级联赛（9）——H - Integers Have Friends 2.0

H-w-H

于 2021-08-31 17:33:57 发布

阅读量140

点赞数 1

分类专栏：题解思维文章标签：算法数据结构

本文链接：https://blog.csdn.net/weixin_45363113/article/details/120022933

版权

题解同时被 2 个专栏收录

82 篇文章 0 订阅

订阅专栏

思维

7 篇文章 0 订阅

订阅专栏

该博客探讨了一道编程题目，涉及数组处理和数论概念。通过证明gcd(∣x−y∣,∣z−y∣)>1在模同余条件下的重要性，提出了一个算法策略，即随机选取数组中的两个数，计算它们差值的质因子，并寻找模数相同的元素个数，从而找到最长的子序列。最后，通过不断优化这个过程来求解问题的最大值。

摘要由CSDN通过智能技术生成

H-Integers_Have_Friends2.0
- 题意：
- 思路：

H-Integers_Have_Friends2.0

题意：

给出一个数组 $a$ ，找到一个最长的子序列 $b$ ，存在 $m (m > 1)$ 使 $b_i~mod~m=b_j~mod~m$ 。

思路：

我们先了解 $x m o d m = y m o d m = z m o d m$ ，则 $g c d (∣ x - y ∣, ∣ z - y ∣) > 1$

证明：

假设 $g c d (∣ x - y ∣, ∣ z - y ∣) = 1$ 存在 $m (m > 1)$ 使得 $x m o d m = y m o d m = z m o d m$ 。

则： $x=k_1m+u,y=k_2m+u,z=k_3m+u$ ，所以 $x-y|=|(k_1-k2)m|，|z-y|=|(k_3-k_2)m|$ ，因为 $m > 1$ ，所以 $g c d (∣ x - y ∣, ∣ z - y ∣) > = m > 1$

所以假设不成立，得证 $x m o d m = y m o d m = z m o d m$ ，则 $g c d (∣ x - y ∣, ∣ z - y ∣) > 1$ 。

现在我们从a数组中随机出两个数字，它们两个同时在最长子序列的概率很大，所以我们就随机选出两个数字，验证它们差值的全部质因子，返回在a数组中模数相同的个数。再取一个最大值。

#include <bits/stdc++.h>

using namespace std;
const int N = 2e6+10;
const int inf = 1e9+10;
typedef long long ll;
// #define int long long
const int K=30;
mt19937_64 gen(time(0));
ll a[N];
int p[N], tot;
bitset<N> st;
void init() {
    for(int i=2; i<N; i++) {
        if(!st[i]) p[tot++] = i;
        for(int j=0; j<tot&&1ll*i*p[j]<N; j++) {
            st[i*p[j]] = 1;
            if(i % p[j] == 0) break;
        }
    }
}
int n;
int C(ll pp, ll y) {
    int cnt = 0;
    for(int i=1; i<=n; i++) {
        if(a[i] % pp == y) cnt++; 
    }
    return cnt;
}
int main() {
#ifndef ONLINE_JUDGE
    freopen("in.txt", "r", stdin);
    freopen("out.txt", "w", stdout);
#endif
    init();
    int t;
    scanf("%d", &t);
    while(t--) { 
        int k = 30, ans = 1;
        
        scanf("%d", &n);
        for(int i=1; i<=n; i++) scanf("%lld", &a[i]);
        while(k--) {
            int x, y;
            while(1) {
                x = gen() % n + 1;
                y = gen() % n + 1;
                if(x != y) break;
            }
            ll tmp = abs(a[x] - a[y]);
            for(int i=0; 1ll*p[i]*p[i]<=tmp; i++) {
                if(tmp % p[i] == 0) {
                    ans = max(ans, C(p[i], a[x]%p[i]));
                    while(tmp % p[i] == 0) tmp /= p[i];
                }
            }
            if(tmp > 1) ans = max(ans, C(tmp, a[x]%tmp));
        }
        printf("%d\n", ans);
    }
}