結果
問題 | No.1966 Median of Divisors |
ユーザー |
![]() |
提出日時 | 2022-06-03 23:19:25 |
言語 | D (dmd 2.109.1) |
結果 |
AC
|
実行時間 | 228 ms / 2,000 ms |
コード長 | 775 bytes |
コンパイル時間 | 2,301 ms |
コンパイル使用メモリ | 204,708 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-06-22 15:23:32 |
合計ジャッジ時間 | 5,330 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 12 |
ソースコード
import std;void main() {long T;readf("%d\n", T);long MOD = 10 ^^ 9 + 7;long inv2 = powMod(2, MOD-2, MOD);long inv3 = powMod(3, MOD-2, MOD);foreach (i; 0 .. T) {long N, M;readf("%d %d\n", N, M);long a = powMod(N, M, MOD), b = powMod(N, M/2, MOD);long u = (a * (a + 1) % MOD) * inv2 % MOD;long v = (((((b * (b + 1) % MOD) * inv2) % MOD) * (b * 2 + 1)) % MOD) * inv3 % MOD;long res = (MOD + u - v) % MOD;res.writeln;}}long powMod(long x, long y, long z) {long res = 1;while (y > 0) {if (y % 2 == 1) {res *= x;if (res > z) res %= z;}x *= x;if (x > z) x %= z;y /= 2;}return res;}