結果
問題 |
No.1035 Color Box
|
ユーザー |
![]() |
提出日時 | 2020-04-29 07:40:02 |
言語 | C#(csc) (csc 3.9.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 2,151 bytes |
コンパイル時間 | 2,508 ms |
コンパイル使用メモリ | 107,648 KB |
実行使用メモリ | 18,048 KB |
最終ジャッジ日時 | 2024-11-27 11:52:46 |
合計ジャッジ時間 | 4,345 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 9 WA * 27 |
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc) Copyright (C) Microsoft Corporation. All rights reserved.
ソースコード
using System; class Modulo { private const int MOD = 1000000007; private readonly int[] m_facs; public int Mul(int a, int b) => (int)(Math.BigMul(a, b) % MOD); public Modulo(int n) { m_facs = new int[n + 1]; m_facs[0] = 1; for (int i = 1; i <= n; ++i) m_facs[i] = Mul(m_facs[i - 1], i); } public int Fac(int n) => m_facs[n]; public int Pow(int a, int m) { switch (m) { case 0: return 1; case 1: return a; default: int p1 = Pow(a, m / 2); int p2 = Mul(p1, p1); return ((m % 2) == 0) ? p2 : Mul(p2, a); } } public int Div(int a, int b) => Mul(a, Pow(b, MOD - 2)); public int Ncr(int n, int r) { if (n < r) return 0; if (n == r) return 1; int res = Fac(n); res = Div(res, Fac(r)); res = Div(res, Fac(n - r)); return res; } } public class Hello { static int MOD = 1000000007; static void Main() { string[] line = Console.ReadLine().Trim().Split(' '); var n = int.Parse(line[0]); var m = int.Parse(line[1]); var md = new Modulo(n); var t1 = md.Ncr(n, m); var t2 = md.Fac(m); var t3 = md.Pow(m, n - m); long ans = (t1 * t2) % MOD; ans *= t3; ans %= MOD; if (m < n) { var t = getInv(2, MOD); ans *= t; ans %= MOD; } Console.WriteLine(ans); } public static long getInv(int a, int m) { var x = 0L; var y = 1L; var gcd = (long)m; var px = 1L; var py = 0L; var pgcd = (long)a; var temp = 0L; while (gcd > 0) { var quotient = pgcd / gcd; temp = x; x = px - quotient * temp; px = temp; temp = y; y = py - quotient * temp; py = temp; temp = gcd; gcd = pgcd - quotient * temp; pgcd = temp; } return px < 0 ? px + m : px; } }