結果

問題 No.575 n! / m / m / m...
ユーザー kotatsugame
提出日時 2017-10-09 21:46:33
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 12 ms / 2,000 ms
コード長 874 bytes
コンパイル時間 621 ms
コンパイル使用メモリ 77,884 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-17 07:05:00
合計ジャッジ時間 1,728 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:20:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type]
   20 | main()
      | ^~~~

ソースコード

diff #
プレゼンテーションモードにする

#include<iostream>
#include<iomanip>
#include<cmath>
#include<algorithm>
using namespace std;
long n,m,c=1e15,d;
double N,M;
double mypow(double a,long b)
{
    if(b==0)return 1;
    else
    {
        while(a>=10)d+=b,a/=10;
        double ret=mypow(a*a,b/2)*(b%2?a:1);
        while(ret>=10)d++,ret/=10;
        while(ret<1)d--,ret*=10;
        return ret;
    }
}
main()
{
    cin>>n>>m;N=n,M=m;
    for(long i=2;i*i<=m;i++)
    {
        if(m%i==0)
        {
            long cnt=0,a=n/i;
            while(m%i==0)m/=i,cnt++;
            for(long j=i;~-n/j+1>=i;a+=n/j)j*=i;
            c=min(c,a/cnt);
        }
    }
    if(m>1)
    {
        long a=n/m;
        for(long j=m;~-n/j+1>=m;a+=n/j)j*=m;
        c=min(c,a);
    }
    N++;
    double ori=(N+1/(12*N-1/(10*N)))/M_E;
    double f=sqrtl(2*M_PI/N)*mypow(ori,N);
    while(f<1)d--,f*=10;
    long nowd=d;d=0;
    M=mypow(M,c);
    d=nowd-d;
    f/=M;
    while(f<1)d--,f*=10;
    while(f>=10)d++,f/=10;
    while(d<0)d++,f/=10;
    cout<<fixed<<setprecision(15)<<f<<"e"<<d<<endl;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0