【数学】数数

题目描述

作为一名普及组选手，小 A 喜欢数数。
一天，小 A 学习了排列相关的知识。定义一个长度为 $n$ 的序列 $p_{1,...,n}$ 是一个 $n$ 阶排列，当且仅当 $p_{1,...,n}$ 都是 $[1,n]$ 中的正整数且它们两两不同。
小 A 想数排列。为了让数排列更有趣，小 A 决定加入一个限制：
对于一个 $n$ 阶排列 $p$，小 A 会构造一个长度为 $n$ 的序列 $Q(p)$，其中 $Q(p{)}_{p_i}=i$。小 A 称排列 $p$ 是优秀的，当且仅当 $p$ 的字典序严格小于 $Q(p)$。即存在一个 $i$ 使得 $\forall 1\le j<i,p_j=Q(p{)}_j$ 且 $p_i<Q(p{)}_i$。
现在，小 A 想了一个数 $n$，他希望对于每个 $[1,n]$ 间的 $m$ ，计算好的 $m$ 阶排列数量，在开始数这样的排列数量前，小 A 给了你一个质数 $\mathrm{mod}$，希望你先求出好的 $m$ 阶排列数量对 $\mathrm{mod}$ 取模的结果。
为了避免极其大的输出，设 $v_n$ 表示好的 $n$ 阶排列数量对 $\mathrm{mod}$ 取模的结果，你只需要输出 ${\oplus }_{i=1}^n{v}_i$，即所有 $v_1,\dots ,v_n$ 的异或和。这样小 A 在自己数错的时候就有大约 $1-\frac{1}{\mathrm{mod}}$ 的概率发现自己数错了，并重新数一遍。

$n⩽10^7$。

题解

$Q$ 显然是 $p$ 的逆排列，所以 $p$ 进行两次变换后一定会变回 $p$，如果 $p\ne Q$，则这两个排列一定有一个小于另一个。

所以我们考虑求出 $p=Q$ 的个数，每个长度为 $i-1$ 的排列都可以加一位 $i$ 来构成新的长度为 $i$ 的排列，同时每一个长度为 $i-2$ 的排列都可以理解为两个数形成一个环，对于每种排列会形成 $i-1$ 种。所以：
$$d{p}_n=d{p}_{n-1}+(n-1)\times d{p}_{n-2}$$

于是总答案显然等于：
$$a_i=\frac{1}{2}(i!-d{p}_i)$$

时间复杂度 $\mathcal{O}(n)$。

#include<bits/stdc++.h>
#define int long long
#define mid ((l+r)>>1)
#define fir first
#define sec second
#define lowbit(i) (i&(-i))
using namespace std;
const int N=1e7+5;
inline int read(){
    char op=getchar();
    int w=0,s=1;
    while(op<'0'||op>'9'){
        if(op=='-') s=-1;
        op=getchar();
    }
    while(op>='0'&&op<='9'){
        w=(w<<1)+(w<<3)+op-'0';
        op=getchar();
    }
    return w*s;
}
int mod;
int Mul(int a,int b){return (a*b)%mod;}
int Add(int a,int b){return (a+b)%mod;}
int Dec(int a,int b){return (a-b+mod)%mod;}
int Pow(int a,int k){
    int ans=1;
    while(k){
        if(k&1) ans=Mul(ans,a);
        a=Mul(a,a);
        k>>=1;
    }
    return ans;
}
int inv(int x){return Pow(x,mod-2);}
int a[N],jc[N];
signed main(){
    int n=read(),m=read();
    mod=m,a[0]=a[1]=1;
    const int inv2=inv(2);
    for(register int i=2;i<=n;i++) a[i]=Add(a[i-1],Mul(a[i-2],i-1));
    jc[0]=1;
    for(register int i=1;i<=n;i++) jc[i]=Mul(jc[i-1],i);
    int ans=0;
    for(register int i=1;i<=n;i++) ans^=Mul(Dec(jc[i],a[i]),inv2);
    printf("%lld",ans);
}