結果

問題 No.413 +5,000,000pts
コンテスト
ユーザー yuppe19 😺
提出日時 2016-06-06 17:07:54
言語 Java
(openjdk 26.0.2.1 + ACL)
コンパイル:
javac -J-Duser.language=en -encoding UTF8 -cp /opt/aclib/ac_library.jar _filename_
実行:
java -ea -Xmx700m -Xss256M -DONLINE_JUDGE=true -cp .:/opt/aclib/ac_library.jar _class_
結果
AC  
実行時間 1,071 ms / 5,000 ms
+ 364µs
コード長 1,213 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,404 ms
コンパイル使用メモリ 86,492 KB
実行使用メモリ 55,816 KB
最終ジャッジ日時 2026-09-14 10:10:09
合計ジャッジ時間 4,442 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 1
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import java.math.*;
import java.util.*;

import static java.lang.Math.*;
import static java.math.BigDecimal.*;

public class Main {
  private static final long MAX_VALUE = (long)pow(10, 18);
  private static final int N = 100000;
  private static final BigDecimal TWO  = BigDecimal.valueOf(2),
                                  FOUR = BigDecimal.valueOf(4);

  private long calc(long d) {
    return (long)((-1 + sqrt(1+4*d)) / 2.0);
  }

  private BigDecimal mysqrt(BigDecimal a) {
    BigDecimal x0 = ZERO,
               x1 = BigDecimal.valueOf(sqrt(a.doubleValue()));
    while(!x0.equals(x1)) {
      x0 = x1;
      x1 = a.divide(x0, 28, ROUND_HALF_UP).add(x0).divide(TWO, 28, ROUND_HALF_UP);
    }
    return x1;
  }

  private long correct(long d) {
    return mysqrt(BigDecimal.valueOf(d).multiply(FOUR).add(BigDecimal.ONE)).subtract(BigDecimal.ONE).divide(TWO).longValue();
  }

  private void solve() {
    for(long k=912345678L, cnt=0; cnt<N; ++k) {
      long d = k*k + k - 1;
      if(!(1 <= d && d <= MAX_VALUE)) { continue; }
      if(calc(d) != correct(d)) {
        ++cnt;
        System.out.println(d);
      }
    }
  }

  public static void main(String[] args) {
    new Main().solve();
  }
}
0