結果
問題 | 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;}}