結果

問題 No.886 Direct
ユーザー 0w10w1
提出日時 2019-09-28 04:51:00
言語 C++17
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 261 ms / 4,000 ms
コード長 1,288 bytes
コンパイル時間 2,203 ms
コンパイル使用メモリ 206,100 KB
実行使用メモリ 99,960 KB
最終ジャッジ日時 2023-10-26 02:35:08
合計ジャッジ時間 6,759 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,348 KB
testcase_01 AC 2 ms
4,348 KB
testcase_02 AC 2 ms
4,348 KB
testcase_03 AC 2 ms
4,348 KB
testcase_04 AC 1 ms
4,348 KB
testcase_05 AC 2 ms
4,348 KB
testcase_06 AC 2 ms
4,348 KB
testcase_07 AC 1 ms
4,348 KB
testcase_08 AC 2 ms
4,348 KB
testcase_09 AC 1 ms
4,348 KB
testcase_10 AC 2 ms
4,348 KB
testcase_11 AC 2 ms
4,348 KB
testcase_12 AC 1 ms
4,348 KB
testcase_13 AC 2 ms
4,348 KB
testcase_14 AC 2 ms
4,348 KB
testcase_15 AC 2 ms
4,348 KB
testcase_16 AC 2 ms
4,348 KB
testcase_17 AC 2 ms
4,348 KB
testcase_18 AC 2 ms
4,348 KB
testcase_19 AC 2 ms
4,348 KB
testcase_20 AC 2 ms
4,348 KB
testcase_21 AC 3 ms
4,400 KB
testcase_22 AC 3 ms
4,348 KB
testcase_23 AC 113 ms
48,732 KB
testcase_24 AC 125 ms
53,744 KB
testcase_25 AC 73 ms
34,264 KB
testcase_26 AC 99 ms
44,012 KB
testcase_27 AC 240 ms
93,596 KB
testcase_28 AC 226 ms
90,332 KB
testcase_29 AC 259 ms
98,740 KB
testcase_30 AC 258 ms
98,040 KB
testcase_31 AC 261 ms
99,960 KB
testcase_32 AC 260 ms
99,960 KB
testcase_33 AC 259 ms
99,960 KB
testcase_34 AC 259 ms
99,960 KB
testcase_35 AC 258 ms
99,920 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