結果

問題 No.1846 Good Binary Matrix
ユーザー boatmusclesboatmuscles
提出日時 2022-02-19 00:10:17
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 353 ms / 2,000 ms
コード長 890 bytes
コンパイル時間 311 ms
コンパイル使用メモリ 43,520 KB
実行使用メモリ 11,008 KB
最終ジャッジ日時 2024-06-29 10:03:15
合計ジャッジ時間 5,148 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 13 ms
11,008 KB
testcase_01 AC 14 ms
11,004 KB
testcase_02 AC 13 ms
10,904 KB
testcase_03 AC 12 ms
10,880 KB
testcase_04 AC 13 ms
10,880 KB
testcase_05 AC 12 ms
11,008 KB
testcase_06 AC 12 ms
10,880 KB
testcase_07 AC 12 ms
10,880 KB
testcase_08 AC 12 ms
10,888 KB
testcase_09 AC 12 ms
10,880 KB
testcase_10 AC 12 ms
10,880 KB
testcase_11 AC 13 ms
11,008 KB
testcase_12 AC 13 ms
10,960 KB
testcase_13 AC 13 ms
10,848 KB
testcase_14 AC 12 ms
11,008 KB
testcase_15 AC 14 ms
10,880 KB
testcase_16 AC 343 ms
10,880 KB
testcase_17 AC 352 ms
10,860 KB
testcase_18 AC 339 ms
10,964 KB
testcase_19 AC 341 ms
10,992 KB
testcase_20 AC 353 ms
10,900 KB
testcase_21 AC 92 ms
10,880 KB
testcase_22 AC 100 ms
10,900 KB
testcase_23 AC 13 ms
11,008 KB
testcase_24 AC 13 ms
10,924 KB
testcase_25 AC 256 ms
10,880 KB
testcase_26 AC 168 ms
10,860 KB
testcase_27 AC 20 ms
10,880 KB
testcase_28 AC 45 ms
10,836 KB
testcase_29 AC 116 ms
10,880 KB
testcase_30 AC 280 ms
10,880 KB
testcase_31 AC 248 ms
10,880 KB
testcase_32 AC 303 ms
11,008 KB
testcase_33 AC 319 ms
11,008 KB
testcase_34 AC 13 ms
10,860 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