結果

問題 No.696 square1001 and Permutation 5
ユーザー PulmnPulmn
提出日時 2018-05-29 23:15:26
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 8,021 ms / 10,000 ms
コード長 2,790 bytes
コンパイル時間 1,886 ms
コンパイル使用メモリ 178,184 KB
実行使用メモリ 37,208 KB
最終ジャッジ日時 2024-06-30 08:06:56
合計ジャッジ時間 37,700 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 8,021 ms
36,468 KB
testcase_01 AC 3 ms
5,376 KB
testcase_02 AC 4 ms
5,376 KB
testcase_03 AC 10 ms
5,376 KB
testcase_04 AC 22 ms
5,376 KB
testcase_05 AC 48 ms
5,376 KB
testcase_06 AC 137 ms
5,376 KB
testcase_07 AC 337 ms
6,460 KB
testcase_08 AC 990 ms
10,008 KB
testcase_09 AC 3,419 ms
19,916 KB
testcase_10 AC 7,980 ms
37,208 KB
testcase_11 AC 5,962 ms
31,748 KB
testcase_12 AC 2 ms
5,376 KB
testcase_13 AC 2 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

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

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