結果

問題 No.767 配られたジャパリまん
コンテスト
ユーザー tails
提出日時 2019-01-23 18:55:54
言語 cLay
(20241019-1)
結果
AC  
実行時間 85 ms / 2,000 ms
コード長 585 bytes
コンパイル時間 6,021 ms
コンパイル使用メモリ 224,740 KB
実行使用メモリ 23,040 KB
最終ジャッジ日時 2024-07-05 13:18:31
合計ジャッジ時間 7,556 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 6
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define MD 100000007

int h,w,k,a[20],b[20],c[20],f[1<<20],g[1<<20];
mint d[1<<20],e[1<<20],l;

{
	combination_mint cm;
	cm.init(3d5);
	
	rd(h,w,k,(a,b)(k));
	rep(i,k)c[i]=i;
	sortA(k,a,b,c);

	e[0]=1;
	d[0]=l=cm.C(h+w,h);
	rep(i,1,1<<k){
		int j0=ffs(i)-1;
		int j1=ffs(i&i-1)-1;
		int y=(j1<0?h:a[j1])-a[j0];
		int x=(j1<0?w:b[j1])-b[j0];
		e[i]=e[i&i-1]*cm.C(y+x,y);
		d[i]=e[i]*cm.C(a[j0]+b[j0],a[j0]);
	}
	d[0]=0;
	MoebiusTransform(1<<k,d);
	rep(i,1,1<<k)f[i]=f[i&i-1]|1<<c[ffs(i)-1];
	rep(i,1<<k)g[i]=l+(__builtin_parity(i)?-d[i]:d[i]);
	sortA(1<<k,f,g);
	rep(i,1<<k)wt(g[i]);
}
0