結果
問題 | No.575 n! / m / m / m... |
ユーザー |
|
提出日時 | 2022-06-28 22:22:13 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,405 bytes |
コンパイル時間 | 2,378 ms |
コンパイル使用メモリ | 194,864 KB |
最終ジャッジ日時 | 2025-01-30 01:32:38 |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 9 WA * 14 |
ソースコード
#include<bits/stdc++.h>using namespace std;#define rep(i,n) for(ll i=0;i<n;i++)#define repl(i,l,r) for(ll i=(l);i<(r);i++)#define per(i,n) for(ll i=(n)-1;i>=0;i--)#define perl(i,r,l) for(ll i=r-1;i>=l;i--)#define fi first#define se second#define ins insert#define pqueue(x) priority_queue<x,vector<x>,greater<x>>#define all(x) (x).begin(),(x).end()#define CST(x) cout<<fixed<<setprecision(x)#define rev(x) reverse(x);using ll=long long;using vl=vector<ll>;using vvl=vector<vector<ll>>;using pl=pair<ll,ll>;using vpl=vector<pl>;using vvpl=vector<vpl>;const ll MOD=1000000007;const ll MOD9=998244353;const int inf=1e9+10;const ll INF=4e18;const ll dy[8]={1,0,-1,0,1,1,-1,-1};const ll dx[8]={0,1,0,-1,1,-1,1,-1};template <typename T> inline bool chmax(T &a, T b) {return ((a < b) ? (a = b, true) : (false));}template <typename T> inline bool chmin(T &a, T b) {return ((a > b) ? (a = b, true) : (false));}int main(){ll n;cin >> n;ll m;cin >> m;ll cnt=0;ll p=n;while(p){cnt+=p/m;p/=m;}CST(10);auto f=log10l(expl(1))*((logl(n)-1)*n+logl(acosl(-1)*2*n)/2)-log10l(m)*cnt;if(n<=1000){f=0;for(ll i=1;i<=n;i++)f+=log10l(i);rep(i,cnt)f-=log10l(m);}ll z=floorl(f);//cout << z << endl;long double ans=0;for(long double i=1.0000;i<10;i+=0.0001){if(log10l(i)<=f-z)ans=i;}cout << ans << "e" << z << endl;}