結果

問題 No.2249 GCDistance
ユーザー kakel-sankakel-san
提出日時 2023-03-17 23:33:35
言語 C#
(.NET 8.0.203)
結果
AC  
実行時間 2,156 ms / 5,000 ms
コード長 2,420 bytes
コンパイル時間 8,300 ms
コンパイル使用メモリ 166,364 KB
実行使用メモリ 389,364 KB
最終ジャッジ日時 2024-09-18 12:34:32
合計ジャッジ時間 36,071 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2,131 ms
210,796 KB
testcase_01 AC 2,098 ms
236,288 KB
testcase_02 AC 2,099 ms
236,544 KB
testcase_03 AC 2,097 ms
236,800 KB
testcase_04 AC 2,074 ms
219,904 KB
testcase_05 AC 2,101 ms
236,544 KB
testcase_06 AC 2,132 ms
236,672 KB
testcase_07 AC 2,100 ms
236,416 KB
testcase_08 AC 2,156 ms
211,200 KB
testcase_09 AC 2,108 ms
236,416 KB
testcase_10 AC 2,068 ms
389,364 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (92 ms)。
MSBuild のバージョン 17.9.6+a4ecab324 (.NET)
  main -> /home/judge/data/code/bin/Release/net8.0/main.dll
  main -> /home/judge/data/code/bin/Release/net8.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