CF 1443 E. Long Permutation 性质分析+逐段处理

很经典的思维套路。

我们知道：长度为i的排列共有 i！ 种，显然20！> 1e5*2e5  

所以其实无论怎么改变，最多修改不到20个数。

我们可以类似于二进制的思想，从前往后修改。

比如1 2 3 4 5，

对于2号位置，其于其后面的4个数总共有4！ = 24种排列，

那么对于1号位置来说：其排列为24*5=120，即这五个数每个数放在最前面，后面都是4个数的排列。

对于字典序从低到高来说：

我们把2放在最前面，其余相对位置不变往后移变为：21345，这个排列是第25字典序小的排列（类比二进制思想，1在首尾时，字典序均小于2在首尾，共24种排列）

对于1号位最多可以控制24*4=96个排列。

然后我们到2号位继续上述操作。

这个过程类比这二进制自己想一下很容易就想出来了。

然后就是实现这个。。预处理阶乘和前缀和即可。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define re register
#define ls (o<<1)
#define rs (o<<1|1)
//#define m (l+r)/2
#define pb push_back
typedef pair<int,int> pii;
const double PI= acos(-1.0);
const int M = 2e5+7;
/*
int head[M],cnt=1;
void init(int n){cnt=1;for(int i=0;i<=n;i++)head[i]=0;}
struct EDGE{int to,nxt,w;}ee[M*2];
void add(int x,int y,int w){ee[++cnt].nxt=head[x],ee[cnt].w=w,ee[cnt].to=y,head[x]=cnt;}
*/
int a[M];
int n,q;
ll sm[M];
ll fac[M];
void gao(int L,ll nm){
	if(!nm||L>n)return ;
/*	cout<<"----------------------------------"<<endl;
	for(int i=1;i<=n;i++)cout<<a[i]<<" ";cout<<endl;
	cout<<" *--  "<<fac[L+1]<<endl;*/
	int tp=0;
	for(int i=L+1;i<=n;i++){
		if((ll)(i-L)*fac[L+1]<=nm)
			tp++;
		else break;
	}
	nm-=tp*fac[L+1];
	int nw=a[L+tp];
	for(int i=L+tp-1;i>=L;i--)
		a[i+1]=a[i];
	a[L]=nw;
	/*cout<<"**********   "<<tp<<"     "<<L<<"  "<<nm<<" "<<endl;
	for(int i=1;i<=n;i++)cout<<a[i]<<" ";
	cout<<endl;
	cout<<"----------------------------------"<<endl;*/
	gao(L+1,nm);
}
ll cal(int l,int r){
	return sm[r]-sm[l-1];
}
int main()
{
	ios::sync_with_stdio(false);
  	cin.tie(0);
	cin>>n>>q;
	for(int i=1;i<=n;i++)a[i]=i,sm[i]=sm[i-1]+a[i];
	fac[n]=1;
	ll up=2e10+7;
	for(int i=n-1;i>=1;i--){
		if(fac[n]>=up)break;
		fac[i]=fac[i+1]*(n-i+1);
	}
	fac[0]=1e18;
	ll nm=0;//总共需要变换多少次 
	for(int i=1;i<=q;i++) {
		int op,l,r,x;
		cin>>op;
		if(op==1){
			cin>>l>>r;
			ll ans=0;
			if(l<=n-20)ans+=cal(l,min(r,n-20));
			for(int j=max(n-19,l);j<=r;j++){
				ans+=a[j];
			}
			cout<<ans<<endl;
		}else{
			if(n>20)for(int j=n-20;j<=n;j++)a[j]=j;
			else for(int j=1;j<=n;j++)a[j]=j;
			cin>>x;
			nm+=x;
			ll ans=0;
			for(int j=n;j>=max(1,n-20);j--){
				if(fac[j]>nm){
					gao(j,nm);
					break;
				}
			}
		}
	}
	return 0;
}