結果

問題 No.907 Continuous Kadomatu
コンテスト
ユーザー HIR180
提出日時 2019-10-11 22:06:20
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 234 ms / 2,000 ms
コード長 4,009 bytes
コンパイル時間 2,743 ms
コンパイル使用メモリ 189,688 KB
実行使用メモリ 69,504 KB
最終ジャッジ日時 2024-11-25 07:33:29
合計ジャッジ時間 8,174 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 5
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
using namespace std;
typedef long long ll;
typedef pair<int,int> P;
typedef pair<int,P> P1;
typedef pair<P,P> P2;
#define pu push
#define pb push_back
#define mp make_pair
#define eps 1e-7
#define INF 1000000000
#define mod 1000000007
#define fi first
#define sc second
#define rep(i,x) for(int i=0;i<x;i++)
#define repn(i,x) for(int i=1;i<=x;i++)
#define SORT(x) sort(x.begin(),x.end())
#define ERASE(x) x.erase(unique(x.begin(),x.end()),x.end())
#define POSL(x,v) (lower_bound(x.begin(),x.end(),v)-x.begin())
#define POSU(x,v) (upper_bound(x.begin(),x.end(),v)-x.begin())
#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
using namespace std;
typedef long long ll;
typedef pair<int,int> P;
typedef pair<int,P> P1;
typedef pair<P,P> P2;
#define pu push
#define pb push_back
#define mp make_pair
#define eps 1e-7
#define INF 1000000000
#define mod 1000000007
#define fi first
#define sc second
#define rep(i,x) for(int i=0;i<x;i++)
#define repn(i,x) for(int i=1;i<=x;i++)
#define SORT(x) sort(x.begin(),x.end())
#define ERASE(x) x.erase(unique(x.begin(),x.end()),x.end())
#define POSL(x,v) (lower_bound(x.begin(),x.end(),v)-x.begin())
#define POSU(x,v) (upper_bound(x.begin(),x.end(),v)-x.begin())
int n;
vector<ll>za;
ll dp[205][405][205],a[205],b[205];
ll modpow(ll x,ll n){
	ll res=1;
	while(n>0){
		if(n&1) res=res*x%mod;
		x=x*x%mod;
		n>>=1;
	}
	return res;
}
ll F[505],R[505];
void make(){
	F[0] = 1;
	for(int i=1;i<505;i++) F[i] = F[i-1]*i%mod;
	for(int i=0;i<505;i++) R[i] = modpow(F[i],mod-2);
}
ll C(int a,int b){
	return F[a]*R[b]%mod*R[a-b]%mod;
}
ll coef[205];
//k = n+1-2*m
//sum_{m=0...n/2} (-1)^m*2^{1-k}*sum_{j=0...k} binomial(k,j)*(-1)^j*(k-2*j)^(n+1)/k - [n=1]
void make2(){
	coef[1] = 1;
	for(int n=2;n<205;n++){
		for(int m=0;m<=n/2;m++){
			ll x = (m%2==0?1:-1);
			int k = n+1-2*m;
			x *= 2LL;
			x = x%mod*modpow(modpow(2,k),mod-2)%mod;
			ll S = 0;
			for(int j=0;j<=k;j++){
				ll c = C(k,j);
				if(j%2 == 1) c *= -1;
				c = c*modpow((k-j-j),n+1)%mod;
				c = c*modpow(k,mod-2)%mod;
				S += c;
			}
			S %= mod;
			coef[n] += x*S%mod;
		}
		coef[n] = (coef[n]%mod+mod)%mod;
	}
	for(int i=1;i<205;i++){
		coef[i] = coef[i] * R[i] % mod;
		if(i > 1) coef[i] = coef[i] * modpow(2LL,mod-2) % mod;
	}
	//for(int i=1;i<=10;i++) cout << coef[i] << endl;
}
ll sum[205][405],rui[205][405];
int main(){
	make();
	make2();
	
	cin >> n;
	rep(i,n) cin >> a[i] >> b[i];
	rep(i,n) { za.pb(a[i]); za.pb(b[i]); }
	SORT(za); ERASE(za);
	int L = POSL(za,a[0]), R = POSL(za,b[0]);
	for(int i=L;i<R;i++){
		ll up = za[i+1]-za[i];
		up = up * modpow(za[R]-za[L],mod-2) % mod;
		
		dp[0][i][1] = up;
	}
	for(int j=0;j<za.size();j++){
		for(int k=1;k<205;k++){
			sum[0][j] += dp[0][j][k]*coef[k]%mod;
		}
		sum[0][j] %= mod;
		rui[0][j] += sum[0][j];
		if(j) rui[0][j] += rui[0][j-1];
		rui[0][j] %= mod;
	}
	for(int i=1;i<n;i++){
		int L = POSL(za,a[i]);
		int R = POSL(za,b[i]);
		ll rev = modpow(b[i]-a[i],mod-2);
		for(int j=L;j<R;j++){
			ll up = za[j+1]-za[j];
			up = up*rev%mod;
			for(int k=1;k<=i+1;k++){
				//go to same
				if(k > 1){
					dp[i][j][k] += dp[i-1][j][k-1]*up%mod;
				}
				//go to dif
				if(k == 1){
    				if(i%2 == 1){
    					ll x = (j==0?0:rui[i-1][j-1]);
    					dp[i][j][k] += x%mod*up%mod;
    				}
    				else{
    					ll x = rui[i-1][za.size()-1];
    					x -= rui[i-1][j];
    					dp[i][j][k] += x%mod*up%mod;
    				}
				}
				if(dp[i][j][k] >= mod) dp[i][j][k] -= mod;
				if(dp[i][j][k] < 0) dp[i][j][k] += mod;
			}
		}
		for(int j=0;j<za.size();j++){
			for(int k=1;k<205;k++){
				sum[i][j] += dp[i][j][k]*coef[k]%mod;
			}
			sum[i][j] %= mod;
			rui[i][j] += sum[i][j];
			if(j) rui[i][j] += rui[i][j-1];
			rui[i][j] %= mod;
		}
	}
	ll ans = 0;
	for(int j=0;j<za.size();j++){
		for(int k=1;k<205;k++){
			ans += dp[n-1][j][k] * coef[k] % mod;
		}
	}
	cout << (ans%mod+mod)%mod << endl;
}
0