結果

問題 No.377 背景パターン
ユーザー EmKjpEmKjp
提出日時 2015-08-11 00:17:48
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 2,513 ms / 5,000 ms
コード長 4,808 bytes
コンパイル時間 1,694 ms
コンパイル使用メモリ 115,688 KB
実行使用メモリ 28,188 KB
最終ジャッジ日時 2024-07-18 06:36:15
合計ジャッジ時間 6,756 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 34 ms
26,404 KB
testcase_01 AC 30 ms
24,016 KB
testcase_02 AC 26 ms
24,312 KB
testcase_03 AC 27 ms
23,984 KB
testcase_04 AC 25 ms
26,140 KB
testcase_05 AC 23 ms
24,244 KB
testcase_06 AC 23 ms
21,808 KB
testcase_07 AC 24 ms
24,024 KB
testcase_08 AC 24 ms
26,068 KB
testcase_09 AC 25 ms
23,980 KB
testcase_10 AC 25 ms
24,284 KB
testcase_11 AC 26 ms
26,064 KB
testcase_12 AC 26 ms
24,028 KB
testcase_13 AC 39 ms
28,188 KB
testcase_14 AC 42 ms
24,352 KB
testcase_15 AC 28 ms
24,156 KB
testcase_16 AC 2,414 ms
26,324 KB
testcase_17 AC 2,513 ms
24,124 KB
testcase_18 AC 27 ms
24,148 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc)
Copyright (C) Microsoft Corporation. All rights reserved.

ソースコード

diff # 

using System;
using System.IO;
using System.Collections.Generic;
using System.Globalization;
using System.Linq;
using System.Text;
using System.Numerics;

partial class Solver {
    static public long Gcd(long s, long t) {
        return t != 0 ? Gcd(t, s % t) : s;
    }

    static public long Lcm(long s, long t) {
        return s / Gcd(s, t) * t;
    }

    static public List<int> Divisors(int x) {
        List<int> ret = new List<int>();
        for (int i = 1; i * i <= x; i++) {
            if (x % i == 0) {
                ret.Add(i);
                if (i != x / i) ret.Add(x / i);
            }
        }
        ret.Sort();
        return ret;
    }
    static public long Phi(long n) {
        var nn = n;
        long phi = n;

        for (long i = 2; i * i <= n && i < n; i++) {
            if (n % i == 0) {
                while (n % i == 0) n /= i;
                phi = phi / i * (i - 1);
            }
        }
        if (n > 1) phi = phi / n * (n - 1);

        return phi;
    }

    static public long ExtendedGcd(long a, long b, ref long x, ref long y) {
        if (b == 0) {
            x = 1; y = 0;
            return a;
        } else {
            long d = ExtendedGcd(b, a % b, ref y, ref x);
            y -= a / b * x;
            return d;
        }
    }

    static public long Inverse(long a, long mod) {
        long x = 0, y = 0;
        if (ExtendedGcd(a, mod, ref x, ref y) == 1)
            return (x + mod) % mod;
        else
            return -1;
    }

    static public long ModPow(long x, long n, long mod) {
        long t = x, ret = 1;
        while (n > 0) {
            if ((n & 1) != 0) {
                ret = (t * ret) % mod;
            }
            t = t * t % mod;
            n >>= 1;
        }
        return ret;
    }

    public void Run() {
        var H = ni();
        var W = ni();
        var K = ni();
        const int mod = 1000000007;
        var dH = Divisors(H);
        var dW = Divisors(W);
        long ans = 0;
        foreach (var h in dH) {
            foreach (var w in dW) {
                ans += (long)(Phi(h) * Phi(w) % mod * ModPow(K, 1L * W * H / w / h * Gcd(w, h), mod) % mod);
                ans %= mod;
            }
        }
        ans = (ans * Inverse(1L * W * H, mod)) % mod;

        cout.WriteLine(ans);
    }
}

// PREWRITEN CODE BEGINS FROM HERE
partial class Solver : Scanner {
    public static void Main(string[] args) {
        new Solver(Console.In, Console.Out).Run();
    }

    TextReader cin;
    TextWriter cout;

    public Solver(TextReader reader, TextWriter writer)
        : base(reader) {
        this.cin = reader;
        this.cout = writer;
    }
    public Solver(string input, TextWriter writer)
        : this(new StringReader(input), writer) {
    }

    public int ni() { return NextInt(); }
    public int[] ni(int n) { return NextIntArray(n); }
    public long nl() { return NextLong(); }
    public long[] nl(int n) { return NextLongArray(n); }
    public double nd() { return NextDouble(); }
    public string ns() { return Next(); }
}

public class Scanner {
    private TextReader Reader;
    private Queue<String> TokenQueue = new Queue<string>();
    private CultureInfo ci = CultureInfo.InvariantCulture;

    public Scanner()
        : this(Console.In) {
    }

    public Scanner(TextReader reader) {
        this.Reader = reader;
    }

    public int NextInt() { return Int32.Parse(Next(), ci); }
    public long NextLong() { return Int64.Parse(Next(), ci); }
    public double NextDouble() { return double.Parse(Next(), ci); }
    public string[] NextArray(int size) {
        var array = new string[size];
        for (int i = 0; i < size; i++) array[i] = Next();
        return array;
    }
    public int[] NextIntArray(int size) {
        var array = new int[size];
        for (int i = 0; i < size; i++) array[i] = NextInt();
        return array;
    }

    public long[] NextLongArray(int size) {
        var array = new long[size];
        for (int i = 0; i < size; i++) array[i] = NextLong();
        return array;
    }

    public String Next() {
        if (TokenQueue.Count == 0) {
            if (!StockTokens()) throw new InvalidOperationException();
        }
        return TokenQueue.Dequeue();
    }

    public bool HasNext() {
        if (TokenQueue.Count > 0)
            return true;
        return StockTokens();
    }

    private bool StockTokens() {
        while (true) {
            var line = Reader.ReadLine();
            if (line == null) return false;
            var tokens = line.Trim().Split(" ".ToCharArray(), StringSplitOptions.RemoveEmptyEntries);
            if (tokens.Length == 0) continue;
            foreach (var token in tokens)
                TokenQueue.Enqueue(token);
            return true;
        }
    }
}
0