HDU 6206 青岛网络赛1001 高精度 简单几何
给出的数据1e12规模,常规判点是否在圆范围内肯定要用到半径,求得过程中无法避免溢出,因此用JAVA自带的浮点大数运算,和个ZZ一样比赛中eclipse出现问题,而且太久没写JAVA语法都不清楚变量忘了申请空间/构造,还不知道为什么写不了类/函数,一个一个把运算符改成JAVA大数对象封装的函数,括号看的头疼,太傻了,应该交给IDE正常的队友写的,而且也没什么难度,JAVA大数还不用考虑什么精度,怎么简单怎么来。
import java.io.*;
import java.math.*;
import java.util.*;
import java.util.concurrent.atomic.DoubleAdder;
import java.util.zip.ZipEntry; import org.omg.PortableInterceptor.DISCARDING; public class Main {
public static BigDecimal dis(BigDecimal ax, BigDecimal ay, BigDecimal bx, BigDecimal by){
BigDecimal res = ((ax.subtract(bx)).multiply(ax.subtract(bx))).add((ay.subtract(by)).multiply(ay.subtract(by)));
return res;
}
public static BigDecimal R(BigDecimal ax, BigDecimal ay, BigDecimal bx, BigDecimal by, BigDecimal cx, BigDecimal cy) {
BigDecimal pa = by.subtract(cy);
BigDecimal pb = cy.subtract(ay);
BigDecimal pc = ay.subtract(by);
BigDecimal PA = ax.multiply(pa);
BigDecimal PB = bx.multiply(pb);
BigDecimal PC = cx.multiply(pc);
return PA.add(PB.add(PC));
}
public static BigDecimal Ux(BigDecimal ax, BigDecimal ay, BigDecimal bx, BigDecimal by, BigDecimal cx, BigDecimal cy){
BigDecimal pa = by.subtract(cy);
BigDecimal pb = cy.subtract(ay);
BigDecimal pc = ay.subtract(by);
BigDecimal PA = ((ax.multiply(ax)).add((ay.multiply(ay)))).multiply(pa);
BigDecimal PB = ((bx.multiply(bx)).add((by.multiply(by)))).multiply(pb);
BigDecimal PC = ((cx.multiply(cx)).add((cy.multiply(cy)))).multiply(pc);
return PA.add(PB.add(PC));
}
public static BigDecimal Uy(BigDecimal ax, BigDecimal ay, BigDecimal bx, BigDecimal by, BigDecimal cx, BigDecimal cy){
BigDecimal pa = cx.subtract(bx);
BigDecimal pb = ax.subtract(cx);
BigDecimal pc = bx.subtract(ax);
BigDecimal PA = ((ax.multiply(ax)).add((ay.multiply(ay)))).multiply(pa);
BigDecimal PB = ((bx.multiply(bx)).add((by.multiply(by)))).multiply(pb);
BigDecimal PC = ((cx.multiply(cx)).add((cy.multiply(cy)))).multiply(pc);
return PA.add(PB.add(PC));
}
public static void main(String[] args) {
int T;
Scanner cin = new Scanner((new BufferedInputStream(System.in)));
BigDecimal[] a = new BigDecimal[2];
BigDecimal[] b = new BigDecimal[2];
BigDecimal[] c = new BigDecimal[2];
BigDecimal[] z = new BigDecimal[2];
BigDecimal r = null;
BigDecimal x = null;
BigDecimal y = null;
BigDecimal dd = null;
BigDecimal eps = BigDecimal.valueOf(1e-20);
T = cin.nextInt();
for(int i = 0; i < T; i++)
{
a[0] = cin.nextBigDecimal(); a[1] = cin.nextBigDecimal();
b[0] = cin.nextBigDecimal(); b[1] = cin.nextBigDecimal();
c[0] = cin.nextBigDecimal(); c[1] = cin.nextBigDecimal();
z[0] = cin.nextBigDecimal(); z[1] = cin.nextBigDecimal();
BigDecimal D = R(a[0], a[1], b[0], b[1], c[0], c[1]).multiply(BigDecimal.valueOf(2));
x = Ux(a[0], a[1], b[0], b[1], c[0], c[1]).divide(D);
y = Uy(a[0], a[1], b[0], b[1], c[0], c[1]).divide(D);
//System.out.println(x+" "+y);
r = dis(a[0], a[1], x, y);
dd = dis(z[0], z[1], x, y);
//System.out.println(r+" "+dd);
if(dd.subtract(r).compareTo(eps) > 0)
System.out.println("Accepted");
else
System.out.println("Rejected"); }
}
}
HDU 6206 青岛网络赛1001 高精度 简单几何的更多相关文章
- 16年青岛网络赛 1001 I Count Two Three
题目链接:http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1001&cid=723 I Count Two ThreeTi ...
- HDU 4759 Poker Shuffle(2013长春网络赛1001题)
Poker Shuffle Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Tot ...
- HDU 4738 Caocao's Bridges (2013杭州网络赛1001题,连通图,求桥)
Caocao's Bridges Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) ...
- HDU 5880 Family View (2016 青岛网络赛 C题,AC自动机)
题目链接 2016 青岛网络赛 Problem C 题意 给出一些敏感词,和一篇文章.现在要屏蔽这篇文章中所有出现过的敏感词,屏蔽掉的用$'*'$表示. 建立$AC$自动机,查询的时候沿着$fa ...
- HDU 6212 Zuma 2017青岛网络赛 区间DP
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6212 解法:看了眼题就发现这个BZOJ 1032不是一毛一样?但是BZOJ上那是个巨坑,数据有错,原来 ...
- HDU 6215 2017Brute Force Sorting 青岛网络赛 队列加链表模拟
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6215 题意:给你长度为n的数组,定义已经排列过的串为:相邻两项a[i],a[i+1],满足a[i]&l ...
- hdu 5881 Tea (2016 acm 青岛网络赛)
原题地址:http://acm.hdu.edu.cn/showproblem.php?pid=5881 Tea Time Limit: 3000/1000 MS (Java/Others) Me ...
- HDU 5884 Sort(2016年青岛网络赛 G 二分+贪心+小优化)
好题 题意:给你n<=100000个数,每个数范围[0,1000],然后给你一个最大的代价T,每次最多合并k个数成为一个数,代价为k个数的总和.问最后合成1个数的总代价不大于T的最小k 题解:我 ...
- HDU 5886 Tower Defence(2016青岛网络赛 I题,树的直径 + DP)
题目链接 2016 Qingdao Online Problem I 题意 在一棵给定的树上删掉一条边,求剩下两棵树的树的直径中较长那的那个长度的期望,答案乘上$n-1$后输出. 先把原来那棵树的 ...
随机推荐
- 第三章 广义线性模型(GLM)
广义线性模型 前面我们举了回归和分类得到例子.在回归的例子中,$y \mid x;\theta \sim N(u,\sigma ^{2})$,在分类例子中,$y\mid x;\theta \sim ...
- YARN中用的作业调度算法:DRF(Dominant Resource Fairness)
在Mesos和YARN中,都用到了dominant resource fairness算法(DRF),它不同于hadoop基于slot-based实现的fair scheduler和capacity ...
- (六)Jmeter重要组件的执行顺序及作用域
一.Jmeter重要组件: 1)配置元件---Config Element: 用于初始化默认值和变量,以便后续采样器使用.配置元件大其作用域的初始阶段处理,配置元件仅对其所在的测试树分支有效,如,在同 ...
- C++解析(23):多态与C++对象模型
0.目录 1.多态 2.C++对象模型 2.1 使用C语言实现封装 3.继承对象模型 4.多态对象模型 4.1 使用C语言实现多态 5.小结 1.多态 面向对象中期望的行为: 根据实际的对象类型判断如 ...
- C++解析(17):操作符重载
0.目录 1.操作符重载 2.完善的复数类 3.小结 1.操作符重载 下面的复数解决方案是否可行? 示例1--原有的解决方案: #include <stdio.h> class Compl ...
- TCP的拥塞控制 (四)
TCP NewReno NewReno是在Reno的基础上,改进了Fast Recovery,主要思想是保证处于network中的packet的总量是饱和的. 在Reno算法中,一个超时会导致相 ...
- 一些noip模拟题一句话题解
Problem A: 序列 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 12 Solved: 9[Submit][Status][Web Boar ...
- Linux学习笔记一:Linux配置java环境变量
一.安装JDK: 1.创建JDK的安装目录: sudo mkdir /usr/jdk 2.将jdk-7u25-linux-x64.tar.gz拷贝至/usr/jdk目录下 sudo cp jdk-7u ...
- Spring Boot系列教程一:Eclipse安装spring-tool-suite插件
一.前言 一直使用eclipse,个人习惯选用Eclipse+spring-tool-suite进行开发,特别注意Eclipse要选用对应的spring-tool-suite进行安装,这 ...
- Classical Binary Search
Find any position of a target number in a sorted array. Return -1 if target does not exist. 与题目 Firs ...