結果

問題 No.1164 GCD Products hard
ユーザー tonegawatonegawa
提出日時 2020-08-15 11:44:11
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 666 bytes
コンパイル時間 698 ms
コンパイル使用メモリ 71,328 KB
実行使用メモリ 88,224 KB
最終ジャッジ日時 2024-10-10 17:30:29
合計ジャッジ時間 57,061 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2,323 ms
64,036 KB
testcase_01 TLE -
testcase_02 AC 1,945 ms
47,872 KB
testcase_03 AC 489 ms
18,688 KB
testcase_04 AC 490 ms
19,584 KB
testcase_05 AC 2,398 ms
56,448 KB
testcase_06 TLE -
testcase_07 TLE -
testcase_08 TLE -
testcase_09 AC 1,999 ms
49,920 KB
testcase_10 AC 512 ms
17,152 KB
testcase_11 AC 2,214 ms
52,992 KB
testcase_12 TLE -
testcase_13 AC 1,758 ms
41,856 KB
testcase_14 AC 1,864 ms
49,152 KB
testcase_15 TLE -
testcase_16 AC 1,952 ms
53,504 KB
testcase_17 AC 2,349 ms
54,784 KB
testcase_18 AC 2,007 ms
50,688 KB
testcase_19 AC 536 ms
20,736 KB
testcase_20 AC 942 ms
28,672 KB
testcase_21 TLE -
testcase_22 AC 2 ms
6,820 KB
testcase_23 TLE -
testcase_24 TLE -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#define vll vector<ll>
using namespace std;
typedef long long int ll;

#define P 1000000007
inline ll mpow(ll a, ll b, ll p){
  ll ret = 1, num = a;
  for(;b;){
    if(b%2) ret = (ret*num)%p;
    num = (num*num)%p;
    b/=2;
  }
  return ret;
}
int main(int argc, char const *argv[]) {
  ll a, b, n, ans = 1;
  scanf("%lld %lld %lld", &a, &b, &n);
  ll cnt[b+1];
  for(int g=b;g>=2;g--){
    ll tmp = (b/g) - ((a-1)/g);
    cnt[g] = mpow(tmp, n, P-1);
    for(int k=2*g;k<=b;k+=g) {
      cnt[g] = (cnt[g] - cnt[k]);
      if(cnt[g]<0) cnt[g] += P-1;
    }
    ans = (ans * mpow(g, cnt[g], P))%P;
  }
  printf("%lld\n", ans);
}
0