結果
| 問題 |
No.767 配られたジャパリまん
|
| コンテスト | |
| ユーザー |
chocorusk
|
| 提出日時 | 2018-12-15 00:31:26 |
| 言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
| 結果 |
AC
|
| 実行時間 | 1,703 ms / 2,000 ms |
| コード長 | 1,728 bytes |
| コンパイル時間 | 1,138 ms |
| コンパイル使用メモリ | 103,564 KB |
| 実行使用メモリ | 14,848 KB |
| 最終ジャッジ日時 | 2024-09-25 05:32:27 |
| 合計ジャッジ時間 | 4,691 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 6 |
| other | AC * 15 |
ソースコード
#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];
for(int i=0; i<k; i++) cin>>a[i]>>b[i];
ll c[1<<20];
c[0]=f[h+w]*invf[h]%MOD*invf[w]%MOD;
for(int i=1; i<(1<<k); i++){
vector<P> v;
for(int j=0; j<k; j++){
if(i&(1<<j)) v.push_back({a[j], b[j]});
}
sort(v.begin(), v.end());
v.push_back({h, w});
P pr=P(0, 0);
c[i]=1;
for(auto p:v){
if(p.first>=pr.first && p.second>=pr.second){
c[i]*=(f[p.first-pr.first+p.second-pr.second]*invf[p.first-pr.first]%MOD*invf[p.second-pr.second]%MOD);
c[i]%=MOD;
pr=p;
}else{
c[i]=0;
break;
}
}
}
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) cout<<(MOD-c[i])%MOD<<endl;
else cout<<c[i]<<endl;
}
return 0;
}
chocorusk