結果

問題 No.575 n! / m / m / m...
コンテスト
ユーザー hos.lyric
提出日時 2019-01-29 08:41:00
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 12 ms / 2,000 ms
コード長 2,130 bytes
コンパイル時間 2,655 ms
コンパイル使用メモリ 161,536 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-13 03:49:34
合計ジャッジ時間 4,046 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます
コンパイルメッセージ
Main.d(59): Deprecation: function `std.math.exponential.log10` is deprecated - `std.math.exponential.log10` called with argument types `(long)` matches both `log10(real)`, `log10(double)`, and `log10(float)`. Cast argument to floating point type instead.
Main.d(64): Deprecation: function `std.math.exponential.log10` is deprecated - `std.math.exponential.log10` called with argument types `(long)` matches both `log10(real)`, `log10(double)`, and `log10(float)`. Cast argument to floating point type instead.

ソースコード

diff # 

import std.conv, std.functional, std.stdio, std.string;
import std.algorithm, std.array, std.bigint, std.container, std.math, std.numeric, std.range, std.regex, std.typecons;
import core.bitop;

class EOFException : Throwable { this() { super("EOF"); } }
string[] tokens;
string readToken() { for (; tokens.empty; ) { if (stdin.eof) { throw new EOFException; } tokens = readln.split; } auto token = tokens.front; tokens.popFront; return token; }
int readInt() { return readToken.to!int; }
long readLong() { return readToken.to!long; }
real readReal() { return readToken.to!real; }

bool chmin(T)(ref T t, in T f) { if (t > f) { t = f; return true; } else { return false; } }
bool chmax(T)(ref T t, in T f) { if (t < f) { t = f; return true; } else { return false; } }

int binarySearch(alias pred, T)(in T[] as) { int lo = -1, hi = cast(int)(as.length); for (; lo + 1 < hi; ) { const mid = (lo + hi) >> 1; (unaryFun!pred(as[mid]) ? hi : lo) = mid; } return hi; }
int lowerBound(T)(in T[] as, T val) { return as.binarySearch!(a => (a >= val)); }
int upperBound(T)(in T[] as, T val) { return as.binarySearch!(a => (a > val)); }


long N, M;

long calc(long p) {
  long ret;
  long n = N;
  for (; (n /= p) > 0; ) {
    ret += n;
  }
  return ret;
}

void main() {
  try {
    for (; ; ) {
      N = readLong();
      M = readLong();
      
      long k = N;
      long m = M;
      for (long d = 2; d * d <= m; ++d) {
        if (m % d == 0) {
          long e;
          do {
            m /= d;
            ++e;
          } while (m % d == 0);
          chmin(k, calc(d) / e);
        }
      }
      if (m > 1) {
        chmin(k, calc(m) / 1);
      }
      debug {
        writeln("k = ", k);
      }
      
      real logAns = 0.0;
      if (N <= 10^^5) {
        foreach (i; 1 .. N + 1) {
          logAns += log10(i);
        }
      } else {
        logAns += log10(2.0L * PI * N) / 2.0L + log10(N / E) * N;
      }
      logAns -= log10(M) * k;
      debug {
        writeln("logAns = ", logAns);
      }
      writefln("%.10fe%.f", 10.0L^^(logAns - floor(logAns)), floor(logAns));
    }
  } catch (EOFException e) {
  }
}
0