結果

問題 No.2249 GCDistance
ユーザー 👑 kakel-sankakel-san
提出日時 2023-03-17 23:33:35
言語 C#
(.NET 8.0.203)
結果
AC  
実行時間 2,332 ms / 5,000 ms
コード長 2,420 bytes
コンパイル時間 7,049 ms
コンパイル使用メモリ 159,808 KB
実行使用メモリ 383,840 KB
最終ジャッジ日時 2023-10-18 16:33:19
合計ジャッジ時間 36,580 ms
ジャッジサーバーID
(参考情報)
judge13 / judge12
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2,219 ms
212,796 KB
testcase_01 AC 2,292 ms
239,044 KB
testcase_02 AC 2,298 ms
239,048 KB
testcase_03 AC 2,332 ms
239,044 KB
testcase_04 AC 2,255 ms
223,940 KB
testcase_05 AC 2,311 ms
239,052 KB
testcase_06 AC 2,312 ms
239,044 KB
testcase_07 AC 2,331 ms
239,048 KB
testcase_08 AC 2,235 ms
213,448 KB
testcase_09 AC 2,314 ms
239,044 KB
testcase_10 AC 2,315 ms
383,840 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (103 ms)。
MSBuild のバージョン 17.7.3+8ec440e68 (.NET)
  main -> /home/judge/data/code/bin/Release/net7.0/main.dll
  main -> /home/judge/data/code/bin/Release/net7.0/publish/

ソースコード

diff #

using System;
using static System.Console;
using System.Linq;
using System.Collections.Generic;
class Program
{
    static int NN => int.Parse(ReadLine());
    static int[] LList(int n) => Enumerable.Repeat(0, n).Select(_ => int.Parse(ReadLine())).ToArray();
    public static void Main()
    {
        Solve();
    }
    static void Solve()
    {
        var t = NN;
        var a = LList(t);
        WriteLine(string.Join("\n", GCDistance(a)));
    }
    static long[] GCDistance(int[] a)
    {
        var max = 10_000_001;
        var phis = Phis(max);
        var cum = new long[max];
        for (var i = 1; i < cum.Length; ++i)
        {
            cum[i] = cum[i - 1] + phis[i];
        }
        var ans = new long[a.Length];
        for (var i = 0; i < a.Length; ++i)
        {
            ans[i] = (long)a[i] * (a[i] - 1) - cum[a[i]];
        }
        return ans;
    }
    static int[] Phis(int max)
    {
        var ans = new int[max + 1];
        var minPDivList = GetMinPDivList(max);
        for (var i = 1; i < ans.Length; ++i)
        {
            var plist = PDiv(i, minPDivList);
            var phi = i;
            for (var j = 0; j < plist.Count; ++j)
            {
                if (j == 0 || plist[j] != plist[j - 1]) phi = phi / plist[j] * (plist[j] - 1);
            }
            ans[i] = phi;
        }
        return ans;
    }
    static int[] GetMinPDivList(int max)
    {
        var minPDivList = new int[max + 1];
        for (var i = 0; i <= max; ++i) minPDivList[i] = i;

        for (var i = 2; i * i <= max; ++i) if (minPDivList[i] == i)
        {
            for (var j = i * i; j <= max; j += i) if (minPDivList[j] == j)
            {
                minPDivList[j] = i;
            }
        }
        return minPDivList;
    }
    static int GetMinPDiv(int n, int min)
    {
        var p = min;
        while ((long)p * p <= n)
        {
            if (n % p == 0) return p;
            ++p;
        }
        return n;
    }
    static List<int> PDiv(int n, int[] minPDivList)
    {
        var p = minPDivList[n];
        if (n == p)
        {
            var dic = new List<int>();
            dic.Add(p);
            return dic;
        }
        else
        {
            var dic = PDiv(n / p, minPDivList);
            dic.Add(p);
            return dic;
        }
    }
}
0