結果

問題 No.644 G L C C D M
コンテスト
ユーザー uafr_cs
提出日時 2018-02-02 23:12:31
言語 Java
(openjdk 25.0.1)
結果
WA  
実行時間 -
コード長 1,267 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,085 ms
コンパイル使用メモリ 77,588 KB
実行使用メモリ 56,096 KB
最終ジャッジ日時 2024-12-31 08:09:43
合計ジャッジ時間 6,451 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 19 WA * 8
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import java.util.Arrays;
import java.util.LinkedList;
import java.util.Scanner;

public class Main {
	
	public static long gcd(long a, long b){
		return b == 0 ? a : gcd(b, a % b);
	}
	public static long lcm(long a, long b){
		return a / gcd(a, b) * b;
	}
	
	public static final long MOD = 1000000007;
	
	public static void main(String[] args) {
		Scanner sc = new Scanner(System.in);
		
		final int N = sc.nextInt();
		final long M = sc.nextLong();
		
		long other_fact = 1;
		for(int i = 1; i <= (N - 2); i++){
			other_fact *= i;
			other_fact %= MOD;
		}
		/*
		for(long i = M; i <= N; i += M){
			for(long j = M; j <= N; j += M){
				if(i == j || gcd(i, j) != M){ continue; }
				System.out.println(i + " " + j);
			}
		}
		*/
		
		long[] pat = new long[N + 1];
		Arrays.fill(pat, N / M);
		if(M <= N){ pat[(int)(M)]--; }
		for(long i = 2 * M; i <= N; i += M){
			final long p = N / i;
			pat[(int)(i)] -= p;
			
			for(long j = 2 * i; j <= N; j += i){ 
				pat[(int)(j)]--;
				//System.out.println("> " + i + " " + j);
			}
		}
		
		long answer = 0;
		
		for(long fst = M; fst <= N; fst += M){
			//System.out.println(fst + " : " + pat[(int)(fst)]);
			answer += (pat[(int)(fst)] * other_fact) % MOD;
			answer %= MOD;
		}
		
		System.out.println(answer);
	}
}
0