結果

問題 No.575 n! / m / m / m...
ユーザー titia
提出日時 2025-03-26 03:55:22
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
WA  
実行時間 -
コード長 740 bytes
コンパイル時間 296 ms
コンパイル使用メモリ 12,288 KB
実行使用メモリ 10,752 KB
最終ジャッジ日時 2025-03-26 03:55:30
合計ジャッジ時間 7,143 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 21 WA * 2
権限があれば一括ダウンロードができます

ソースコード

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

n,m=map(int,input().split())
F=[]
x=m
for i in range(2,10**6):
    if x%i==0:
        F.append(i)
        while x%i==0:
            x//=i
if x!=1:
    F.append(x)
A=[]
for f in F:
    s=0
    for i in range(1,100):
        k=n//(f**i)
        if k==0:
            break
        else:
            s+=k
    A.append(s)
a=min(A)
if n<=100:
    ANS=1
    for i in range(1,n+1):
        ANS*=i
    ANS//=(m**a)
    #print(ANS)
    S=str(ANS)
    S=S[0]+"."+S[1:]
    print(S,end="")
    print("e",end="")
    print(len(S)-2)
else:
    from math import pi,e,sqrt,log10
    ANS=log10(sqrt(2*pi*n))+n*log10(n/e) - a* log10(m)
    u=int(ANS)
    ANS-=u
    S=10**ANS
    print(S,end="")
    print("e",end="")
    print(u)
    
    
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0