标签：int res ll x% Stability Modular inv mod

题意：

求出满足要求的数组个数：

\(1\leq a[i] \leq n，1\leq a[1]<a[2]<...<a[k]\leq n\)，元素个数为：\(k\)

分析：

问题取决于最小的元素 \(a[1]\) 。
设 \(x=n*a[1]+m\)，则\(x%a[1]=m\)，即 \(x%a[1]%a[2]%...%a[n]=m\)。
要使的 \(x%a[i]=m\)，即 \(x=a[i]*n'+m\)，所以 \(a[i]\) 是 \(a[1]\) 的倍数。
所以，可以枚举 \(a[1]\) 的取值，剩下的 \(k-1\) 个数在其倍数中选择，最终的结果为：\(\sum_{i=1}^{n}{C(\frac{n}{i}-1,k-1)}\)。
需要预处理出阶乘逆元：\(inv[i]=inv[i+1]*(i+1)\)

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=998244353;
const int N=5e5+5;
ll inv[N],fac[N];
ll power(ll a,ll b)
{
    ll res=1;
    while(b)
    {
        if(b&1) res=res*a%mod;
        b>>=1;
        a=a*a%mod;
    }
    return res%mod;
}
void init()
{
    fac[0]=1;
    for(int i=1;i<=5e5;i++)
        fac[i]=(fac[i-1]*i)%mod;
    int t=5e5;
    inv[t]=power(fac[t],1LL*(mod-2));
    for(int i=t-1;i>=0;i--)//注意要到0,会用到inv[0]
        inv[i]=inv[i+1]*(i+1)%mod;
}
ll C(int x,int y)
{
    if(x<y) return 0;
    ll res=1;
    y=min(x-y,y);
    for(int i=x;i>=x-y+1;i--)
        res=res*i%mod;
    res=res*inv[y]%mod;
    return res;
}
int main()
{
    int n,k;
    scanf("%d%d",&n,&k);
    ll ans=0;
    init();
    for(int i=1;i<=n;i++)
        ans=(ans+C(n/i-1,k-1))%mod;
    printf("%lld\n",ans);
    return 0;
}

标签：int,res,ll,x%,Stability,Modular,inv,mod
来源： https://www.cnblogs.com/1024-xzx/p/12994919.html

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