数学专题训练3 数论2

zhezhidashi

于 2023-07-16 17:57:44 发布

阅读量156

点赞数

分类专栏： ACM题目整理文章标签：算法

本文链接：https://blog.csdn.net/qq_45812711/article/details/131753063

版权

ACM题目整理专栏收录该内容

76 篇文章 6 订阅

订阅专栏

1. J-Product of GCDs_2021牛客暑期多校训练营2 (nowcoder.com)

给定一个集合，求其所有大小为 $k$ 的子集的 $g c d$ 之积. 子集的 $g c d$ 定义为子集所有元素的最大公因数.

思路很简单，每个质因数都分开考虑. 因为 $p$ 不一定是质数所以要求欧拉函数降幂。由于 $n * k$ 很小所以可以 $O (nk)$ 递推法求组合数.

这个题卡线性，有一个做法是保存至少含有 $p^c$ 的数字有多少个，记为 $c n t (p, c)$ ，然后这样的方案数是 $C_{cnt(p, c)}^{k}$ ，每个质数把幂次从大到小枚举，容斥一下。这样质数的个数是 $\frac{n}{\log n}$ ，质数的幂是 $\log n$ . 因此枚举一遍是 $O (n)$

这里放一个超时的板子，希望能够改一改卡过去.

#include<bits/stdc++.h>
using namespace std;
const int N = 80010;
typedef long long ll;
int prime[N], cnt, st[N];
ll C[40][40010], f[N];
int n, k, a[N];
ll phi, ans, P;

typedef pair<int, int> PII;

#define x first
#define y second

vector<PII> divs[N];

ll mod_pow(ll x, ll n, ll mod)
{
    ll res = 1;
    while(n)
    {
        if(n & 1) res = (__int128)res * x % mod;
        x = (__int128)x * x % mod;
        n >>= 1;
    }
    return res;
}

void sieve(int n)
{
    st[1] = 1;
    for(int i = 2; i <= n; i++)
    {
        if(!st[i]) st[i] = i, prime[cnt++] = i;
        for(int j = 0; prime[j] <= n / i; j++)
        {
            st[i * prime[j]] = prime[j];
            if(i % prime[j] == 0) break;
        }
    }
    for(int i = 1; i <= n; i++)
    {
        int m = i;
        while(m > 1)
        {
            ll cnt = 0, p = st[m];
            while(m % p == 0)
            {
                m /= p;
                cnt++;
            }
            divs[i].push_back({p, cnt});
        }
    }
}

ll get_phi(ll n)
{
    ll res = n;
    for(ll i = 2; i <= n / i; i++)
    {
        if(n % i == 0)
        {
            res = res / i * (i - 1);
            while(n % i == 0) n /= i;
        }
    }
    if(n > 1) res = res / n * (n - 1);
    return res;
}

int CNT[N][20];
void solve(ll p)
{
    CNT[p][0] = n;
    for(int i = 1; i < 20; i++) CNT[p][0] -= CNT[p][i];
    if(CNT[p][0] == n) return;

    int remain = n;
    for(int i = 0; i < 20; i++)
    {
        if(remain < k) break;
        ans = (__int128)ans * mod_pow(p, (__int128)i * (f[remain - 1] -
                                        ((remain - CNT[p][i] - 1) < 0 ? 0 : (f[remain - CNT[p][i] - 1])) + phi) % phi, P) % P;
        remain -= CNT[p][i];
    }
}

int main()
{
    sieve(N - 1);
    int T;
    scanf("%d", &T);
    while(T--)
    {
        scanf("%d%d%lld", &n, &k, &P);

        ans = 1;
        fill(f, f + n + 1, 0);
        for(int i = 0; i < cnt; i++) fill(CNT[prime[i]], CNT[prime[i]] + 20, 0);
        
        for(int i = 1; i <= n; i++) scanf("%d", &a[i]);

        phi = get_phi(P);

        for(int i = 0; i < n; i++)
        {
            int sz = min(i, k);
            for(int j = 0; j <= sz; j++)
            {
                if(j == 0) C[j][i] = 1;
                else C[j][i] = (C[j][i - 1] + C[j - 1][i - 1]) % phi;
            }
        }

        f[k - 1] = 1;
        for(int i = k; i <= n; i++) f[i] = (f[i - 1] + C[k - 1][i]) % phi;

        for(int i = 1; i <= n; i++)
        {
            for(auto p : divs[a[i]]) CNT[p.x][p.y]++;
        }

        for(int i = 0; i < cnt; i++)
        {
            solve(prime[i]);
        }

        printf("%lld\n", ans);
    }
    return 0;
}

2. H-Convolution_2021牛客暑期多校训练营4 (nowcoder.com)

比较难

3. G-Greater Integer, Better LCM_2021牛客暑期多校训练营5 (nowcoder.com)

比较难

4. Pass! - HDU 6956 - Virtual Judge (vjudge.net)

BSGS+卡常

5. I love permutation - HDU 6970 - Virtual Judge (vjudge.net)

6. Problem - 1548B - Codeforces

在一个序列上找到最长的连续子区间，使得存在 $\ge 2$ ，子区间内所有数字模 $m$ 同余.

如果每个数字模 $m$ 的值都一样，设模数为 $d$ ，那么说明 $a_i - d$ 模 $m$ 都为 $0$ ，说明 $\gcd(a-d, b-d) =\gcd(a-d, |a - b|) > 1$ . 就是在差分数组上找到最长的一段子区间他们的最大公约数不为 $1$ . 然后答案就是这个区间的长度加1. 考虑到最大公约数的单调性用双指针即可.