ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

[Luogu4916]魔力环[Burnside引理、组合计数、容斥]

2019-03-31 16:52:08  阅读:382  来源: 互联网

标签:Luogu4916 phi return int LL 容斥 long Burnside mod


题意

题目链接

分析

orz yyb

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
#define go(u) for(int i = head[u], v = e[i].to; i; i=e[i].lst, v=e[i].to)
#define rep(i, a, b) for(int i = a; i <= b; ++i)
#define pb push_back
#define re(x) memset(x, 0, sizeof x)
inline int gi() {
    int x = 0,f = 1;
    char ch = getchar();
    while(!isdigit(ch)) { if(ch == '-') f = -1; ch = getchar();}
    while(isdigit(ch)) { x = (x << 3) + (x << 1) + ch - 48; ch = getchar();}
    return x * f;
}
template <typename T> inline bool Max(T &a, T b){return a < b ? a = b, 1 : 0;}
template <typename T> inline bool Min(T &a, T b){return a > b ? a = b, 1 : 0;}
const int N = 2e5 + 7, mod = 998244353;
int n, m, k, pc;
LL fac[N], inv[N], invfac[N];
int phi[N], low[N], pri[N];
void sieve(int n) {
    phi[1] = 1;int to;
    for(int i = 2; i <= n; ++i) {
        if(!low[i]) low[i] = i, phi[i] = i - 1, pri[++pc] = i;
        for(int j = 1; j <= pc && (to = i * pri[j]) <= n; ++j) {
            low[to] = pri[j];
            if(i % pri[j] == 0) {
                phi[to] = phi[i] * pri[j];
                break;
            }
            phi[to] = phi[i] * (pri[j] - 1);
        }
    }
}
LL C(int n, int m) {
    return fac[n] * invfac[m] % mod * invfac[n - m] % mod;
}
void add(LL &a, LL b) {
    a += b;
    if(a >= mod) a -= mod;
}
LL put(int a, int b) {
    if(!b) return 0;
    LL res = 0;
    for(int i = 0; i <= a / (k + 1); ++i)
        add(res, (i & 1 ? mod - 1 : 1) * C(b, i) % mod * C(a - i * (k + 1) + b - 1, b - 1) % mod);
    return res;
}
LL solve(int d) {
    int c = m / (n / d); LL res = 0;
    if(c <= k) return C(d, c);
    for(int i = 0; i <= k; ++i)
        add(res, (i + 1) * put(c - i, d - c - 1) % mod);
    return res;
}
int main() {
    n = gi(), m = gi(), k = gi();
    if(n == m)
        return puts(k < m ? "0" : "1"), 0;
    sieve(n);
    inv[1] = fac[0] = invfac[0] = 1;
    rep(i, 1, 2 * n) {
        if(i ^ 1) inv[i] = (mod - mod / i) * inv[mod % i] % mod;
        fac[i] = fac[i - 1] * i % mod;
        invfac[i] = invfac[i - 1] * inv[i] % mod;
    }
    LL ans = 0;
    for(int i = 1; i <= n; ++i) if(n % i == 0 && m % (n / i) == 0) 
        add(ans, solve(i) * phi[n / i] % mod);
    printf("%lld\n", ans * inv[n] % mod);
    return 0;
}

标签:Luogu4916,phi,return,int,LL,容斥,long,Burnside,mod
来源: https://www.cnblogs.com/yqgAKIOI/p/10631918.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有