結果
問題 | No.1383 Numbers of Product |
ユーザー |
![]() |
提出日時 | 2021-02-07 21:08:00 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 566 ms / 2,000 ms |
コード長 | 1,896 bytes |
コンパイル時間 | 1,746 ms |
コンパイル使用メモリ | 175,244 KB |
実行使用メモリ | 68,352 KB |
最終ジャッジ日時 | 2024-11-22 23:11:51 |
合計ジャッジ時間 | 16,370 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 51 |
ソースコード
#include <bits/stdc++.h>using namespace std;bool ok(long long K, long long X){long long tv = 1000000000;long long fv = 0;while (tv - fv > 1){long long mid = (tv + fv) / 2;double prod = 1;prod *= mid;prod *= mid + K;if (prod > X * 2){tv = mid;} else {long long prod2 = 1;prod2 *= mid;prod2 *= mid + K;if (prod2 >= X){tv = mid;} else {fv = mid;}}}return tv * (tv + K) == X;}int main(){long long N, K, M;cin >> N >> K >> M;map<long long, int> mp;int B = 2;while (true){long long A = 1;while (true){double prod = 1;for (int i = 0; i <= B; i++){prod *= A + K * i;}if (prod > N * 2){break;}long long prod2 = 1;for (int i = 0; i <= B; i++){prod2 *= A + K * i;}if (prod2 > N){break;}mp[prod2]++;A++;}if (A == 1){break;}B++;}if (M >= 2){int ans = 0;for (auto P : mp){long long x = P.first;int cnt = P.second;if (cnt == M - 1 && ok(K, x) || cnt == M && !ok(K, x)){ans++;}}cout << ans << endl;} else {long long tv = 0;long long fv = 1000000000;while (fv - tv > 1){long long mid = (tv + fv) / 2;double prod = 1;prod *= mid;prod *= mid + K;if (prod > N * 2){fv = mid;} else {long long prod2 = 1;prod2 *= mid;prod2 *= mid + K;if (prod2 <= N){tv = mid;} else {fv = mid;}}}long long ans = tv;for (auto P : mp){long long x = P.first;int cnt = P.second;if (ok(K, x)){ans--;} else if (cnt == 1){ans++;}}cout << ans << endl;}}