結果
| 問題 |
No.1966 Median of Divisors
|
| コンテスト | |
| ユーザー |
umezo
|
| 提出日時 | 2022-06-04 23:32:35 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 190 ms / 2,000 ms |
| コード長 | 622 bytes |
| コンパイル時間 | 1,729 ms |
| コンパイル使用メモリ | 193,532 KB |
| 最終ジャッジ日時 | 2025-01-29 18:13:07 |
|
ジャッジサーバーID (参考情報) |
judge1 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 |
| other | AC * 12 |
ソースコード
#define rep(i,n) for(int i=0;i<(int)(n);i++)
#define ALL(v) v.begin(),v.end()
typedef long long ll;
#include<bits/stdc++.h>
using namespace std;
const int MOD=1e9+7;
ll modpow(ll x,ll n){
ll ans=1;
while(n){
if(n&1) ans=ans*x%MOD;
x=x*x%MOD;
n/=2;
}
return ans;
}
int main(){
ios::sync_with_stdio(false);
std::cin.tie(nullptr);
ll e=modpow(2,MOD-2),f=modpow(6,MOD-2);
int t;
cin>>t;
while(t--){
ll n,m;
cin>>n>>m;
ll c=modpow(n,m),d=modpow(n,m/2);
ll a=c*(c+1)%MOD*e%MOD;
ll b=d*(d+1)%MOD*(2*d%MOD+1)%MOD*f%MOD;
cout<<(a-b+MOD)%MOD<<endl;
}
return 0;
}
umezo