E. Modular Stability（2000，数学）

题目传送门

题意： 给你一个n，一个k，在n个数里挑选k个数，使得对于任意非负整数x，对于这k个数的任何排列顺序，然后用x对这个排列一次取模，如果最后取模结果不变，那么称它为稳定的，求稳定数组的个数。

思路： 首先很容易想到，如果数组里存在1，那么其他随便选就行。如果我们有x,y，约定x<y且y%x！=0，那么y%x%y!=y%y%x，所以并不能保证对于所有非负正整数都会成立，即y%x!=0的情况都要排除。所以我们必须保证数组中最小的数是其他所有数的因子，这样对于所有数取模之后，就是一样的。这里要注意在计算组合数的时候先算a[x]-kk+1到a[x]的阶乘再处理逆元，这样总的时间复杂度是nlogn，另外一种计算方法不能保证这个时间复杂度。
这种方法每次都要计算一遍阶乘再处理阶乘逆元，我们可以预处理阶乘和阶乘逆元来节省时间。
与原来代码不同的地方：

void pp()
{
    for(int i=1;i*k<=n;i++)
        a[i]=n/i-1;
    j[0]=1;
    for(int i=1;i<=n;i++)
        j[i]=j[i-1]*i%mod;
    inv[n]=qpow(j[n],mod-2);
    for(int i=n-1;i>=0;i--)
        inv[i]=inv[i+1]*(i+1)%mod;
}
int fun(int x,int kk)
{
    return j[a[x]]*inv[kk]%mod*inv[a[x]-kk]%mod;
}

代码：

#include<bits/stdc++.h>
#define endl '\n'
#define null NULL
#define ls p<<1
#define rs p<<1|1
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define ll long long
#define int long long
#define pii pair<int,int>
#define ull unsigned long long
#define all(x) x.begin(),x.end()
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define ct cerr<<"Time elapsed:"<<1.0*clock()/CLOCKS_PER_SEC<<"s.\n";
char *fs,*ft,buf[1<<20];
#define gc() (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<20,stdin),fs==ft))?0:*fs++;
inline int read(){int x=0,f=1;char ch=gc();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=gc();}
return x*f;}
using namespace std;
const int N=1e6+5;
const int inf=0x7fffffff;
const int mod=998244353;
const double eps=1e-6;
const double PI=acos(-1);
int a[N],n,k;
void pp()
{
    for(int i=1;i*k<=n;i++)//小于n的i的倍数的个数
        a[i]=n/i-1;
}
ll qpow(ll a,ll b)
{
    ll res=1;
    while(b)
    {
        if(b&1)
            res=res*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return res;
}
int fun(int x,int kk)
{
    int res=1;
    for(int i=a[x]-kk+1;i<=a[x];i++)
        res=res*i%mod;
    for(int i=1;i<=kk;i++)
        res=res*qpow(i,mod-2)%mod;//求组合数的时候处理逆元
    return res;
}
signed main()
{
    cin>>n>>k;
    if(n<k)
        cout<<0<<endl;
    else
    {
        pp();int res=0;
        for(int i=1;i*k<=n;i++)
            res=(res+fun(i,k-1))%mod;
        cout<<res<<endl;
    }
}