結果

問題 No.644 G L C C D M
ユーザー rickythetarickytheta
提出日時 2018-02-02 23:52:53
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 77 ms / 2,000 ms
コード長 1,501 bytes
コンパイル時間 1,493 ms
コンパイル使用メモリ 158,196 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-10 01:37:05
合計ジャッジ時間 2,001 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,812 KB
testcase_01 AC 2 ms
6,944 KB
testcase_02 AC 2 ms
6,944 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 1 ms
6,944 KB
testcase_05 AC 1 ms
6,940 KB
testcase_06 AC 1 ms
6,944 KB
testcase_07 AC 2 ms
6,944 KB
testcase_08 AC 2 ms
6,944 KB
testcase_09 AC 1 ms
6,944 KB
testcase_10 AC 2 ms
6,944 KB
testcase_11 AC 1 ms
6,944 KB
testcase_12 AC 1 ms
6,940 KB
testcase_13 AC 4 ms
6,940 KB
testcase_14 AC 2 ms
6,944 KB
testcase_15 AC 1 ms
6,940 KB
testcase_16 AC 2 ms
6,940 KB
testcase_17 AC 1 ms
6,940 KB
testcase_18 AC 1 ms
6,940 KB
testcase_19 AC 1 ms
6,940 KB
testcase_20 AC 1 ms
6,940 KB
testcase_21 AC 2 ms
6,944 KB
testcase_22 AC 77 ms
6,944 KB
testcase_23 AC 32 ms
6,940 KB
testcase_24 AC 19 ms
6,940 KB
testcase_25 AC 1 ms
6,944 KB
testcase_26 AC 2 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;

typedef int _loop_int;
#define REP(i,n) for(_loop_int i=0;i<(_loop_int)(n);++i)
#define FOR(i,a,b) for(_loop_int i=(_loop_int)(a);i<(_loop_int)(b);++i)
#define FORR(i,a,b) for(_loop_int i=(_loop_int)(b)-1;i>=(_loop_int)(a);--i)

#define DEBUG(x) cout<<#x<<": "<<x<<endl
#define DEBUG_VEC(v) cout<<#v<<":";REP(i,v.size())cout<<" "<<v[i];cout<<endl
#define ALL(a) (a).begin(),(a).end()

#define CHMIN(a,b) a=min((a),(b))
#define CHMAX(a,b) a=max((a),(b))

// mod
const ll MOD = 1000000007ll;
#define FIX(a) ((a)%MOD+MOD)%MOD

int main(){
  int n,m;
  cin>>n>>m;
  // gcd(a[0], a[1]) = m
    // after operation,
    // b[0] = gcd(a[0], a[1]) = a[0] / s
    // b[1] = lcm(a[1], a[2]) = a[1] * t
    // and next operation,
    // b2[0] = gcd(b[0], b[1]) = gcd(a[0]/s, a[1]*t) = a[0]/s = gcd(a[0], a[1])
  // a[0] = x * m
  // a[1] = y * m
  // x \perp y
  // 1 <= x*m < y*m <= n
  // answer = count (x,y) * 2 * (n-2)! mod M
  ll sum = 0;
  for(int y=2;y*m<=n;y++){
    ll tmp = y;
    ll add = y;
    for(int x=2;x*x<=tmp;x++){
      if(tmp%x==0){
        while(tmp%x==0)tmp/=x;
        add /= x;
        add *= (x-1);
      }
    }
    if(tmp>1){
      add /= tmp;
      add *= (tmp-1);
    }
    sum += add;
    sum %= MOD;
  }
  sum = sum * 2 % MOD;
  FOR(i,1,n-1){
    sum = sum * i % MOD;
  }
  printf("%lld\n",sum);
  return 0;
}
0