結果
問題 | No.1846 Good Binary Matrix |
ユーザー |
|
提出日時 | 2023-04-10 21:57:09 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 507 ms / 2,000 ms |
コード長 | 811 bytes |
コンパイル時間 | 1,892 ms |
コンパイル使用メモリ | 193,044 KB |
最終ジャッジ日時 | 2025-02-12 04:51:43 |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 35 |
ソースコード
#include<bits/stdc++.h> using namespace std; const int mod=1e9+7; long long fac[2000000]; long long inv[2000000]; long long invf[2000000]; long long modpow(long long a,long long n){ long long ret=1,t=a; for(int i=0;i<50;i++){ if(n>>i&1)ret=ret*t%mod; t=t*t%mod; } return ret; } long long modcomb(int n,int r){ return fac[n]*invf[r]%mod*invf[n-r]%mod; } int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); long long H,W; cin>>H>>W; fac[0]=inv[0]=invf[0]=1; fac[1]=inv[1]=invf[1]=1; for(int i=2;i<2e6;i++){ fac[i]=fac[i-1]*i%mod; inv[i]=inv[mod%i]*(mod-mod/i)%mod; invf[i]=invf[i-1]*inv[i]%mod; } long long ans=0; for(int i=0;i<=H;i++){ ans=(ans+modpow(modpow(2,H-i)-1,W)*(i%2==0?1:-1)*modcomb(H,i))%mod; } if(ans<0)ans+=mod; cout<<ans<<'\n'; }