康托展开
考虑如何求出比原方案字典序小的方案总数
那
么
对
于
a
[
i
]
位
置
上
的
字
母
,
我
们
假
设
前
i
−
1
个
字
母
都
相
同
(
这
样
就
可
以
保
证
不
重
复
)
,
我
们
发
现
,
后
面
比
a
[
i
]
字
典
序
小
的
字
母
放
在
a
[
i
]
位
置
一
定
可
以
比
当
前
放
置
方
案
字
典
序
小
,
然
后
后
面
其
它
字
母
随
便
放
就
行
了
那么对于a[i]位置上的字母,我们假设前i-1个字母都相同(这样就可以保证不重复),我们发现,后面比a[i]字典序小的字母放在a[i]位置一定可以比当前放置方案字典序小,然后后面其它字母随便放就行了
那么对于a[i]位置上的字母,我们假设前i−1个字母都相同(这样就可以保证不重复),我们发现,后面比a[i]字典序小的字母放在a[i]位置一定可以比当前放置方案字典序小,然后后面其它字母随便放就行了
所
以
总
答
案
为
∑
i
=
1
n
∑
j
=
i
n
[
a
[
j
]
<
a
[
i
]
]
所以总答案为\sum_{i=1}^{n}\sum_{j=i}^{n}[ a[j]<a[i]]
所以总答案为∑i=1n∑j=in[a[j]<a[i]] 我们发现这个逆序对个数是可以用树状数组维护,所以时间复杂度为(nlogn)
#include<bits/stdc++.h>
using namespace std;
const int N=1e6+5;
typedef long long ll;
ll ansout=1,n,mod=998244353;
ll jc[N];
int a[N],tree[N];
ll lowbit(ll x){
return (x&(-x));
}
void add(int val){for(int rt=val;rt<=n;rt+=lowbit(rt)) (++tree[rt])%mod; }
ll sum(int val){
ll ans=0;
for(int rt=val;rt;rt-=lowbit(rt)) ans=(tree[rt]+ans)%mod;
return ans;
}
int main(){
scanf("%d",&n);
jc[0]=1;
for(int i=1;i<=n;i++){
jc[i]=1ll*jc[i-1]*i%mod;
scanf("%d",&a[i]);
}
add(a[n]);
for(int i=n-1;i>=1;i--){
ll sumtmp=sum(a[i])%mod;
ansout=(ansout+(sumtmp*jc[n-i])%mod)%mod;
add(a[i]);
}
printf("%lld",ansout);
}