CF1359E Modular Stability 组合

分析一波发现如果一个序列合法,那么必须满足所有数都是最小值的倍数.  

加上组合数就行了. 

code: 

#include <bits/stdc++.h>  
#define ll long long 
#define mod 998244353 
#define N 500006   
#define setIO(s) freopen(s".in","r",stdin)   
using namespace std;     
int fac[N],inv[N];    
int qpow(int x,int y) 
{
    int tmp=1; 
    for(;y;y>>=1,x=(ll)x*x%mod)  if(y&1) tmp=(ll)tmp*x%mod;  
    return tmp; 
}
int INV(int x) { return qpow(x,mod-2); }   
void init() 
{
    fac[0]=inv[0]=1;  
    for(int i=1;i<N;++i) fac[i]=(ll)fac[i-1]*i%mod,inv[i]=INV(fac[i]);  
}
int C(int x,int y) { return (x<y||x<0||y<0)?0:(ll)fac[x]*inv[y]%mod*inv[x-y]%mod; }   
int main() 
{ 
    // setIO("input"); 
    int n,k;  
    init();   
    scanf("%d%d",&n,&k);   
    int ans=0; 
    for(int i=1;i<=n;++i) 
    { 
        (ans+=(ll)C(n/i-1,k-1)%mod)%=mod; 
    }        
    printf("%d\n",ans); 
    return 0; 
}

  

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值