結果

問題 No.767 配られたジャパリまん
コンテスト
ユーザー chocorusk
提出日時 2018-12-15 01:16:30
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 260 ms / 2,000 ms
コード長 1,936 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 847 ms
コンパイル使用メモリ 99,636 KB
実行使用メモリ 14,976 KB
最終ジャッジ日時 2024-09-25 05:39:06
合計ジャッジ時間 2,632 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 6
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;
const ll MOD=1e8+7;
ll powmod(ll a, ll k){
    ll ap=a, ans=1;
    while(k>0){
        if(k%2==1){
            ans*=ap;
            ans%=MOD;
        }
        ap=ap*ap;
        ap%=MOD;
        k/=2;
    }
    return ans;
}
ll inv(ll a){
	return powmod(a, MOD-2);
}
int main()
{
	ll f[200001]; f[0]=1;
	for(ll i=1; i<=200000; i++) f[i]=f[i-1]*i%MOD;
	ll invf[200001];
	invf[200000]=inv(f[200000]);
	for(ll i=200000-1; i>=0; i--) invf[i]=invf[i+1]*(i+1)%MOD;
	int h, w, k;
	cin>>h>>w>>k;
	int a[20], b[20]; P q[20]; int ind[20];
	for(int i=0; i<k; i++) cin>>a[i]>>b[i], q[i]=P(a[i], b[i]), ind[i]=i;
	sort(ind, ind+k, [&](int i, int j){ return q[i]<q[j];});
	ll c0, c[1<<20];
	c0=f[h+w]*invf[h]%MOD*invf[w]%MOD;
	c[0]=c0;
	for(int i=1; i<(1<<k); i++){
		P pr=P(0, 0);
		c0=1;
		int i0=0;
		for(int j=0; j<k; j++){
			if((i&(1<<j))==0) continue;
			i0|=(1<<(ind[j]));
			if(c0==0) continue;
			P p=q[ind[j]];
			if(p.first>=pr.first && p.second>=pr.second){
				c0*=(f[p.first-pr.first+p.second-pr.second]*invf[p.first-pr.first]%MOD*invf[p.second-pr.second]%MOD);
				c0%=MOD;
				pr=p;
			}else{
				c0=0;
			}
		}
		if(c0!=0){
		    c0=(c0*f[h-pr.first+w-pr.second]%MOD*invf[h-pr.first]%MOD*invf[w-pr.second]%MOD);
		}
		c[i0]=c0;
	}
	for(int i=0; i<k; i++){
		for(int j=0; j<(1<<k); j++){
			if(j&(1<<i)){
			    c[j]-=c[j^(1<<i)];
			    c[j]=(c[j]+MOD)%MOD;
			}
		}
	}
	for(int i=0; i<(1<<k); i++){
	    int ct=0;
	    for(int j=0; j<k; j++) if(i&(1<<j)) ct++;
	    if(ct&1) printf("%lld\n", (MOD-c[i])%MOD);
	    else printf("%lld\n", c[i]);
	}
	return 0;
}
0