結果

問題 No.1383 Numbers of Product
ユーザー chocorusk
提出日時 2021-02-07 21:25:40
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,672 bytes
コンパイル時間 1,551 ms
コンパイル使用メモリ 131,752 KB
最終ジャッジ日時 2025-01-18 14:26:22
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 41 WA * 10
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
#include <fstream>
#include <utility>
#include <functional>
#include <stack>
#include <array>
#include <list>
#define popcount __builtin_popcount
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;

int main()
{
    ll n, k, M; cin>>n>>k>>M;
    using lll=__int128_t;
    map<ll, int> mp;
    for(int i=2; i<=19; i++){
        int myon=0;
        for(lll a=1; ; a++){
            lll p=1;
            for(lll j=0; j<=i; j++){
                p*=(a+j*k);
                p=min(p, (lll)(n+1));
            }
            mp[p]^=(1<<i);
            if(p==n+1){
                myon=a;
                break;
            }
        }
        if(myon==1) break;
    }
    int c=0;
    for(auto &p:mp){
        ll x=p.first;
        ll l=0, r=1e9;
        while(r-l>1){
            lll m=(l+r)/2;
            if(m*(m+k)>=x) r=m;
            else l=m;
        }
        if((lll)r*(r+k)==x){
            p.second^=2;
            c++;
        }
    }
    if(M>=2){
        int ans=0;
        for(auto p:mp){
            if(popcount(p.second)==M) ans++;
        }
        cout<<ans<<endl;
        return 0;
    }
    ll l=0, r=1e9;
    while(r-l>1){
        lll m=(l+r)/2;
        if(m*(m+k)<=n) l=m;
        else r=m;
    }
    ll ans=l-c;
    for(auto p:mp){
        if(popcount(p.second)==1) ans++;
    }
    cout<<ans<<endl;
    return 0;
}
0