結果

問題 No.696 square1001 and Permutation 5
コンテスト
ユーザー Pulmn
提出日時 2018-05-30 19:32:41
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 2,888 bytes
コンパイル時間 1,824 ms
コンパイル使用メモリ 180,940 KB
実行使用メモリ 21,828 KB
最終ジャッジ日時 2024-06-30 08:20:10
合計ジャッジ時間 11,650 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 8 WA * 4
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<ll> vl;

const ll p=1012924417,r=5,mod=1000;

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);
	}
};

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%mod;
		t/=mod;
	}
	c[N]=t;
	while(!c.empty()&&!c.back()) c.pop_back();
	return c;
}

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

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

void Out(vi a){
	reverse(a.begin(),a.end());
	for(int i=0;i<a.size();i++){
		string s=to_string(a[i]);
		while(i&&s.size()<3) s='0'+s;
		cout<<s;
	}
	cout<<endl;
}

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[i*2+1],c[i*2]),b[i*2]);
		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[i*2],c[i*2+1]);
		if(N%2==1) d[M-1]=c.back();
		c=d;
	}
	Out(Sum(b[0],vi(1,1)));
}
0