結果

問題 No.696 square1001 and Permutation 5
ユーザー PulmnPulmn
提出日時 2018-05-29 23:12:30
言語 C++11
(gcc 13.3.0)
結果
AC  
実行時間 7,909 ms / 10,000 ms
コード長 3,369 bytes
コンパイル時間 2,094 ms
コンパイル使用メモリ 178,292 KB
実行使用メモリ 37,120 KB
最終ジャッジ日時 2024-06-30 08:06:06
合計ジャッジ時間 37,361 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define syosu(x) fixed<<setprecision(x)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> P;
typedef pair<double,double> pdd;
typedef pair<ll,ll> pll;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<double> vd;
typedef vector<vd> vvd;
typedef vector<ll> vl;
typedef vector<vl> vvl;
typedef vector<string> vs;
typedef vector<P> vp;
typedef vector<vp> vvp;
typedef vector<pll> vpll;
typedef pair<int,P> pip;
typedef vector<pip> vip;
const int inf=1<<30;
const ll INF=1ll<<60;
const double pi=acos(-1);
const double eps=1e-8;
const ll mod=1e9+7;
const int dx[4]={-1,0,1,0},dy[4]={0,-1,0,1};

const ll p=924844033,r=5;

ll Pow(ll n,ll m){
	ll res=1;
	while(m>0){
		if(m&1) (res*=n)%=p;
		(n*=n)%=p;
		m>>=1;
	}
	return res;
}

ll Inv(ll x){return Pow(x,p-2);}

void ntt(vl& a,bool B){
	int n=a.size(),h=0;
	while(1<<h<n) h++;
	for(int i=0;i<n;i++){
		int j=0;
		for(int k=0;k<h;k++) if(i&1<<k) j+=1<<(h-k-1);
		if(i<j) swap(a[i],a[j]);
	}
	for(int i=1;i<n;i*=2){
		ll w=Pow(r,(p-1)/(i*2));
		if(B) w=Inv(w);
		for(int j=0;j<n;j+=i*2){
			ll wn=1;
			for(int k=0;k<i;k++){
				ll s=a[j+k],t=a[i+j+k]*wn%p;
				a[j+k]=(s+t)%p;
				a[i+j+k]=(s-t+p)%p;
				(wn*=w)%=p;
			}
		}
	}
	if(B){
		ll v=Inv(n);
		for(int i=0;i<n;i++) (a[i]*=v)%=p;
	}
}

vl Conv(vi a_,vi b_){
	int A=a_.size(),B=b_.size();
	vl a(A),b(B);
	for(int i=0;i<A;i++) a[i]=a_[i];
	for(int i=0;i<B;i++) b[i]=b_[i];
	int s=A+B-1,t=1;
	while(t<s) t*=2;
	a.resize(t);
	b.resize(t);
	ntt(a,0);ntt(b,0);
	for(int i=0;i<t;i++) (a[i]*=b[i])%=p;
	ntt(a,1);
	a.resize(s);
	return a;
}

class Segment_Tree{
	private:
	int n;
	vi date;
	int Rec(int a,int b,int k,int l,int r){
		if(r<=a||b<=l) return 0;
		if(a<=l&&r<=b) return date[k];
		int m=(l+r)/2;
		return Rec(a,b,k*2+1,l,m)+Rec(a,b,k*2+2,m,r);
	}
	public:
	Segment_Tree(int n_){
		n=1;
		while(n<n_) n*=2;
		date=vi(2*n-1);
	}
	void Update(int k,int x){
		k+=n-1;
		date[k]+=x;
		while(k>0){
			k=(k-1)/2;
			date[k]=date[k*2+1]+date[k*2+2];
		}
	}
	int Query(int a,int b){
		return Rec(a,b,0,0,n);
	}
	int Open(int k){
		return date[k+n-1];
	}
};

vi Sum(vi a,vi b){
	int A=a.size(),B=b.size(),N=max(A,B),t=0;
	vi c(N+1);
	for(int i=0;i<N;i++){
		t+=(i<A?a[i]:0)+(i<B?b[i]:0);
		c[i]=t%10;
		t/=10;
	}
	c[N]=t;
	while(!c.empty()&&!c.back()) c.pop_back();
	return c;
}

vi Mul(vi a,vi b){
	vl c=Conv(a,b);
	int N=c.size(),t=0;
	vi d(N);
	for(int i=0;i<N;i++){
		t+=c[i];
		d[i]=t%10;
		t/=10;
	}
	while(t){
		d.push_back(t%10);
		t/=10;
	}
	while(!d.empty()&&!d.back()) d.pop_back();
	return d;
}

vi Int(int t){
	vi a;
	while(t){
		a.push_back(t%10);
		t/=10;
	}
	return a;
}

int n;
vi a;
vvi b,c;

int main(){
	ios::sync_with_stdio(0);
	cin.tie(0);
	cin>>n;
	a=vi(n);
	b=vvi(n);
	c=vvi(n-1);
	Segment_Tree st(n);
	for(int i=0;i<n;i++) cin>>a[i];
	for(int i=0;i<n;i++){
		int A=a[n-i-1]-1;
		st.Update(A,1);
		b[i]=Int(st.Query(0,A));
	}
	for(int i=1;i<n;i++) c[i-1]=Int(i);
	while(b.size()>1){
		int N=b.size(),M=(N+1)/2;
		vvi d(M);
		for(int i=0;i<N/2;i++) d[i]=Sum(Mul(b[2*i+1],c[2*i]),b[2*i]);
		if(N%2==1) d[M-1]=b.back();
		b=d;
		N=c.size(),M=(N+1)/2;
		d=vvi(M);
		for(int i=0;i<N/2;i++) d[i]=Mul(c[2*i],c[2*i+1]);
		if(N%2==1) d[M-1]=c.back();
		c=d;
	}
	vi res=b[0];
	res=Sum(res,vi(1,1));
	reverse(res.begin(),res.end());
	for(auto i:res) cout<<i;
	cout<<endl;
}
0