結果

問題 No.3004 ヤング図形
ユーザー kakel-san
提出日時 2025-01-17 22:29:35
言語 C#
(.NET 8.0.404)
結果
MLE  
実行時間 -
コード長 4,489 bytes
コンパイル時間 15,871 ms
コンパイル使用メモリ 170,588 KB
実行使用メモリ 853,116 KB
最終ジャッジ日時 2025-01-17 22:31:24
合計ジャッジ時間 98,510 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 8 TLE * 15 MLE * 2
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (287 ミリ秒)。
  main -> /home/judge/data/code/bin/Release/net8.0/main.dll
  main -> /home/judge/data/code/bin/Release/net8.0/publish/

ソースコード

diff # 

using System;
using static System.Console;
using System.Linq;
using System.Collections.Generic;

class Program
{
    static int NN => int.Parse(ReadLine());
    static int[] NList => ReadLine().Split().Select(int.Parse).ToArray();
    static int[][] NArr(long n) => Enumerable.Repeat(0, (int)n).Select(_ => NList).ToArray();
    public static void Main()
    {
        Solve();
        // Test();
    }
    static void Test()
    {
        var r = new Random();
        var count = 0;
        while (true)
        {
            var k = 1000;
            var map = Make(r, k);
            // 愚直解
            // var s1 = All(n, t, map);
            // 作成解
                // WriteLine(string.Join("\n", map.Select(m => string.Join(" ", m))));
            var s2 = Y(k, map);
            // if (s1 != s2)
            {
                // WriteLine(n + " " + t);
                // WriteLine(string.Join(" ", a));
                // WriteLine(s1);
                // WriteLine(s2);
                // return;
            }
            ++count;
            if (count % 1 == 0) WriteLine(count);
        }
    }
    static int[][] Make(Random r, int k)
    {
        var ans = new int[k][];
        var sum = 0L;
        for (var i = 0; i < k; ++i)
        {
            ans[i] = new int[] { i + 1, r.Next(1, 1000) };
            sum += ans[i][0] * ans[i][1];
        }
        while (sum > 100_000_000)
        {
            var pos = r.Next(k);
            if (ans[pos][1] > 1)
            {
                sum -= ans[pos][0] * (ans[pos][1] - 1);
                ans[pos][1] = 1;
            }
        }
        return ans;
    }
    static void Solve()
    {
        var k = NN;
        var map = NArr(k);
        WriteLine(Y(k, map));
    }
    static long Y(int k, int[][] map)
    {
        var sum = 0;
        foreach (var y in map)
        {
            sum += y[0] * y[1];
        }
        var mod = 998_244_353;
        var ncr = new NCR(sum / 2 + 1, mod);
        var ans = 1L;
        var rest = sum;
        foreach (var y in map)
        {
            var len = y[0];
            var count = y[1];
            ans = ans * ncr.Calc2(rest, len * count) % mod;
            var r2 = len * count;
            for (var i = 0; i < count; ++i)
            {
                ans = ans * ncr.Calc2(r2 - 1, len - 1) % mod;
                r2 -= len;
            }
            rest -= len * count;
        }
        return ans;
    }
    class NCR
    {
        int[] facts;
        int[] revFacts;
        int mod;
        public NCR(int n, int mod)
        {
            facts = new int[n + 1];
            revFacts = new int[n + 1];
            this.mod = mod;
            facts[0] = 1;
            var tmp = 1L;
            for (var i = 1; i <= n; ++i)
            {
                tmp = tmp * i % mod;
                facts[i] = (int)tmp;
            }
            tmp = Exp(facts[n], mod - 2);
            revFacts[n] = (int)tmp;
            for (var i = n; i > 1; --i)
            {
                tmp = tmp * i % mod;
                revFacts[i - 1] = (int)tmp;
            }
            revFacts[0] = 1;
        }
        public long Exp(long n, long k)
        {
            n = n % mod;
            if (k == 0) return 1;
            if (k == 1) return n;
            var half = Exp(n, k / 2);
            var result = half * half % mod;
            return ((k % 2) == 0) ? result : (result * n % mod);
        }
        public long Calc(int n, int r)
        {
            return (long)facts[n] * revFacts[r] % mod * revFacts[n - r] % mod;
        }
        /// <summary>nが大きくrが小さい場合の計算</summary>
        public long Calc2(int n, int r)
        {
            if (r > n - r) r = n - r;
            var tmp = 1L;
            for (var i = 0; i < r; ++i)
            {
                tmp = tmp * (n - i) % mod;
            }
            return tmp * revFacts[r] % mod;
        }
        public long NPR(int n, int r)
        {
            return (long)facts[n] * revFacts[r] % mod;
        }
        public long Fact(int n)
        {
            return facts[n];
        }
        public long RevFact(int n)
        {
            return revFacts[n];
        }
        public long Inverse(int n)
        {
            return (long)revFacts[n] * facts[n - 1] % mod;
        }
    }
}
0