結果

問題 No.1966 Median of Divisors
ユーザー hiro71687khiro71687k
提出日時 2023-03-07 13:16:22
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 593 ms / 2,000 ms
コード長 773 bytes
コンパイル時間 4,196 ms
コンパイル使用メモリ 262,272 KB
実行使用メモリ 4,348 KB
最終ジャッジ日時 2023-10-18 05:29:22
合計ジャッジ時間 10,185 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,348 KB
testcase_01 AC 2 ms
4,348 KB
testcase_02 AC 4 ms
4,348 KB
testcase_03 AC 5 ms
4,348 KB
testcase_04 AC 393 ms
4,348 KB
testcase_05 AC 584 ms
4,348 KB
testcase_06 AC 593 ms
4,348 KB
testcase_07 AC 47 ms
4,348 KB
testcase_08 AC 284 ms
4,348 KB
testcase_09 AC 379 ms
4,348 KB
testcase_10 AC 381 ms
4,348 KB
testcase_11 AC 213 ms
4,348 KB
testcase_12 AC 285 ms
4,348 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
	}
}
0