結果

問題 No.3296 81-like number
ユーザー GrayCoder
提出日時 2025-10-05 14:53:28
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 3 ms / 2,000 ms
コード長 1,389 bytes
コンパイル時間 2,795 ms
コンパイル使用メモリ 281,416 KB
実行使用メモリ 7,716 KB
最終ジャッジ日時 2025-10-05 14:53:51
合計ジャッジ時間 3,576 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = int64_t;
#define all(a) begin(a), end(a)
#define rep(i, s, e) for (int i = s; i < e; ++i)
#define asort(a) stable_sort(all(a))
#define rsort(a) stable_sort(rbegin(a), rend(a))
#define sum(a) reduce(all(a))
#define min_(a) reduce(all(a), *begin(a), [](auto x, auto y) { return min(x, y); })
#define max_(a) reduce(all(a), *begin(a), [](auto x, auto y) { return max(x, y); })
constexpr double degrees(double a) { return (a * 180.0) / M_PI; }
constexpr double radians(double a) { return (a * M_PI) / 180.0; }

vector<ll> Eratosthenes(const ll n) {
  vector<bool> is_prime(n + 1);
  for (ll i = 0; i <= n; ++i) {
    is_prime[i] = true;
  }
  vector<ll> p;
  for (ll i = 2; i <= n; ++i) {
    if (is_prime[i]) {
      for (ll j = 2 * i; j <= n; j += i) {
        is_prime[j] = false;
      }
      p.emplace_back(i);
    }
  }
  return p;
}

int main(int, [[maybe_unused]] char *argv[]) {
#ifdef DEBUG
  freopen(argv[1], "r", stdin);
#endif

  ios_base::sync_with_stdio(false);
  cin.tie(0);
  cout << fixed << setprecision(15);

  ll N;
  cin >> N;
  vector<ll> p = Eratosthenes(100000);
  ll ans = 0;
  ll x;
  for (ll i = 2; i < N; ++i) {
    if (pow(2, i) > N) break;
    for (ll j : p) {
      x = pow(j, i);
      if (x <= N) {
        ans += x;
      } else {
        break;
      }
    }
  }
  cout << ans << endl;
  return 0;
}
0