結果

問題 No.767 配られたジャパリまん
ユーザー Hello-World-hello-worldHello-World-hello-world
提出日時 2024-10-17 19:49:42
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 244 ms / 2,000 ms
コード長 2,181 bytes
コンパイル時間 1,840 ms
コンパイル使用メモリ 173,196 KB
実行使用メモリ 24,716 KB
最終ジャッジ日時 2024-10-17 19:49:46
合計ジャッジ時間 3,835 ms
ジャッジサーバーID
(参考情報)
judge2 / judge3
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 7 ms
10,444 KB
testcase_01 AC 7 ms
11,716 KB
testcase_02 AC 7 ms
10,792 KB
testcase_03 AC 7 ms
11,048 KB
testcase_04 AC 7 ms
11,340 KB
testcase_05 AC 21 ms
12,104 KB
testcase_06 AC 7 ms
11,756 KB
testcase_07 AC 6 ms
10,372 KB
testcase_08 AC 7 ms
10,856 KB
testcase_09 AC 7 ms
11,768 KB
testcase_10 AC 13 ms
10,820 KB
testcase_11 AC 8 ms
11,224 KB
testcase_12 AC 7 ms
10,612 KB
testcase_13 AC 14 ms
10,992 KB
testcase_14 AC 6 ms
11,312 KB
testcase_15 AC 6 ms
11,108 KB
testcase_16 AC 8 ms
11,452 KB
testcase_17 AC 7 ms
11,308 KB
testcase_18 AC 7 ms
10,908 KB
testcase_19 AC 8 ms
11,824 KB
testcase_20 AC 244 ms
24,716 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

//Date: 2024-10-17 17:57:39
#include<bits/stdc++.h>
using namespace std;
#define int long long
inline int rd() {
    int s=0,m=0;char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-')m=1;ch=getchar();}
    while( isdigit(ch)) s=(s<<3)+(s<<1)+(ch^48),ch=getchar();
    return m?-s:s;
}
bool MBE;
namespace SOLVER {
const int M=1e8+7;
int fr[200005],ifr[200005],inv[200005];int C(int n,int m) {return n<m||m<0?0:fr[n]*ifr[m]%M*ifr[n-m]%M;}
int H,W,n,x[25],y[25],id[25],f[(1<<20)+5],ans[(1<<20)+5];bool cmp(int i,int j) {return x[i]+y[i]<x[j]+y[j];}
void OR(int f[],int op) {
    for(int i=1;i<(1<<n);i<<=1) for(int j=0;j<(1<<n);j+=(i<<1)) 
        for(int k=0;k<i;k++) f[i+j+k]=(f[i+j+k]+f[j+k]*op+M)%M;
}
void AND(int f[],int op) {
    for(int i=1;i<(1<<n);i<<=1) for(int j=0;j<(1<<n);j+=(i<<1)) 
        for(int k=0;k<i;k++) f[j+k]=(f[j+k]+f[i+j+k]*op+M)%M;
}
void XOR(int f[],int op) {
    for(int i=1;i<(1<<n);i<<=1) for(int j=0;j<(1<<n);j+=(i<<1)) 
        for(int k=0;k<i;k++) {
            int X=f[j+k],Y=f[i+j+k];
            f[j+k]=((X+Y)*op%M+M)%M,f[i+j+k]=((X-Y)*op%M+M)%M;
        }
}
void MAIN() {
    inv[1]=1;fr[0]=ifr[0]=1;
    for(int i=2;i<=200000;i++) inv[i]=inv[i]=(M-M/i*inv[M%i]%M)%M;
    for(int i=1;i<=200000;i++) fr[i]=fr[i-1]*i%M,ifr[i]=ifr[i-1]*inv[i]%M;
    cin>>H>>W>>n;for(int i=0;i<n;i++) x[i]=rd(),y[i]=rd(),id[i]=i;
    sort(id,id+n,cmp);f[0]=C(H+W,H);for(int i=0;i<n;i++) f[1<<i]=C(x[id[i]]+y[id[i]],x[id[i]]);
    for(int i=1;i<(1<<n);i++) if(__builtin_popcount(i)!=1) {
        int t=__lg(i),s=__lg(i-(1<<t));
        f[i]=f[i-(1<<t)]*C(x[id[t]]-x[id[s]]+y[id[t]]-y[id[s]],x[id[t]]-x[id[s]])%M;
    }
    for(int i=1;i<(1<<n);i++) {int s=__lg(i);f[i]=f[i]*C(H-x[id[s]]+W-y[id[s]],H-x[id[s]])%M;}
    AND(f,-1);OR(f,1);
    for(int i=0;i<(1<<n);i++) {
        int sum=0;
        for(int j=0;j<n;j++) if(i>>j&1) sum|=(1<<id[j]);
        ans[sum]=f[i];
    } 
    for(int i=0;i<(1<<n);i++) printf("%lld\n",ans[(1<<n)-i-1]);
}
}
bool MED;
signed main() {
    //freopen(".in","r",stdin);freopen(".out","w",stdout);
    for(int tt=1;tt;tt--) SOLVER::MAIN();
    // cerr<<(&MBE-&MED)/1024<<" KB, "<<1000*clock()/CLOCKS_PER_SEC<<" ms\n";
    return 0;
}
0