codeforces E. Modular Stability

题目

题意：

给你$n,k$表示$1\le a_1,a_2,.....a_k\le n$，然后是否存在这样的一个$a$序列满足$$(((x\%a_1)\%a_2)\dots \%a_{k-1})\%a_k=(((x\%a_{p1})\%a_{p2})\dots \%a_{pk-1})\%a_p$$
$x$是一个正整数。

思路：

要满足上述公式的话，那么最后的值一定要为$0$，如果不为$0$的话，那么$x$肯定会出现使这个公式不成立的，比如$2,3$和$3,2$，那么当$x=3$的时候，就不同了，所以只有整除后才会无论怎么取模都为$0$，都符合公式，那么最后能够取模为$0$的话，那么就是倍数的关系，比如说$3,6,9$最后无论怎么取都相当于$\%3$,所以最后能取得就是倍数得关系，我们可以遍历一下$1$~$n$，那么此时对于每一个$i$总共能够取$n/i$个数字，然后这其中要有一个最小$i$这样无论怎么取都是$i$得倍数，所以$an{s}_i=C_{(n/i-1)}^{k-1}$，最后的公式就是${\sum }_{i=1}^{i\le n}{C}_{(n/i-1)}^{k-1}$。

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <vector>
#include <string>
#include <cmath>
#include <set>
#include <map>
#include <deque>
#include <stack>
#include <cctype>
using namespace std;
typedef long long ll;
typedef vector<int> veci;
typedef vector<ll> vecl;
typedef pair<int, int> pii;
template <class T>
inline void read(T &ret) {
    char c;
    int sgn;
    if (c = getchar(), c == EOF) return ;
    while (c != '-' && (c < '0' || c > '9')) c = getchar();
    sgn = (c == '-') ? -1:1;
    ret = (c == '-') ? 0:(c - '0');
    while (c = getchar(), c >= '0' && c <= '9') ret = ret * 10 + (c - '0');
    ret *= sgn;
    return ;
}
int xxxxxxxxx = 1;
inline void outi(int x) {if (x > 9) outi(x / 10);putchar(x % 10 + '0');}
inline void outl(ll x) {if (x > 9) outl(x / 10);putchar(x % 10 + '0');}
inline void debug(ll x) {cout << xxxxxxxxx++ << " " << x << endl;}
inline void debugs(string s) {cout << s << endl;}
const int maxn = 5e5 + 10;
const int mod = 998244353;
ll fac[maxn] = {0};
ll quickpow(ll a, ll b) {
    ll ans = 1, base = a;
    while (b) {
        if (b & 1) ans = ans * base % mod;
        base = base * base % mod;
        b >>= 1;
    }
    return ans;
}
ll C(ll n, ll m) {
    if (n < m) return 0;
    return fac[n] * quickpow(fac[m] * fac[n - m] % mod, mod - 2) % mod;
}
int main() {
    fac[0] = fac[1] = 1;
    for (int i = 1; i <= maxn - 1; i++) fac[i] = fac[i - 1] * i % mod;
    ll n, k, sum = 0;
    read(n), read(k);
    for (int i = 1; i <= n; i++) {
        if (n / i - 1 >= k - 1) sum = (sum + C((ll)(n / i - 1), k - 1)) % mod;
        else break;
    }
    outl(sum);
    return 0;
}