結果

問題 No.1022 Power Equation
コンテスト
ユーザー batsumaru
提出日時 2020-04-10 23:29:23
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,158 bytes
コンパイル時間 1,713 ms
コンパイル使用メモリ 173,592 KB
実行使用メモリ 9,472 KB
最終ジャッジ日時 2024-09-16 02:23:01
合計ジャッジ時間 9,769 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample WA * 1
other WA * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<ll, ll> P;
#define DUMP(x) cout << #x << " = " << (x) << endl;
#define FOR(i, m, n) for (ll i = m; i < n; i++)
#define IFOR(i, m, n) for (ll i = n - 1; i >= m; i--)
#define REP(i, n) FOR(i, 0, n)
#define IREP(i, n) IFOR(i, 0, n)
#define FOREACH(x, a) for (auto&(x) : (a))
#define ALL(v) (v).begin(), (v).end()
#define SZ(x) ll(x.size())

ll gcd(ll a, ll b) {
  if (b == 0) return a;
  return gcd(b, a % b);
}

ll lcm(ll a, ll b) { return a / gcd(a, b) * b; }

void solve() {
  ll n;
  cin >> n;
  map<ll, bool> seen;
  ll bd = floor(sqrt(n));
  ll rem = max(n - bd, 0LL);
  ll ans = n * n;  // a = c = 1

  for (ll i = 2; i <= bd; i++) {
    if (seen[i]) {
      continue;
    }
    seen[i] = true;
    ll t = i;
    ll cnt = 0;
    while (t <= n) {
      t *= i;
      seen[t] = true;
      if (t > bd) {
        rem = max(rem - 1, 0LL);
      }
      cnt++;
    }
    FOR(j, 1, cnt + 1) FOR(k, 1, cnt + 1) {
      ll u = max(j, k) / gcd(j, k);
      ans += n / u;
    }
  }
  ans += rem * n;
  cout << ans << endl;
}

int main() {
  ll t;
  cin >> t;
  while (t--) solve();
}
0