結果
問題 | No.864 四方演算 |
ユーザー |
|
提出日時 | 2023-12-17 18:08:03 |
言語 | C++17(clang) (17.0.6 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 13 ms / 1,000 ms |
コード長 | 1,437 bytes |
コンパイル時間 | 1,128 ms |
コンパイル使用メモリ | 132,144 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-09-27 08:01:43 |
合計ジャッジ時間 | 2,167 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 27 |
ソースコード
#include <iostream>#include <vector>#include <algorithm>using namespace std;//約数列挙のアルゴリズムvector<long long> enum_divisors(long long N) {vector<long long> res;for (long long i = 1; i * i <= N; ++i) {if (N % i == 0) {res.push_back(i);// 重複しないならば i の相方である N/i も pushif (N / i != i) res.push_back(N / i);}}// 小さい順に並び替えるsort(res.begin(), res.end());return res;}int main(void){long long N,K,ans=0;cin >> N >> K;//step1: Kの約数列挙vector<long long> divisors = enum_divisors(K);long long a[2],sum=0;//step2: 個数のカウントfor (int i = 0; i < divisors.size() / 2; i++){sum = 0;a[0] = divisors[i];a[1] = divisors[divisors.size() - 1 - i];if (N < (a[0] - 1)) sum = max((long long)0, (a[0] - 1 - 2 * (a[0] - 1 - N)));else sum = a[0] - 1;if (N < (a[1] - 1)) sum *= max((long long)0, (a[1] - 1 - 2 * (a[1] - 1 - N)));else sum *= a[1] - 1;ans += (2*sum);}if (divisors.size() % 2){a[0] = divisors[divisors.size() / 2];long long tmp = max((long long)0, (a[0] - 1 - 2 * (a[0] - 1 - N)));if (N < (a[0] - 1)) sum = tmp * tmp;else sum = (a[0] - 1) * (a[0] - 1);ans += sum;}cout << ans << endl;}