正題

---

題目大意

求n個(gè)數(shù)的全排列的第k個(gè)。

---

解題思路

首先康拓逆展開
${\sum }_i^{i<=n}{x}_i(n?i)!$
求每個(gè)時(shí)候第$x$大的數(shù)
然后因?yàn)?span id="ozvdkddzhkzd" class="katex--inline">$n(n?1)!=n!$
so我們可以直接用余數(shù)

這是$n=3$時(shí)是序列，我們可以發(fā)現(xiàn)我們需要求第$x_i$個(gè)。這時(shí)候我們可以用二分加樹狀數(shù)組來(lái)$lo{g}^2$的求答案。
之后我們又會(huì)發(fā)現(xiàn)
${\sum }_i^{i<=n}{x}_i(n?i)!$
這個(gè)直接計(jì)算時(shí)間復(fù)雜度是$n^{2}lo{g}^{2}nw$(w是高進(jìn)度位數(shù))。
這時(shí)我們引入:
$x\%yz/z=x/z\%y$
證明:

這時(shí)候我們可以枚舉$1～n$，然后每次$kmodi$，這時(shí)候的余數(shù)就是上面式子$n?i+1$次的余數(shù)。

---

code

#include<iostream> #include<cstring> #define lobit(x) x&-x #define N 100010 #define ll long long using namespace std; ll a[N*2],n,t[N],l,mo[N]; char k[N*2]; void change(ll x,ll num)//改變 {while(x<=n){t[x]+=num;x+=lobit(x);} } ll ask(ll x)//詢問 {ll sum=0;while(x){sum+=t[x];x-=lobit(x);}return sum; } void read()//輸入——高精度 {scanf("%s",k);l=strlen(k)-1;ll L=0;for(ll i=l;i>=0;i--) L+=((l-i)%13)==0,mo[i]=L-1;//壓行之后的位置for(ll i=0;i<=l;i++) a[mo[i]]=a[mo[i]]*10+k[i]-48;//求值a[0]--;for(ll i=0;i<=l;i++) if(a[i]<0) a[i+1]--,a[i]+=(1e13);//要先-1l=L-1; } ll div(ll x)//高精除 { ll g=0;for (ll i=l;i>=0;i--){ll s;s=g*(1e13)+a[i];a[i]=s/x;g=s%x;}while(!a[l]) l--;return g; } int main() {scanf("%lld",&n);read();for(ll i=1;i<=n;i++)mo[n-i+1]=div(i);//計(jì)算余數(shù)for(ll i=1;i<=n;i++) {ll ans=0;for(ll l=1,r=n;l<=r;) {ll m=(l+r)/2;ll s=0,x=m;s=x-ask(x);if(s<=mo[i]) l=m+1;else ans=m,r=m-1;}//二分位置ll x=ans;change(x,1);//去掉這個(gè)數(shù)printf("%lld ",ans);}return 0; }

總結(jié)

以上是生活随笔為你收集整理的nssl1156-今天你AK了吗?【康托展开,高精度,二分答案,树状数组】的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。