結果

問題 No.1706 Many Bus Stops (hard)
コンテスト
ユーザー kakel-san
提出日時 2023-12-30 21:57:37
言語 C#
(.NET 8.0.404)
結果
AC  
実行時間 59 ms / 2,000 ms
コード長 2,508 bytes
コンパイル時間 8,232 ms
コンパイル使用メモリ 166,404 KB
実行使用メモリ 186,364 KB
最終ジャッジ日時 2024-09-27 16:53:14
合計ジャッジ時間 11,343 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 41
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (102 ms)。
MSBuild のバージョン 17.9.6+a4ecab324 (.NET)
  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();
    }
    static void Solve()
    {
        var p = NList;
        var (c, n, m) = (p[0], p[1], p[2]);
        var r = Exp(c, mod - 2, mod);
        var r1 = (c - 1) * r % mod;
        var r2 = (c - 2) * r % mod;
        var x = new long[][]
        {
            new long[] { 0, 1, 0, 0 },
            new long[] { 0, r, r1, 0 },
            new long[] { 0, 0, 0, 1 },
            new long[] { r, 0, r2, r },
        };
        var a = new long[][] { new long[] { 0, 1, 0, 0 } };
        var b = Matrix.Mul(a, Matrix.Pow(x, n));
        var revb = (mod + 1 - b[0][1]) % mod;
        var ex = Exp(revb, m, mod);
        WriteLine((mod + 1 - ex) % mod);
    }
    static int mod = 1_000_000_007;
    static long Exp(long n, long p, int mod)
    {
        long _n = n % mod;
        var _p = p;
        var result = 1L;
        if ((_p & 1) == 1) result *= _n;
        while (_p > 0)
        {
            _n = _n * _n % mod;
            _p >>= 1;
            if ((_p & 1) == 1) result = result * _n % mod;
        }
        return result;
    }
    class Matrix
    {
        // 行列の累乗
        public static long[][] Pow(long[][] m, long k)
        {
            var multi = m;
            var r = new long[m.Length][];
            for (var i = 0; i < m.Length; ++i)
            {
                r[i] = new long[m.Length];
                r[i][i] = 1;
            }
            while (k > 0)
            {
                if ((k & 1) == 1) r = Mul(r, multi);
                multi = Mul(multi, multi);
                k >>= 1;
            }
            return r;
        }
        // 行列の積
        public static long[][] Mul(long[][] x, long[][] y)
        {
            var r = new long[x.Length][];
            for (var i = 0; i < x.Length; ++i) r[i] = new long[y[0].Length];
            for (var i = 0; i < x.Length; ++i) for (var k = 0; k < x[0].Length; ++k)
            {
                for (var j = 0; j < y[0].Length; ++j) r[i][j] = (r[i][j] + x[i][k] * y[k][j]) % mod;
            }
            return r;
        }
    }
}
0