Description
给你一个序列s,你把这个序列的所有不同排列按字典序排列后,求s的排名mod m
题解:
完全不会啊……代码都几乎是照着栋老师打的……这题思维难度不算太大(但我依然没有想出来),就是一位一位固定,算每一位上比当前排列小的有多少个,然后就可以推出式子(具体的式子可以看别人的blog),然后用中国剩余定理合并就行了,但这个东西我只在学扩展lucas的时候用过,已经忘了,于是又粗略地学了一次,要注意当要乘或除的数与当前的模数有公因子的话,就一直除,直到没有公因子,统计现在还需要乘或除多少个
p
p
<script type="math/tex" id="MathJax-Element-17">p</script>,然后再乘或除。
代码:
#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
const int Maxn=300010;
const int inf=2147483647;
LL read()
{
LL x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
return x*f;
}
struct lsh{int v,id;}A[Maxn];
bool cmp(lsh a,lsh b){return a.v<b.v;}
int n,a[Maxn];LL mod;
int s[Maxn];
void add(int x,int t){for(;x<=n;x+=(x&-x))s[x]=(s[x]+t)%mod;}
int query(int x){int re=0;for(;x;x-=(x&-x))re=(re+s[x])%mod;return re;}
LL exgcd(LL a,LL b,LL &x,LL &y)
{
if(!a){x=0,y=1;return b;}
LL tx,ty,d=exgcd(b%a,a,tx,ty);
x=ty-(b/a)*tx;
y=tx;
return d;
}
LL Pow(LL x,int y,LL MOD)
{
if(!y)return 1%MOD;
if(y==1)return x%MOD;
LL t=Pow(x,y>>1,MOD),re=t*t%MOD;
if(y&1)re=re*x%MOD;return re;
}
int p1[Maxn],p2[Maxn],lp=0;
int p,pm,inv[Maxn];
void pre()
{
memset(inv,0,sizeof(inv));
inv[1]=1;int phi=pm/p*(p-1);
for(int i=2;i<=n;i++)if(i%p)inv[i]=Pow(i,phi-1,pm);
}
int tot[Maxn],mx;
LL get()
{
LL v=1,re=0;int cnt=0;pre();
for(int i=1;i<n;i++)
{
tot[a[i]]++;
int t=i;
while(t%p==0)t/=p,cnt++;
v=v*(LL)(t)%pm;
}tot[a[n]]++;
for(int i=1;i<=mx;i++)
{
add(i,tot[i]);
for(int j=2;j<=tot[i];j++)
{
int t=j;
while(t%p==0)t/=p,cnt--;
v=v*(LL)(inv[t])%pm;
}
}
for(int i=1;i<=n;i++)
{
LL t1=query(a[i]-1);
LL t2;if(cnt<0)t2=Pow(p,-cnt,inf);else t2=Pow(p,cnt,pm);
LL t3;if(cnt<0)t3=t1/t2*v%pm;else t3=t1*t2%pm*v%pm;
re=(re+t3)%pm;
if(i!=n)
{
int t;
t=n-i;while(t%p==0)t/=p,cnt--;
v=v*inv[t]%pm;
t=tot[a[i]];while(t%p==0)t/=p,cnt++;
v=v*t%pm;
}
tot[a[i]]--;add(a[i],-1);
}
return re;
}
LL Inv(LL a,LL mod)
{
LL x,y;exgcd(a,mod,x,y);
return (x%mod+mod)%mod;
}
LL solve()
{
int t=mod;LL re=0;
for(int i=2;i*i<=mod;i++)
if(t%i==0)
{
p1[++lp]=i;p2[lp]=0;
while(t%i==0)t/=i,p2[lp]++;
}
if(t!=1)p1[++lp]=t,p2[lp]=1;
for(int i=1;i<=lp;i++)
{
p=p1[i];pm=1;
for(int j=1;j<=p2[i];j++)pm*=p;
LL c=get();
re=(re+Inv(mod/pm,pm)*c%mod*(LL)(mod/pm)%mod)%mod;
}return (re+1)%mod;
}
int main()
{
n=read(),mod=read();
for(int i=1;i<=n;i++)A[i].v=read(),A[i].id=i;
sort(A+1,A+1+n,cmp);
a[A[1].id]=1;int cnt=1;
for(int i=2;i<=n;i++)a[A[i].id]=(cnt+=((A[i].v==A[i-1].v)?0:1));
mx=cnt;
printf("%lld",solve());
}