# 挑战程序竞赛系列（85）：3.6极限情况（2）

## 挑战程序竞赛系列（85）：3.6极限情况（2）

```import java.io.BufferedReader;
import java.io.File;
import java.io.FileInputStream;
import java.io.IOException;
import java.io.InputStream;
import java.io.PrintWriter;
import java.util.Arrays;
import java.util.StringTokenizer;

public class Main{

String INPUT = "./data/judge/201709/P1418.txt";

public static void main(String[] args) throws IOException {
new Main().run();
}

static final double PI  = Math.acos(-1);
static final double EPS = 5E-13;
static final int MAX_N  = 102;

class P{

double x;
double y;

P(double x, double y){
this.x = x;
this.y = y;
}
}

P[]   o = new P[MAX_N];                 // 圆心坐标
double[] r = new double[MAX_N];         // 圆半径
boolean[] visible = new boolean[MAX_N]; // 对应圆是否可见
int N;

double[] angle = new double[2 * MAX_N];
int tot;

double distance(P a, P b) {
double dx = a.x - b.x;
double dy = a.y - b.y;
return Math.sqrt(dx * dx + dy * dy);
}

double norm(double ang) {
while (ang < 0) {
ang += 2 * PI;
}
while (ang > 2 * PI) {
ang -= 2 * PI;
}
return ang;
}

void solve() {
Arrays.fill(visible, false);
for (int i = 0; i < N; ++i) {
tot = 0;
angle[tot++] = 0;
angle[tot++] = 2 * PI;

for (int j = 0; j < N; ++j) {
if (i == j) continue;
double d = distance(o[i], o[j]);
if (r[i] + r[j] < d || d < r[i] - r[j] || d < r[j] - r[i]) continue; // 包含或者不相交
double phi = Math.atan2(o[j].y - o[i].y, o[j].x - o[i].x);
double the = Math.acos((r[i] * r[i] + d * d - r[j] * r[j]) / (2 * r[i] * d));
angle[tot++] = norm(phi - the);
angle[tot++] = norm(phi + the);
}

Arrays.sort(angle, 0, tot);

for (int j = 0; j < tot - 1; ++j) {
double mid = (angle[j] + angle[j + 1]) / 2;

double nx = 0;
double ny = 0;

for (int t = -1; t < 2; t += 2) {
nx = o[i].x + (r[i] + t * EPS) * Math.cos(mid);
ny = o[i].y + (r[i] + t * EPS) * Math.sin(mid);

P np = new P(nx, ny);

int k = 0;
for (k = N - 1; k >= 0; --k) {
if (distance(np, o[k]) < r[k]) {
break;
}
}
if (k != -1)
visible[k] = true;
}
}
}

int ans = 0;
for (int i = 0; i < N; ++i) {
if (visible[i]) ans ++;
}
out.println(ans);
}

while (true) {
N = ni();
if (N == 0) break;

for (int i = 0; i < N; ++i) {
o[i] = new P(nd(), nd());
r[i] = nd();
}

solve();
}
}

FastScanner in;
PrintWriter out;

void run() throws IOException {
boolean oj;
try {
oj = ! System.getProperty("user.dir").equals("F:\\java_workspace\\leetcode");
} catch (Exception e) {
oj = System.getProperty("ONLINE_JUDGE") != null;
}

InputStream is = oj ? System.in : new FileInputStream(new File(INPUT));
in = new FastScanner(is);
out = new PrintWriter(System.out);
long s = System.currentTimeMillis();
out.flush();
if (!oj){
System.out.println("[" + (System.currentTimeMillis() - s) + "ms]");
}
}

public boolean more(){
return in.hasNext();
}

public int ni(){
return in.nextInt();
}

public long nl(){
return in.nextLong();
}

public double nd(){
return in.nextDouble();
}

public String ns(){
return in.nextString();
}

public char nc(){
return in.nextChar();
}

class FastScanner {
StringTokenizer st;
boolean hasNext;

public FastScanner(InputStream is) throws IOException {
hasNext = true;
}

public String nextToken() {
while (st == null || !st.hasMoreTokens()) {
try {
} catch (Exception e) {
hasNext = false;
return "##";
}
}
return st.nextToken();
}

String next = null;
public boolean hasNext(){
next = nextToken();
return hasNext;
}

public int nextInt() {
if (next == null){
hasNext();
}
String more = next;
next = null;
return Integer.parseInt(more);
}

public long nextLong() {
if (next == null){
hasNext();
}
String more = next;
next = null;
return Long.parseLong(more);
}

public double nextDouble() {
if (next == null){
hasNext();
}
String more = next;
next = null;
return Double.parseDouble(more);
}

public String nextString(){
if (next == null){
hasNext();
}
String more = next;
next = null;
return more;
}

public char nextChar(){
if (next == null){
hasNext();
}
String more = next;
next = null;
return more.charAt(0);
}
}
}```

98 篇文章33 人订阅

0 条评论

## 相关文章

### 等价类方法和边界值分析方法

NextDate是一个有三个变量（月份、日期和年）的函数。函数返回输入日期后面的那个日期。变量月份、日期和年都是整数值，并满足以下条件：

622

34710

### [五]java函数式编程归约reduce概念原理 stream reduce方法详解 reduce三个参数的reduce方法如何使用

java8 流相关的操作中,我们把它理解 "累加器",之所以加引号是因为他并不仅仅是加法

1593

1876

### Python：lambda表达式的两种应用场景

python书写简单，功能强大， 迅速发展成为 AI ，深度学习的主要语言。介绍Python中的lambda表达式，注意到，它只是一个表达式，不是语句啊。

941

892

792

### CUDA PTX ISA阅读笔记（二）

8. 第八章 指令集 这一章占了整个手册的一大半(百十来页吧)，主要介绍各种指令，虽然页数很多，但是大多数指令都很简单。 8.1. 指令的形式和语义描述 这章就...

3775

### 每周学点大数据 | No.8基础数据结构之线性表

No.8期 基础数据结构之线性表 Mr. 王：为了以后的知识描述方便，这里简单介绍一下数据结构的概念。数据结构是一个广泛存在于计算机科学中的概念。曾经有一位计...

33111

852