AOJ本の読書メモ   AOJ    次のAOJの問題へ    前のAOJの問題へ

NTL_1_D: Euler's Phi Function


問題へのリンク


C#のソース

using System;
using System.Collections.Generic;
using System.Linq;

// Q087 オイラーのトーシェント関数 https://judge.u-aizu.ac.jp/onlinejudge/description.jsp?id=NTL_1_D&lang=ja
class Program
{
    static string InputPattern = "InputX";

    static List<string> GetInputList()
    {
        var WillReturn = new List<string>();

        if (InputPattern == "Input1") {
            WillReturn.Add("6");
            //2
        }
        else if (InputPattern == "Input2") {
            WillReturn.Add("1000000");
            //400000
        }
        else {
            string wkStr;
            while ((wkStr = Console.ReadLine()) != null) WillReturn.Add(wkStr);
        }
        return WillReturn;
    }

    static void Main()
    {
        List<string> InputList = GetInputList();

        long N = long.Parse(InputList[0]);
        Dictionary<long, long> SoinsuuDict = DeriveSoinsuuDict(N);

        long Answer = N;
        foreach (var EachPair in SoinsuuDict) {
            Answer /= EachPair.Key;
        }
        foreach (var EachPair in SoinsuuDict) {
            Answer *= EachPair.Key - 1;
        }

        Console.WriteLine(Answer);
    }

    // 素因数分解し、指数[素数]なDictを返す
    static Dictionary<long, long> DeriveSoinsuuDict(long pTarget)
    {
        var WillReturn = new Dictionary<long, long>();

        long CurrVal = pTarget;
        // ルートより大きい数を、素因素に持つとしても、1つだけのため
        // ルートまで試し割りを行い、素因数が残っていたら、追加する
        for (long I = 2; I * I <= pTarget; I++) {
            if (CurrVal % I > 0) continue;

            WillReturn[I] = 0;
            while (CurrVal % I == 0) {
                WillReturn[I]++;
                CurrVal /= I;
            }
            if (CurrVal == 1) break;
        }

        // 素因数が残っている場合
        if (CurrVal > 1) {
            WillReturn[CurrVal] = 1;
        }
        return WillReturn;
    }
}


解説

オイラーのトーシェント関数を使ってます。