#69222 (C++11) No.109 N! mod M

• 問題: No.109 N! mod M
• ユーザー: piyoko_212
• 提出日時: 2015-12-28 19:44:57
• 言語: C++11 (gcc 11.4.0)
• 結果: AC
• 実行時間: 1,186 ms / 5,000 ms
• コード長: 791 bytes
• コンパイル時間: 338 ms
• コンパイル使用メモリ: 39,400 KB
• 実行使用メモリ: 4,380 KB
• 最終ジャッジ日時: 2023-09-04 05:32:32
• 合計ジャッジ時間: 2,336 ms
• ジャッジサーバーID （参考情報）: judge15 / judge13

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 4,376 KB
testcase_01	AC	77 ms; 4,380 KB
testcase_02	AC	98 ms; 4,380 KB
testcase_03	AC	1 ms; 4,380 KB
testcase_04	AC	22 ms; 4,380 KB
testcase_05	AC	1,186 ms; 4,380 KB
testcase_06	AC	5 ms; 4,380 KB
testcase_07	AC	11 ms; 4,376 KB
testcase_08	AC	1 ms; 4,376 KB

ソースコード

#include<stdio.h>
#include<algorithm>
using namespace std;
long long extgcd(long long a,long long b,long long&x,long long&y){
	for(long long u=y=1,v=x=0;a;){
		long long q=b/a;swap(x-=q*u,u);swap(y-=q*v,v);swap(b-=q*a,a);
	}
	return b;
}
long long mod_inv(long long a,long long m){
	long long x,y;
	extgcd(a,m,x,y);
	return (m+x%m)%m;
}
int main(){
	int T;
	scanf("%d",&T);
	while(T--){
		int a,b;scanf("%d%d",&a,&b);
		if(a>=b){
			printf("0\n");continue;
		}
		if(a<300000){
			long long ret=1%b;
			for(int i=1;i<=a;i++)ret=ret*i%b;
			printf("%lld\n",ret);continue;
		}
		bool hp=false;
		for(int i=2;i*i<=b;i++){
			if(b%i==0)hp=true;
		}
		if(hp){
			printf("0\n");continue;
		}
		long long now=b-1;
		for(int i=b-1;i>a;i--){
			now=now*mod_inv(i,b)%b;
		}
		printf("%lld\n",now);
	}
}