ural 1682

最新推荐文章于 2022-02-09 12:01:41 发布

xcszbdnl

最新推荐文章于 2022-02-09 12:01:41 发布

阅读量611

点赞数

CC 4.0 BY-SA版权

分类专栏： ACM_并查集

本文链接：https://blog.youkuaiyun.com/xcszbdnl/article/details/38520557

ACM_并查集专栏收录该内容

3 篇文章

订阅专栏

本文介绍了一种基于模k图构建最小环的方法，通过枚举数字并利用连边条件a+b^2==0modk或b+a^2==0modk来形成环状结构，实现了寻找特定条件下最小环的算法。

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

从1写到无穷，每写1个数可以和前面的数连边，连边有两个条件a + b ^ 2 == 0 mod k || b + a ^ 2 == 0 mod k

首先写下的数是b，满足条件a + b ^ 2 的话，那么可以和k - b ^ 2, 2 * k - b ^ 2一直连下去

如果满足条件是b + a ^ 2，那么要找到符合条件的a，对每插入的一个数把它的平方对k取模保留下来，那么就可以迅速查找到哪些数满足b + a ^ 2 == 0 mod k

因为1, k - 1, 2k - 1一定是一个环，所以数只要枚举到2k - 1就可以了。

#include <cstdio>
#include <cstring>
#include <vector>

using namespace std;
const int maxn = 200007;

int fa[maxn], pre[maxn];
vector<int> mod[maxn];

int findFa(int x) {
    if (x == fa[x]) {
        return x;
    }
    return fa[x] = findFa(fa[x]);
}
int main() {
//    freopen("in.txt", "r", stdin);
    int k;
    scanf("%d", &k);
    if (k == 1) {
        printf("3\n");
        return 0;
    }
    else if (k == 2) {
        printf("5\n");
        return 0;
    }
    for (int i = 0; i < 2 * k; i++) {
        fa[i] = i;
        mod[i].clear();
    }
    int ans = 0;
    for (int i = 1; i < 2 * k; i++) {
        long long temp = (long long)i *(long long)i;
        temp %= k;
        int cnt_mod = (k - temp);
        int flag = 0;
        while (cnt_mod < i) {
            int x = findFa(i);
            int y = findFa(cnt_mod);
            if (x != y) {
                fa[x] = y;
            }
            else {
                flag = 1;
                break;
            }
            pre[cnt_mod] = i;
            cnt_mod += k;
        }
        cnt_mod = (((k - i) % k) + k) % k;
        for (int j = 0; j < mod[cnt_mod].size(); j++) {
            int s = mod[cnt_mod][j];
            if (pre[s] == i) {
                continue;
            }
            int x = findFa(i);
            int y = findFa(s);
            if (x != y) {
                fa[x] = y;
            }
            else {
                flag = 1;
                break;
            }
        }
        if (flag) {
            ans = i;
            break;
        }
        temp = (long long)i * (long long)i;
        temp %= k;
        cnt_mod = temp;
        mod[cnt_mod].push_back(i);
    }
    printf("%d\n", ans);
    return 0;
}