結果
問題 | No.1383 Numbers of Product |
ユーザー | SSRS |
提出日時 | 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; } }