結果
問題 | No.1844 Divisors Sum Sum |
ユーザー |
👑 |
提出日時 | 2022-02-18 21:59:13 |
言語 | C (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 544 ms / 3,000 ms |
コード長 | 833 bytes |
コンパイル時間 | 941 ms |
コンパイル使用メモリ | 31,872 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-06-29 08:51:35 |
合計ジャッジ時間 | 9,334 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 38 |
ソースコード
#include <stdio.h> const int Mod = 1000000007; long long div_mod(long long x, long long y, long long z) { if (x % y == 0) return x / y; else return (div_mod((1 + x / y) * y - x, (z % y), y) * z + x) / y; } long long pow_mod(int n, long long k) { long long N, ans = 1; for (N = n; k > 0; k >>= 1, N = N * N % Mod) if (k & 1) ans = ans * N % Mod; return ans; } int main() { int i, L, P[200001], e[200001]; scanf("%d", &L); for (i = 1; i <= L; i++) scanf("%d %d", &(P[i]), &(e[i])); long long ans = 1; for (i = 1; i <= L; i++) ans = ans * (((e[i] + 1) * div_mod(pow_mod(P[i], e[i] + 1) - 1, P[i] - 1, Mod) + div_mod(pow_mod(P[i], e[i]) - 1, pow_mod((P[i] - 1), 2), Mod) * P[i] % Mod + Mod - div_mod(e[i] * pow_mod(P[i], e[i] + 1) % Mod, P[i] - 1, Mod)) % Mod) % Mod; printf("%lld\n", ans); fflush(stdout); return 0; }