CF1359E Modular Stability

Jennie

这个东西需要证明,不过可以考虑先拿出两个来,$a_1,a_2\quad a_1<a_2$

这有什么用处?我们要找到对于任意的x $x\%a_1=x\%a_2\%a_1$

为什么呢,考虑一下他们的大小关系,然后在%意义下,设

$k*a_2+y\equiv\quad y$

而k是任意的,所以式子要恒成立,前面的那一部分%意义下为0就要恒成立

所以$a_2$一定是$a_1$的倍数。推广一下,整个序列都是$a_1$的倍数

所以说就是个枚举$a_1$并且求组合数的问题

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define int long long
using namespace std;
template<class T>
void read(T &now){
    now=0;
    char c=getchar();
    int f=1;
    while((!isdigit(c))){
        if(c=='-') f=-1;
    //  cout<<isdigit(c)<<" "<<c<<" ";
        c=getchar();
    }
    while(isdigit(c)){
        now=(now<<1)+(now<<3)+c-'0';
        c=getchar();
    }
    now*=f;
}
int n,k;
int fac[500005];
int inv[500005];
int mod=998244353;
void cal(){
    fac[0]=1;
    inv[0]=inv[1]=1;
    for(int i=1;i<=n;++i){
        fac[i]=(fac[i-1]*i)%mod;
    }
    for(int i=2;i<=n;++i){
        inv[i]=(mod-mod/i)*inv[mod%i]%mod;
    }
    for(int i=2;i<=n;++i){
        inv[i]=(inv[i-1]*inv[i])%mod;
    }
}
int ans;
signed main(){
    read(n);read(k);
    cal();
    for(int i=1;i<=n;++i){
        int tot=n/i;
        if(tot<k) break;
        if(k==1){
            ans+=1;
            ans%=mod;
        }else{
            tot--;
            ans+=(fac[tot]*inv[k-1]%mod*inv[tot-k+1]%mod)%mod;
            ans%=mod;
        }
    }
    cout<<ans;
    return 0;
}
暂无评论

发送评论 编辑评论


				
|´・ω・)ノ
ヾ(≧∇≦*)ゝ
(☆ω☆)
(╯‵□′)╯︵┴─┴
 ̄﹃ ̄
(/ω\)
∠( ᐛ 」∠)_
(๑•̀ㅁ•́ฅ)
→_→
୧(๑•̀⌄•́๑)૭
٩(ˊᗜˋ*)و
(ノ°ο°)ノ
(´இ皿இ`)
⌇●﹏●⌇
(ฅ´ω`ฅ)
(╯°A°)╯︵○○○
φ( ̄∇ ̄o)
ヾ(´・ ・`。)ノ"
( ง ᵒ̌皿ᵒ̌)ง⁼³₌₃
(ó﹏ò。)
Σ(っ °Д °;)っ
( ,,´・ω・)ノ"(´っω・`。)
╮(╯▽╰)╭
o(*////▽////*)q
>﹏<
( ๑´•ω•) "(ㆆᴗㆆ)
😂
😀
😅
😊
🙂
🙃
😌
😍
😘
😜
😝
😏
😒
🙄
😳
😡
😔
😫
😱
😭
💩
👻
🙌
🖕
👍
👫
👬
👭
🌚
🌝
🙈
💊
😶
🙏
🍦
🍉
😣
Source: github.com/k4yt3x/flowerhd
颜文字
Emoji
小恐龙
花!
上一篇
下一篇