結果

問題 No.1846 Good Binary Matrix
ユーザー boatmusclesboatmuscles
提出日時 2022-02-19 00:10:17
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 391 ms / 2,000 ms
コード長 890 bytes
コンパイル時間 305 ms
コンパイル使用メモリ 43,124 KB
実行使用メモリ 10,868 KB
最終ジャッジ日時 2023-09-11 20:18:34
合計ジャッジ時間 6,063 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 14 ms
10,860 KB
testcase_01 AC 13 ms
10,860 KB
testcase_02 AC 15 ms
10,844 KB
testcase_03 AC 13 ms
10,860 KB
testcase_04 AC 13 ms
10,864 KB
testcase_05 AC 14 ms
10,732 KB
testcase_06 AC 14 ms
10,856 KB
testcase_07 AC 14 ms
10,852 KB
testcase_08 AC 13 ms
10,844 KB
testcase_09 AC 14 ms
10,852 KB
testcase_10 AC 14 ms
10,864 KB
testcase_11 AC 13 ms
10,832 KB
testcase_12 AC 13 ms
10,848 KB
testcase_13 AC 13 ms
10,804 KB
testcase_14 AC 13 ms
10,796 KB
testcase_15 AC 14 ms
10,732 KB
testcase_16 AC 372 ms
10,792 KB
testcase_17 AC 385 ms
10,848 KB
testcase_18 AC 368 ms
10,792 KB
testcase_19 AC 361 ms
10,848 KB
testcase_20 AC 391 ms
10,848 KB
testcase_21 AC 100 ms
10,736 KB
testcase_22 AC 108 ms
10,828 KB
testcase_23 AC 14 ms
10,800 KB
testcase_24 AC 13 ms
10,848 KB
testcase_25 AC 274 ms
10,844 KB
testcase_26 AC 181 ms
10,848 KB
testcase_27 AC 21 ms
10,860 KB
testcase_28 AC 49 ms
10,796 KB
testcase_29 AC 123 ms
10,828 KB
testcase_30 AC 299 ms
10,800 KB
testcase_31 AC 262 ms
10,796 KB
testcase_32 AC 318 ms
10,868 KB
testcase_33 AC 345 ms
10,860 KB
testcase_34 AC 14 ms
10,792 KB
権限があれば一括ダウンロードができます

ソースコード

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