結果

問題 No.2249 GCDistance
コンテスト
ユーザー kakel-san
提出日時 2023-03-17 23:25:40
言語 C#(csc)
(csc 3.9.0)
結果
TLE  
実行時間 -
コード長 2,312 bytes
コンパイル時間 2,697 ms
コンパイル使用メモリ 112,864 KB
実行使用メモリ 186,868 KB
最終ジャッジ日時 2024-09-18 12:29:42
合計ジャッジ時間 15,681 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample TLE * 1
other TLE * 10
権限があれば一括ダウンロードができます
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc)
Copyright (C) Microsoft Corporation. All rights reserved.

ソースコード

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 p = PDiv(i, minPDivList).Distinct();
            var phi = i;
            foreach (var pi in p) phi = phi / pi * (pi - 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