結果

問題 No.2326 Factorial to the Power of Factorial to the...
ユーザー AckvyAckvy
提出日時 2023-05-28 15:26:14
言語 Java21
(openjdk 21)
結果
WA  
実行時間 -
コード長 2,250 bytes
コンパイル時間 3,271 ms
コンパイル使用メモリ 84,184 KB
実行使用メモリ 41,588 KB
最終ジャッジ日時 2024-06-08 07:45:45
合計ジャッジ時間 6,220 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 118 ms
41,416 KB
testcase_01 WA -
testcase_02 AC 122 ms
41,288 KB
testcase_03 AC 116 ms
39,872 KB
testcase_04 WA -
testcase_05 AC 120 ms
41,396 KB
testcase_06 AC 124 ms
41,220 KB
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 AC 119 ms
41,172 KB
testcase_11 AC 122 ms
41,348 KB
testcase_12 AC 120 ms
41,224 KB
testcase_13 WA -
testcase_14 AC 121 ms
41,504 KB
testcase_15 WA -
testcase_16 AC 118 ms
41,080 KB
testcase_17 WA -
testcase_18 WA -
testcase_19 AC 103 ms
41,580 KB
testcase_20 AC 106 ms
40,864 KB
testcase_21 AC 121 ms
41,212 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import java.util.*;

public class Main{
	static final int MOD = 1000000007;

	public static void main(String args[]){
		Scanner sc = new Scanner(System.in);
		int N = sc.nextInt();
		int P = sc.nextInt();

		int primes[] = new int[1000];
		int primeNum = 1;
		primes[0] = 2;
		for(int i = 3; i < 1000; i++){
			boolean isPrime = true;
			for(int j = 0; j < primeNum; j++){
				if(i % primes[j] == 0){
					isPrime = false;
					break;
				}
			}
			if(isPrime){
				primes[primeNum] = i;
				primeNum++;
			}
		}

		if(N < P){
			System.out.println(0);
			return;
		}

		int dum = P;
		int dividePrime[] = new int[primeNum+1];
		int productPrime[] = new int[primeNum+1];
		for(int i = 0; i < primeNum; i++){
			dividePrime[i] = 0;
			productPrime[i] = 1;
			while(dum % primes[i] == 0){
				dividePrime[i]++;
				productPrime[i] *= primes[i];
				dum /= primes[i];
			}
		}
		if(dum != 1){
			primes[primeNum] = dum;
			dividePrime[primeNum] = 1;
			productPrime[primeNum] = dum;
			primeNum++;
		}

		int maxProdInd = 0;
		for(int i = 0; i < primeNum; i++){
			if(productPrime[maxProdInd] < productPrime[i]){
				maxProdInd = i;
			}
		}

		int maxPrime = primes[maxProdInd];
		int divide = dividePrime[maxProdInd];
//		System.out.println("["+maxPrime+","+divide+"]");



		Long numP = 0l;
		for(int i = 1; i <= N; i++){
			int n = i;
			while(n % maxPrime == 0){
				numP++;
				n /= maxPrime;
			}
		}

		long amariP = numP % divide;
		numP /= divide;
//		System.out.println("["+numP+","+amariP+"]");


		long nM = 1;
		long amariNM = 1;
		for(int i = 1; i <= N; i++){
			if(i != divide){
				amariNM *= i;
			}
			nM *= i;
			nM %= MOD;
			amariNM %= MOD;
		}
//		System.out.println(nM);

		char[] S = Long.toBinaryString(nM-1).toCharArray();
		long binM = nM;
		long ans = 1;

		for(int i = S.length-1; i >= 0; i--){
//			System.out.print(S[i]);
			if(S[i] == '1'){
				ans *= binM;
				ans %= MOD;
			}
			binM *= binM;
			binM %= MOD;
		}
//		System.out.println();
		long amariAns = amariP * amariNM;
		amariAns %= MOD;
		amariAns *= ans;
		amariAns %= MOD;
//		System.out.println(amariAns);

		ans *= numP;
		ans %= MOD;
		ans *= nM;
		ans %= MOD;
//		System.out.println(ans);
		ans += amariAns;
		ans %= MOD;

		System.out.println(ans);
	}
}
0