結果
問題 | No.1383 Numbers of Product |
ユーザー | chocorusk |
提出日時 | 2021-02-07 21:32:47 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 703 ms / 2,000 ms |
コード長 | 1,731 bytes |
コンパイル時間 | 1,247 ms |
コンパイル使用メモリ | 131,468 KB |
最終ジャッジ日時 | 2025-01-18 14:34:11 |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 51 |
ソースコード
#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); if(p>n){ p=n+1; break; } } if(p==n+1){ myon=a; break; } mp[p]^=(1<<i); } 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; }