結果

問題 No.1966 Median of Divisors
ユーザー asaringoasaringo
提出日時 2022-06-03 22:41:49
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 93 ms / 2,000 ms
コード長 1,064 bytes
コンパイル時間 1,981 ms
コンパイル使用メモリ 194,076 KB
最終ジャッジ日時 2025-01-29 17:52:39
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std ;
#define fast_input_output ios::sync_with_stdio(false); cin.tie(nullptr);
typedef long long ll ;
typedef long double ld ;
typedef pair<ll,ll> P ;
typedef tuple<ll,ll,ll> TP ;
#define chmin(a,b) a = min(a,b)
#define chmax(a,b) a = max(a,b)
#define bit_count(x) __builtin_popcountll(x)
#define gcd(a,b) __gcd(a,b)
#define lcm(a,b) a / gcd(a,b) * b
#define rep(i,n) for(int i = 0 ; i < n ; i++)
#define rrep(i,a,b) for(int i = a ; i < b ; i++)
#define endl "\n"

const int mod = 1000000007 ;
ll powmod(ll x , ll n){
    ll res = 1 ;
    while(n > 0){
        if(n & 1) (res *= x) %= mod ;
        (x *= x) %= mod ;
        n >>= 1 ;
    }
    return res ;
}

int t ;

void solve(){
    ll n , m ;
    cin >> n >> m ;
    ll v = powmod(n,m/2) ;
    ll u = powmod(n,m) ;
    ll sv = v * (v + 1) % mod * (2 * v + 1) % mod * powmod(6,mod-2) % mod ;
    ll su = u * (u + 1) % mod * powmod(2,mod-2) % mod ;
    cout << (su - sv + mod) % mod << endl ;
}

int main(){
    fast_input_output
    cin >> t ;
    rep(i,t) solve() ;
}
0