結果

問題 No.575 n! / m / m / m...
コンテスト
ユーザー ryoissy
提出日時 2017-10-07 15:07:15
言語 C++11
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=gnu++11 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 5 ms / 2,000 ms
コード長 870 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,093 ms
コンパイル使用メモリ 183,188 KB
実行使用メモリ 7,972 KB
最終ジャッジ日時 2026-05-12 05:56:18
合計ジャッジ時間 2,824 ms
ジャッジサーバーID
（参考情報） 		judge3_1 / tmp-judge_1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
#define MOD 1000000007LL
using namespace std;
typedef long long ll;
typedef pair<int,int> P;

ll n,m;
vector<ll> vec;
vector<int> cnts;

int main(void){
	scanf("%lld%lld",&n,&m);
	ll mp=m;
	for(ll i=2;i*i<=n;i++){
		if(mp%i==0LL){
			vec.push_back(i);
			cnts.push_back(0);
			while(mp%i==0LL){
				cnts[cnts.size()-1]++;
				mp/=i;
			}
		}
	}
	if(mp!=1LL){
		vec.push_back(mp);
		cnts.push_back(1);
	}
	ll cnt=0;
	ll np=n;
	int pp=1;
	for(int i=0;i<vec.size();i++){
		ll now=vec[i];
		ll tmp=0;
		while(1){
			tmp+=n/now;
			if(n/now<vec[i])break;
			now*=vec[i];
		}
		tmp/=(ll)cnts[i];
		if(i==0)cnt=tmp;
		else cnt=min(cnt,tmp);
	}
	long double sum=lgammal(n+1LL);
	sum/=log(10);
	long double nex=log10l(m);
	sum-=(long double)nex*cnt;
	ll resp=sum;
	double rest=(long double)sum-resp;
	printf("%.5fe%lld\n",pow(10,rest),resp);
	return 0;
}
0