結果

問題 No.644 G L C C D M
ユーザー rickytheta
提出日時 2018-02-02 23:52:53
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 91 ms / 2,000 ms
コード長 1,501 bytes
コンパイル時間 1,417 ms
コンパイル使用メモリ 158,900 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-12-31 08:19:18
合計ジャッジ時間 2,512 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

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