結果

問題 No.3156 Count That Day's N
コンテスト
ユーザー srjywrdnprkt
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
}
0