結果

問題 No.93 ペガサス
ユーザー AreTrash
提出日時 2019-02-07 18:05:45
言語 C#(csc)
(csc 3.9.0)
結果
AC  
実行時間 622 ms / 5,000 ms
コード長 3,454 bytes
コンパイル時間 1,947 ms
コンパイル使用メモリ 113,784 KB
実行使用メモリ 36,096 KB
最終ジャッジ日時 2024-06-25 06:54:11
合計ジャッジ時間 6,117 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 16
権限があれば一括ダウンロードができます
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc)
Copyright (C) Microsoft Corporation. All rights reserved.

ソースコード

diff # 

using System;

namespace No93
{
    public class Program
    {
        static void Main()
        {
            var N = int.Parse(Console.ReadLine());
            var dp = new ModInt[N + 2, N + 2, 2, 2];
            dp[1, 0, 0, 0] = 1;
            dp[2, 0, 0, 0] = 2;

            for (var n = 3; n <= N; n++)
            for (var k = 0; k <= n - 2; k++)
            {
                dp[n, k + 1, 1, 0] += dp[n - 1, k, 0, 0] * 2;
                dp[n, k + 1, 1, 0] += dp[n - 1, k, 0, 1] * 1;
                dp[n, k + 1, 1, 1] += dp[n - 1, k, 1, 0] * 2;
                dp[n, k + 1, 1, 1] += dp[n - 1, k, 1, 1] * 1;

                dp[n, k, 0, 0] += dp[n - 1, k, 0, 0] * (n - k - 2);
                dp[n, k, 0, 0] += dp[n - 1, k, 0, 1] * (n - k - 1);
                dp[n, k, 0, 1] += dp[n - 1, k, 1, 0] * (n - k - 2);
                dp[n, k, 0, 1] += dp[n - 1, k, 1, 1] * (n - k - 1);
                dp[n, k, 1, 0] += dp[n - 1, k, 0, 1] * 1;
                dp[n, k, 1, 1] += dp[n - 1, k, 1, 1] * 1;

                if (k <= 0) continue;

                dp[n, k - 1, 0, 0] += dp[n - 1, k, 0, 0] * k;
                dp[n, k - 1, 0, 0] += dp[n - 1, k, 1, 0] * 1;
                dp[n, k - 1, 0, 0] += dp[n - 1, k, 0, 1] * (k - 1);
                dp[n, k - 1, 0, 0] += dp[n - 1, k, 1, 1] * 1;
                dp[n, k - 1, 0, 1] += dp[n - 1, k, 1, 0] * (k - 1);
                dp[n, k - 1, 0, 1] += dp[n - 1, k, 1, 1] * (k - 2);
            }

            Console.WriteLine(dp[N, 0, 0, 0]);
        }
    }

    public struct ModInt
    {
        public static int Mod = 1000000007;

        readonly long value;
        public int Value { get { return (int)value; } }
        public ModInt Inverse { get { return GetInverse(); } }

        public ModInt(long value)
        {
            value %= Mod;
            this.value = value < 0 ? value + Mod : value;
        }

        ModInt GetInverse()
        {
            if (value == 0) throw new InvalidOperationException("value must NOT equal 0");
            return Pow(Mod - 2);
        }

        public ModInt Square()
        {
            return this * this;
        }

        public ModInt Pow(long exp)
        {
            if (exp == 0) return 1;
            if (exp % 2 == 0) return Pow(exp / 2).Square();
            return this * Pow(exp - 1);
        }

        public static ModInt Parse(string str)
        {
            return long.Parse(str);
        }

        public static implicit operator ModInt(long value)
        {
            return new ModInt(value);
        }

        public static ModInt operator +(ModInt left, ModInt right)
        {
            return left.value + right.value;
        }

        public static ModInt operator -(ModInt left, ModInt right)
        {
            return left.value - right.value;
        }

        public static ModInt operator *(ModInt left, ModInt right)
        {
            return left.value * right.value;
        }

        public static ModInt operator /(ModInt left, ModInt right)
        {
            return left * right.Inverse;
        }

        public override bool Equals(object obj)
        {
            if (!(obj is ModInt)) return false;
            return value == ((ModInt)obj).value;
        }

        public override int GetHashCode()
        {
            return value.GetHashCode();
        }

        public override string ToString()
        {
            return value.ToString();
        }
    }
}
0