結果

問題 No.886 Direct
ユーザー 0w10w1
提出日時 2019-09-28 04:51:00
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 450 ms / 4,000 ms
コード長 1,288 bytes
コンパイル時間 2,588 ms
コンパイル使用メモリ 206,284 KB
実行使用メモリ 97,116 KB
最終ジャッジ日時 2024-09-25 10:36:46
合計ジャッジ時間 7,496 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 2 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 2 ms
5,376 KB
testcase_13 AC 2 ms
5,376 KB
testcase_14 AC 2 ms
5,376 KB
testcase_15 AC 2 ms
5,376 KB
testcase_16 AC 2 ms
5,376 KB
testcase_17 AC 3 ms
5,376 KB
testcase_18 AC 3 ms
5,376 KB
testcase_19 AC 3 ms
5,376 KB
testcase_20 AC 2 ms
5,376 KB
testcase_21 AC 4 ms
5,376 KB
testcase_22 AC 3 ms
5,376 KB
testcase_23 AC 133 ms
47,608 KB
testcase_24 AC 156 ms
51,752 KB
testcase_25 AC 84 ms
33,120 KB
testcase_26 AC 113 ms
42,396 KB
testcase_27 AC 356 ms
90,296 KB
testcase_28 AC 328 ms
87,428 KB
testcase_29 AC 369 ms
97,040 KB
testcase_30 AC 394 ms
96,984 KB
testcase_31 AC 409 ms
97,116 KB
testcase_32 AC 450 ms
96,988 KB
testcase_33 AC 389 ms
97,108 KB
testcase_34 AC 449 ms
96,984 KB
testcase_35 AC 413 ms
97,104 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

const int MOD = 1e9 + 7;

int main() {
  int H, W;
  cin >> H >> W;

  vector<int> f(H);
  for (int i = 1; i < H; ++i) f[i] = H - i;

  vector<int> g(W);
  for (int i = 1; i < W; ++i) g[i] = W - i;

  int n = max(H, W) + 1;
  f.resize(n);
  g.resize(n);

  auto convolve = [](vector<int> a) {
    vector<int> b = a;
    vector<int> np(a.size());
    for (int i = 2; i < a.size(); ++i) {
      if (np[i]) continue;
      for (int j = (a.size() - 1) / i; j; --j) {
        np[j * i] = 1;
        (b[j] += b[j * i]) %= MOD;
      }
    }
    return b;
  };

  auto inv_convolve = [](vector<int> a) {
    vector<int> b = a;
    vector<int> np(a.size());
    for (int i = 2; i < a.size(); ++i) {
      if (np[i]) continue;
      for (int j = 1; j * i < a.size(); ++j) {
        np[j * i] = 1;
        (b[j] -= b[j * i]) %= MOD;
      }
    }
    return b;
  };

  vector<int> cf = convolve(f);
  vector<int> cg = convolve(g);

  vector<int> cfg(n);
  for (int i = 1; i < n; ++i) {
    cfg[i] = 1LL * cf[i] * cg[i] % MOD;
  }

  vector<int> fg = inv_convolve(cfg);

  int ans = 2 * fg[1] % MOD;
  (ans += 1LL * W * (H - 1) % MOD) %= MOD;
  (ans += 1LL * H * (W - 1) % MOD) %= MOD;

  if (ans < 0) ans += MOD;
  cout << ans << endl;

  return 0;
}
0