結果

問題 No.864 四方演算
コンテスト
ユーザー mitas11110
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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 も push
            if (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;
}
0