提出 #20008581

---

ソースコード 拡げる

#include <bits/stdc++.h>
using namespace std;

#define PI acos(-1)
#define pb push_back
#define int long long int
#define double long double
#define pi pair<int, int>
#define pii pair<pi, int>
#define fir first
#define sec second
#define MAXN 2000005
#define mod 1000000007

struct pt
{
  double x, y;
  pt operator+(pt p) { return {x + p.x, y + p.y}; } // soma de pontos
  pt operator-(pt p) { return {x - p.x, y - p.y}; } // subtração de pontos
  pt operator*(double d) { return {x * d, y * d}; } // multiplicação por um double
  pt operator/(double d) { return {x / d, y / d}; } // divisão por um double
};
double dist(pt a, pt b)
{
  pt c = a - b;
  return sqrt(c.x * c.x + c.y * c.y);
}
signed main()
{
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  double x, y, raio;
  cin >> x >> y >> raio;
  pt a = {x, y};
  int ans = 0;
  for (int i = -MAXN; i <= MAXN; i++)
  {
    pt b;
    int l = ceil(y), r = MAXN, kappa = ceil(y);
    while (l < r)
    {
      int mid = (l + r + 1) >> 1;
      b = {(double)i, (double)mid};
      (dist(a, b) <= raio) ? l = mid : r = mid - 1;
    }
    b = {(double)i, (double)l};
    if (dist(a, b) <= raio)
      ans += ((l - kappa) + 1);
    l = -MAXN, r = ceil(y) - 1;
    while (l < r)
    {
      int mid = (l + r) >> 1;
      pt b = {(double)i, (double)mid};
      (dist(a, b) <= raio) ? r = mid : l = mid + 1;
    }
    b = {(double)i, (double)l};
    if (dist(a, b) <= raio)
      ans += (kappa - l);
  }
  cout << ans << endl;
  return 0;
}
// pontos com coordenadas inteiras tal que dist(i, j) <= r

提出情報

ジャッジ結果