結果

問題 No.8072 Sum of sqrt(x)
ユーザー 👑 hos.lyric
提出日時 2020-09-12 16:00:45
言語 D
(dmd 2.109.1)
結果
WA  
実行時間 -
コード長 2,487 bytes
コンパイル時間 2,618 ms
コンパイル使用メモリ 158,684 KB
実行使用メモリ 22,244 KB
最終ジャッジ日時 2024-06-22 09:02:45
合計ジャッジ時間 13,662 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1 WA * 2
other AC * 3 WA * 2 TLE * 1 -- * 21
権限があれば一括ダウンロードができます

ソースコード

diff # 

import std.conv, std.functional, std.range, std.stdio, std.string;
import std.algorithm, std.array, std.bigint, std.bitmanip, std.complex, std.container, std.math, std.mathspecial, std.numeric, 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 readLong() {
  int s = 1, c;
  for (; ; ) {
    c = getchar();
    if ('0' <= c && c <= '9') break;
    if (c == -1) throw new EOFException;
    if (c == '-') s = -1;
  }
  long x = c - '0';
  for (; ; ) {
    c = getchar();
    if (!('0' <= c && c <= '9')) return s * x;
    x = x * 10 + (c - '0');
  }
}

// floor(sqrt(a))
long floorSqrt(long a) {
  import core.bitop : bsr;
  import std.algorithm : min;
  long b = a, x = 0, y = 0;
  for (int e = bsr(a) & ~1; e >= 0; e -= 2) {
    x <<= 1;
    y <<= 1;
    if (b >= (y | 1) << e) {
      b -= (y | 1) << e;
      x |= 1;
      y += 2;
    }
  }
  return x;
}


void main() {
  try {
    for (; ; ) {
      const N = cast(int)(readLong());
      auto X = new long[N];
      foreach (i; 0 .. N) {
        X[i] = readLong();
      }
      
      long upp;
      real low = 0.0L;
      foreach (i; 0 .. N) {
        if (X[i] > 0) {
          const k = floorSqrt(X[i]);
          upp += k;
          /*
            sqrt(X[i]) - k
            = (X[i] - k^2) / (sqrt(X[i]) + k)
          */
          low += (X[i] - k^^2) / (sqrt(cast(real)(X[i])) + k);
        }
        const ans = format("%.18f", upp + low);
        foreach (c; ans[0 .. min($, 16)]) {
          putchar(c);
        }
        putchar('\n');
      }
    }
  } catch (EOFException e) {
  }
}
0