結果
問題 |
No.3156 Count That Day's N
|
ユーザー |
![]() |
提出日時 | 2025-07-25 18:34:59 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 993 ms / 3,000 ms |
コード長 | 950 bytes |
コンパイル時間 | 3,112 ms |
コンパイル使用メモリ | 280,240 KB |
実行使用メモリ | 82,560 KB |
最終ジャッジ日時 | 2025-07-25 18:35:29 |
合計ジャッジ時間 | 25,571 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 32 |
ソースコード
#include <bits/stdc++.h> //#include <atcoder/modint> using namespace std; //using namespace atcoder; using ll = long long; //using mint = modint998244353; ll isqrt(ll a){ if (a == 0) return 0; ll l=1, r=3037000500, c; while(r-l>1){ c = (l+r)/2; if (c*c <= a) l=c; else r=c; } return l; } int main(){ cin.tie(nullptr); ios_base::sync_with_stdio(false); /* x^6+y^4 とかけるNは高々1e15^(1/4)*1e15^(1/6) = 1.7e6程度 */ ll N, K; cin >> K >> N; const ll a = 316, b = 5623; set<ll> st; for (ll x=1; x<=a; x++){ for (ll y=1; y<=b; y++){ ll n = x*x*x*x*x*x+y*y*y*y; if (n <= N) st.insert(n); } } ll ans=0; for (auto n : st){ if (n % K == 0){ ll z2 = n / K; ll z = isqrt(z2); if (z*z == z2) ans++; } } cout << ans << endl; return 0; }