結果
問題 | No.1706 Many Bus Stops (hard) |
ユーザー |
|
提出日時 | 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/
ソースコード
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;}}}