結果

問題 No.1001 注文の多い順列
コンテスト
ユーザー vjudge1
提出日時 2026-10-05 15:27:21
言語 C++17(gcc12)
(gcc 12.4.0 + boost 1.92.0 + ACL)
コンパイル:
g++-12 -O2 -lm -std=c++17 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 64 ms / 2,000 ms
+ 375µs
コード長 1,235 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,615 ms
コンパイル使用メモリ 203,080 KB
実行使用メモリ 49,792 KB
最終ジャッジ日時 2026-10-05 15:27:27
合計ジャッジ時間 4,990 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
using namespace std;
#define int long long
#define pi pair<int,int>
#define ull unsigned long long
#define add(a,b) (a+=b,a=a>=mod?a-mod:a)
const int mod=1e9+7;
const int N=3010;
int cnt0[N],cnt1[N];
int dp[N][N];
int fct[N],inv[N];
int fp(int a,int b){
	int ans=1;
	while(b){
		if(1&b)ans=ans*a%mod;
		a=a*a%mod;
		b>>=1;
	}
	return ans;
}
int C(int a,int b){return (a<b||b<0)?0:fct[a]*inv[a-b]%mod;}
signed main(){
	ios::sync_with_stdio(false);
	cin.tie(0),cout.tie(0);
//	freopen("mod.in","r",stdin);
//	freopen("ans2.out","w",stdout);
	fct[0]=1;for(int i=1;i<N;i++)fct[i]=fct[i-1]*i%mod;
	inv[N-1]=fp(fct[N-1],mod-2)%mod;
	for(int i=N-1;i>=1;i--)inv[i-1]=inv[i]*i%mod;
	inv[0]=1;
	int n;cin>>n;int s0=0,s1=0;
	for(int i=1;i<=n;i++){
		int t,x;cin>>t>>x;
		if(t)cnt1[x]++,s1++;
		else cnt0[x]++,s0++;
	}
	for(int i=1;i<=n;i++)cnt1[i]+=cnt1[i-1];
	for(int i=1;i<=n;i++)cnt0[i]+=cnt0[i-1];
	dp[0][0]=1;
	for(int i=0;i<n;i++){
		for(int c0=0;c0<=i;c0++){
			int c1=i-c0;if(!dp[i][c1])continue;
			if(cnt1[i+1]-c1>0)
				add(dp[i+1][c1+1],dp[i][c1]*(cnt1[i+1]-c1)%mod*C(c0-cnt0[i],cnt0[i+1]-cnt0[i])%mod);
			add(dp[i+1][c1],dp[i][c1]*C(c0-cnt0[i]+1,cnt0[i+1]-cnt0[i])%mod);
		}
	}
	cout<<dp[n][s1];
	return 0;
}
0