結果

問題 No.1846 Good Binary Matrix
ユーザー boatmusclesboatmuscles
提出日時 2022-02-19 00:10:17
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 409 ms / 2,000 ms
コード長 890 bytes
コンパイル時間 389 ms
コンパイル使用メモリ 39,612 KB
最終ジャッジ日時 2025-01-28 00:48:51
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 35
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:9:10: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
    9 |     scanf("%d %d", &H, &W);
      |     ~~~~~^~~~~~~~~~~~~~~~~

ソースコード

diff # 

#include<cstdio>
#include<atcoder/modint>
using namespace atcoder;
using Mint = modint1000000007;
constexpr int MAX = 1000000;
using ll = long long;
int main(){
    int H, W;
    scanf("%d %d", &H, &W);
    Mint fact[MAX+1];
    fact[0] = Mint::raw(1);
    for (int i = 1; i <= MAX; i++) fact[i] = fact[i-1]*Mint::raw(i);
    Mint ifact[MAX+1];
    ifact[MAX] = fact[MAX].inv();
    for (int i = MAX; (i); i--) ifact[i-1] = ifact[i]*Mint::raw(i);
    Mint answer = ifact[H]*((Mint::raw(1) - ifact[2].pow(H)).pow(W) - Mint::raw(1));
    for (int i = 1; i < H; i++)
    {
        if(i&1) answer -= ifact[i]*ifact[H-i]*ifact[2].pow((ll)W*i)*(Mint::raw(1) - ifact[2].pow(H-i)).pow(W);
        else answer += ifact[i]*ifact[H-i]*ifact[2].pow((ll)W*i)*(Mint::raw(1) - ifact[2].pow(H-i)).pow(W);
    }
    printf("%u\n", ((answer*fact[H] + Mint::raw(1))*Mint::raw(2).pow((long long)H*W)).val());
}
0