結果
問題 | No.575 n! / m / m / m... |
ユーザー |
|
提出日時 | 2021-06-19 18:45:38 |
言語 | PyPy3 (7.3.15) |
結果 |
RE
|
実行時間 | - |
コード長 | 1,311 bytes |
コンパイル時間 | 316 ms |
コンパイル使用メモリ | 82,432 KB |
実行使用メモリ | 66,048 KB |
最終ジャッジ日時 | 2024-06-22 22:15:59 |
合計ジャッジ時間 | 2,671 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 RE * 1 |
other | AC * 17 RE * 6 |
ソースコード
import sys input = lambda : sys.stdin.readline().rstrip() sys.setrecursionlimit(2*10**5+10) write = lambda x: sys.stdout.write(x+"\n") debug = lambda x: sys.stderr.write(x+"\n") writef = lambda x: print("{:.12f}".format(x)) n,m = list(map(int, input().split())) import math def factor(n, m=None): # mを与えると、高々その素因数まで見て、残りは分解せずにそのまま出力する f = {} tmp = n M = int(-(-n**0.5//1))+1 if m is not None: M = min(m+1, M) for i in range(2, M+1): if tmp<i: break if tmp%i==0: cnt=0 while tmp%i==0: cnt+=1 tmp //= i f[i] = cnt if tmp!=1: f[tmp] = 1 if not f: f[n] = 1 return f def divnum(v,k): """v!がkで割れる回数 """ assert v>=k ans = 0 kk = k while kk<=v: ans += v//kk kk *= k return ans f2 = factor(m) val = float("inf") for k in f2: val = min(val, divnum(n,k)//f2[k]) if n<10**6: ans = sum((math.log10(i) for i in range(1, n+1))) else: ans = math.log(2*math.pi*n)/2 + n*math.log(n) - n ans /= math.log(10) ans -= val * math.log10(m) v = math.pow(10, ans%1) # print(v) e = ans//1 print(f"{v}e{int(e)}") # print("{:e}".format(p))