結果

問題 No.1846 Good Binary Matrix
コンテスト
ユーザー boatmuscles
提出日時 2022-02-19 00:10:17
言語 C++17
(gcc 15.2.0 + boost 1.90.0)
コンパイル:
g++-15 -O2 -lm -std=c++17 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 288 ms / 2,000 ms
コード長 890 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 344 ms
コンパイル使用メモリ 52,348 KB
実行使用メモリ 11,776 KB
最終ジャッジ日時 2026-06-25 08:11:35
合計ジャッジ時間 4,964 ms
ジャッジサーバーID
（参考情報） 		judge1_1 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 35
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#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