結果
問題 |
No.1966 Median of Divisors
|
ユーザー |
![]() |
提出日時 | 2023-03-07 13:16:22 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 620 ms / 2,000 ms |
コード長 | 773 bytes |
コンパイル時間 | 4,292 ms |
コンパイル使用メモリ | 251,440 KB |
最終ジャッジ日時 | 2025-02-11 06:11:14 |
ジャッジサーバーID (参考情報) |
judge1 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 12 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/all> using namespace std; using namespace atcoder; using ll=long long; using ld=long double; ld pie=3.14159265359; ll inf=1001; ll mod=1000000007; ll modpow(ll x, ll n) { x = x%mod; if(n==0) return 1; //再帰の終了条件 else if(n%2==1) { return (x*modpow(x, n-1))%mod; //nが奇数ならnを1ずらす } else return modpow((x*x)%mod, n/2)%mod; //nが偶数ならnが半分になる } int main(){ ll t; cin >> t; for (ll i = 0; i < t; i++) { ll n,m; cin >> n >> m; ll x=modpow(n,m); ll ans=(x*(x+1))%mod; ans*=inv_mod(2,mod); ans%=mod; ll z=modpow(n,m/2); ll y=(z*(z+1))%mod; y*=(2*z+1)%mod; y%=mod; y*=inv_mod(6,mod); y%=mod; ans+=mod-y; ans%=mod; cout << ans << endl; } }