c - Code Jam 第 1A 轮 ProblemA 的 C++ 中的二分搜索

Question

我认为比赛分析中的代码片段不正确，while(right - left >= 1)可能会陷入无限循环。所以我使用 > 而不是 >=。但它没有给出正确的答案。

但是，代码通过了断言，这意味着左边是结果，右边是界限..我认为我的代码做了正确的事情，但是......

这是我的代码:

#define _USE_MATH_DEFINES

#include<cstdio>
#include<cassert>
#include<cmath>

int main()
{
  int cases; scanf("%d", &cases);
  for(int c=1;c<=cases;c++) {
    double r, t;
    //long long t;
    scanf("%lf %lf", &r, &t);
    long long maxn = (long long)((sqrt((2*r-1)*(2*r-1)+8*t)-2*r+1)/4)+1;
    long long left=0, right=1;
    long long re= (long long)t;
    while ((2*r-1)*right+right*right*2 <= t) {
      left = right; right *= 2;
    }
    //printf("%lld\n",maxn);
    while(left+1<right) {
      long long m = left + (right - left)/2;
      double tt = (2*r-1)*1.0*m+2*m*m;
      if (tt<=t)
        left=m;
      else
        right=m;
    }
    assert((2*r-1)*1.0*right+2*right*right>t);
    assert((2*r-1)*1.0*left+2*left*left<=t);
    printf("Case #%d: %lld\n", c, left); 
  }
  return 0;
}

Answer 1

很可能，double tt = (2*r-1)*1.0*m+2*m*m; 正在失去精度。 double 只有大约 52 位整数精度。