結果

問題 No.1846 Good Binary Matrix
コンテスト
ユーザー 👑 potato167
提出日時 2021-08-21 17:05:45
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 635 ms / 2,000 ms
コード長 1,396 bytes
コンパイル時間 2,106 ms
コンパイル使用メモリ 198,104 KB
最終ジャッジ日時 2025-01-24 01:23:41
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 35
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using std::cout;
using std::cin;
using std::endl;
using ll=long long;
ll mod=1e9+7;
#define rep(i,a) for (ll i=0;i<a;i++)

//return x^y (mod z)
ll my_pow(ll x,ll y,ll z){
	if(x==0) return 0;
	if(y==0) return 1;
	if(y%2==1) return (x*my_pow(x,y-1,z))%z;
	return my_pow((x*x)%z,y/2,z);
}

namespace po167{
struct combination{
	std::vector<long long> fact;
	std::vector<long long> rev;
	std::vector<long long> fact_rev;
	long long MOD;
	combination(int max,long long mod):MOD(mod),fact(max+1),rev(max+1),fact_rev(max+1){
		for(long long i=0;i<=max;i++){
			if(i<2){
				fact[i]=1;
				fact_rev[i]=1;
				rev[i]=1;
				continue;
			}
			fact[i]=(fact[i-1]*i)%mod;
			rev[i]=mod-((mod/i)*rev[mod%i])%mod;
			fact_rev[i]=(fact_rev[i-1]*rev[i])%mod;
		}
	}
	long long Comb(int x,int y){
		if (x<y||y<0) return 0;
		return (((fact_rev[y]*fact_rev[x-y])%MOD)*fact[x])%MOD;
	}
};
}
using po167::combination;

void solve(int H,int W);

//おちこんだりもしたけれど、私はげんきです。
int main() {

	int H,W;
	cin>>H>>W;
	assert(1<=H&&H<=(int)(1e6));
	assert(1<=W&&W<=(int)(1e6));
	solve(H,W);
}

void solve(int H,int W){
	combination C(H,mod);
	ll ans=0;
	ll tmp=1;
	ll sign=1;
	if((H+W)%2) sign=-1;
	rep(i,H+1){
		ans+=sign*(C.Comb(H,i)*my_pow(1-tmp,W,mod))%mod;
		ans%=mod;
		tmp*=2;
		tmp%=mod;
		sign*=-1;
	}
	cout<<(ans+mod)%mod<<"\n";
}
0