ural 1682

从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;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值