2022 icpc杭州站 A. Modulo Ruins the Legend

最新推荐文章于 2025-01-08 11:30:25 发布

yi_heng123

最新推荐文章于 2025-01-08 11:30:25 发布

阅读量1.5k

点赞数 6

分类专栏：数论文章标签：算法 c++

本文链接：https://blog.csdn.net/yi_heng123/article/details/128229153

版权

数论专栏收录该内容

1 篇文章

订阅专栏

2022 icpc杭州站 A. Modulo Ruins the Legend

Problem - A - Codeforces

在这里插入图片描述

思路

题目可以转换为求 $\cdot n + d_t \cdot \frac{n(n + 1)}{2} + sum) \% m （1）$ 的最小值

$gcd(n,\frac{n(n + 1)}{2})$

于是 $（ 1 ）$ 转换为
$\begin{aligned}(k \cdot d + sum) \% m &= (k\cdot d \%m + sum \%m) \%m \\&=(k\cdot d + t\cdot m + sum \%m) \%m \\&g = gcd(d, m), k\cdot d + t\cdot m = z\cdot g, \ sum_2 = sum \% m\\&=(z \cdot g + sum_2) \% m\end{aligned}$
因为 $\le m, sum_2 < m$ , 所以 $\cdot g \ge m - sum_2$ 时候能得到最小值， $\lceil \frac{m - sum_2}{g} \rceil$

证明：

$z\cdot g + sum_2 < m$ 时， $\cdot g + sum_2) \% m$ 中最小都会有一个 $g$

$z\cdot g + sum_2 > m$ 时， $\cdot g + sum_2) \% m$ 最小值能够小于 $g$

求出 z 之后，再一步一步的代回去，用 exgcd 求出 $s, d_t$ 即可

最小值就是 $z\cdot g + sum_2 - m$

#include <bits/stdc++.h>
#define endl "\n"
#define INF 0x3f3f3f3f
#define LLINF 0x3f3f3f3f3f3f3f3f
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0)
using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;

ll exgcd(ll a, ll b, ll& x, ll& y) {
    if (!b) {
        x = 1, y = 0;
        return a;
    }
    ll d = exgcd(b, a % b, x, y);
    ll tx = x;
    x = y, y = tx - y * (a / b);
    return d;
}

signed main() {
    IOS;
    ll n, mod; cin >> n >> mod;
    ll sum = 0;
    for (int i = 1; i <= n; i++) {
        ll t; cin >> t;
        sum += t;
    }

    ll a = n, b = n * (n + 1) / 2;
    ll s, dt;
    ll d = exgcd(a, b, s, dt);
    sum %= mod;

    ll k, t;
    ll g = exgcd(d, mod, k, t);

    ll z = (mod - sum + g - 1) / g;
    (k *= z) %= mod;
    
    s = ((s % mod * k) % mod + mod) % mod, dt = ((dt % mod * k) % mod + mod) % mod;

    cout << (z * g + sum - mod) << endl;
    cout << s << " " << dt << endl;
}