题解 \(by\;zj\varphi\)

发现当把 \(\rm scale×cos\theta,scale×sin\theta,dx,dy\) 当作变量时只有四个,两个方程就行。

当 \(\rm n\le 500\) 时,可以选取两组进行高斯消元,解出答案后回带。

但当 \(n\) 极大时,采用随机化的做法,每次随机选取两个,这样每次选取不正确的概率为 \(\frac{3}{4}\),选取 50 次后基本就会出答案了。

记得判断 \(\rm sin\) 的正负

Code
#include<bits/stdc++.h>
#define ri register signed
#define p(i) ++i
namespace IO{
char buf[1<<21],*p1=buf,*p2=buf;
#define gc() p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?(-1):*p1++
struct nanfeng_stream{
template<typename T>inline nanfeng_stream &operator>>(T &x) {
ri f=0;x=0;register char ch=gc();
while(!isdigit(ch)) {f|=ch=='-';ch=gc();}
while(isdigit(ch)) {x=(x<<1)+(x<<3)+(ch^48);ch=gc();}
return x=f?-x:x,*this;
}
}cin;
}
using IO::cin;
namespace nanfeng{
#define FI FILE *IN
#define FO FILE *OUT
template<typename T>inline T cmax(T x,T y) {return x>y?x:y;}
template<typename T>inline T cmin(T x,T y) {return x>y?y:x;}
typedef double db;
static const int N=1e5+7;
static const db eps=1e-6;
struct node{db x1,y1,x2,y2;}pnt[N];
int p[N],n,hl;
db mt[51][51],x1,x2,x3,x4;
inline void Guass() {
for (ri i(1);i<=4;p(i)) {
ri k=i;
for (ri j(i+1);j<=4;p(j)) if (fabs(mt[j][i])>fabs(mt[k][i])) k=j;
if (k!=i) std::swap(mt[k],mt[i]);
for (ri j(1);j<=4;p(j)) {
if (i==j) continue;
db tmp=mt[j][i]/mt[i][i];
for (ri l(1);l<=5;p(l)) mt[j][l]-=tmp*mt[i][l];
}
}
}
inline bool judge() {
ri cnt(0);
x1=mt[1][5]/mt[1][1];
x2=mt[2][5]/mt[2][2];
x3=mt[3][5]/mt[3][3];
x4=mt[4][5]/mt[4][4];
for (ri i(1);i<=n;p(i))
if (fabs(pnt[i].x1*x1-pnt[i].y1*x2+x3-pnt[i].x2)<=eps&&fabs(pnt[i].x1*x2+pnt[i].y1*x1+x4-pnt[i].y2)<=eps) p(cnt);
return cnt>=hl;
}
inline int solve() {
register db sc=sqrt(x1*x1+x2*x2);
if (x2/sc<=eps) printf("%.10lf\n",acos(-1)*2.0-acos(x1/sc));
else printf("%.10lf\n",acos(x1/sc));
printf("%.10lf\n%.10lf %.10lf\n",sc,x3,x4);
return 0;
}
inline int main() {
//FI=freopen("nanfeng.in","r",stdin);
//FO=freopen("nanfeng.out","w",stdout);
srand(unsigned(time(0)));
std::cin >> n;
hl=(n>>1)+(n&1);
for (ri i(1);i<=n;p(i)) {
scanf("%lf%lf%lf%lf",&pnt[i].x1,&pnt[i].y1,&pnt[i].x2,&pnt[i].y2);
p[i]=i;
}
if (n>500) {
std::random_shuffle(p+1,p+n+1);
for (ri i(1);i<=n;p(i)) {
ri cur1(p[i]),cur2(p[i+1]);
mt[1][1]=pnt[cur1].x1,mt[1][2]=-pnt[cur1].y1,mt[1][3]=1.0,mt[1][4]=0.0,mt[1][5]=pnt[cur1].x2;
mt[2][1]=pnt[cur1].y1,mt[2][2]=pnt[cur1].x1,mt[2][3]=0.0,mt[2][4]=1.0,mt[2][5]=pnt[cur1].y2;
mt[3][1]=pnt[cur2].x1,mt[3][2]=-pnt[cur2].y1,mt[3][3]=1.0,mt[3][4]=0.0,mt[3][5]=pnt[cur2].x2;
mt[4][1]=pnt[cur2].y1,mt[4][2]=pnt[cur2].x1,mt[4][3]=0.0,mt[4][4]=1.0,mt[4][5]=pnt[cur2].y2;
Guass();
if (judge()) return solve();
}
} else {
for (ri i(1);i<=n;p(i))
for (ri j(i+1);j<=n;p(j)) {
ri cur1(i),cur2(j);
mt[1][1]=pnt[cur1].x1,mt[1][2]=-pnt[cur1].y1,mt[1][3]=1.0,mt[1][4]=0.0,mt[1][5]=pnt[cur1].x2;
mt[2][1]=pnt[cur1].y1,mt[2][2]=pnt[cur1].x1,mt[2][3]=0.0,mt[2][4]=1.0,mt[2][5]=pnt[cur1].y2;
mt[3][1]=pnt[cur2].x1,mt[3][2]=-pnt[cur2].y1,mt[3][3]=1.0,mt[3][4]=0.0,mt[3][5]=pnt[cur2].x2;
mt[4][1]=pnt[cur2].y1,mt[4][2]=pnt[cur2].x1,mt[4][3]=0.0,mt[4][4]=1.0,mt[4][5]=pnt[cur2].y2;
Guass();
if (judge()) return solve();
}
}
return 0;
}
}
int main() {return nanfeng::main();}

NOIP 模拟 $32\; \rm Walker$的更多相关文章

  1. NOIP 模拟 $32\; \rm Six$

    题解 二维状压. 第一维直接压选不同质因子的方案,第二位压方案. 分两种讨论,显然一种方案最多出现两次,否则就不合法了,所以一种是出现了一次的,另一种是出现了两次的,这样可以减小状态数. 实现可以用 ...

  2. NOIP 模拟 $32\; \rm Smooth$

    题解 \(by\;zj\varphi\) 很简单的贪心题. 开 \(B\) 个队列,每个队列存最后一次乘上的数为当前队列编号的数. 每次去所有队列中队首的最小值,不用开堆,因为开堆用于将所有数排序,但 ...

  3. noip模拟32[好数学啊]

    noip模拟32 solutions 真是无语子,又没上100,无奈死了 虽然我每次都觉得题很难,但是还是有好多上100的 战神都200多了,好生气啊啊啊 从题开始变难之后,我的时间分配越来越不均匀, ...

  4. 20190902+0903合集-NOIP模拟

    一直没时间写QwQ 于是补一下. Day 1 晚饭吃的有点恶心…… $1s\,2s\,5s$ 还开 -O2 ?? 有点恐怖. T1 猛的一想: 把外面设成一个点, 向入口连一条权为排队时间的边 从出口 ...

  5. 2021.5.22 noip模拟1

    这场考试考得很烂 连暴力都没打好 只拿了25分,,,,,,,,好好总结 T1序列 A. 序列 题目描述 HZ每周一都要举行升旗仪式,国旗班会站成一整列整齐的向前行进. 郭神作为摄像师想要选取其中一段照 ...

  6. NOIP模拟

    1.要选一个{1,2,...n}的子集使得假如a和b在所选集合里且(a+b)/2∈{1,2,...n}那么(a+b)/2也在所选集合里 f[i]=2*f[i-1]-f[i-2]+g[i] g[n]:选 ...

  7. NOIP模拟3

    期望得分:30+90+100=220 实际得分:30+0+10=40 T1智障错误:n*m是n行m列,硬是做成了m行n列 T2智障错误:读入三个数写了两个%d T3智障错误:数值相同不代表是同一个数 ...

  8. 7.22 NOIP模拟7

    又是炸掉的一次考试 T1.方程的解 本次考试最容易骗分的一道题,但是由于T2花的时间太多,我竟然连a+b=c都没判..暴力掉了40分. 首先a+b=c,只有一组解. 然后是a=1,b=1,答案是c-1 ...

  9. NOIP模拟 1

    NOIP模拟1,到现在时间已经比较长了.. 那天是6.14,今天7.18了 //然鹅我看着最前边缺失的模拟1,还是终于忍不住把它补上,为了保持顺序2345重新发布了一遍.. #   用  户  名   ...

随机推荐

  1. Java和C#语法对比(转)

        Java C# 访问修饰符 public 修饰类.接口.变量.方法. 对所有类可见. 修饰类.接口.变量.方法. 对所有类可见.   internal 无. 修饰类.接口.变量.方法. 类,接 ...

  2. 基于腾讯云Serverless的HTTP服务探活函数

    本文基于 Golang 开发了一款简单易用的拨测云函数,入口函数与腾讯云 Serverless SDK 绑定.与目前腾讯云中默认的拨测函数不同的是, url-tester-func 支持非 200 响 ...

  3. 锐捷路由器 RSR20-X-28

    学习帮助视频 Ruijie#show run | in natip nat outsideip nat insideip nat pool NAT_POOL prefix-length 24ip na ...

  4. dataX windows10安装

    按照视频课程,从Github上下载文件:https://github.com/alibaba/DataX 然后将下载的压缩包解压即可,不过需要的前提Python环境是要求python2,于是在pyth ...

  5. 夏壹队——nabcd

    我们团队项目名称是TD校园通,无课表查询你作为主要功能,是一个亮眼的功能. 我们的无课表查询功能主要解决大家平时上自习还要到教学楼拍照,教室课程情况也有断更的时候,而这个功能就解决大家的这个功能,能够 ...

  6. ArrayList 深入浅出

    ArrayList 特点:按添加顺序排列.可重复.非线程安全: 底层实现:数组 扩容原理:初始化集合时,默认容量为 0,第一次添加元素时扩容为 10,容量不够时扩容为原来容量的 1.5 倍. 这里扩容 ...

  7. Spring RestTemplate 之post请求

    ●post请求:在RestTemplate中,POST请求可以通过如下三个方法来发起,但post提交方式又有两种 formData 和 payLoad,而且接口设计与传统的浏览器使用的提交方式又有差异 ...

  8. windows程序快速启动的方式:WIN键+R

    WIN键+R是windows快速启动程序的一种方式,一般能独立运行的程序都能以这种方式启动.如notepad.calc.explorer等程序. 在命令行方式下explorer加上不同的参数,会得到不 ...

  9. 微信小程序云开发-添加数据

    一.数据的添加 使用add方法添加数据 添加完成后,在数据库中查询,可以看到数据库中添加了1条数据,此时添加的数据系统自动添加了_openid 将[添加]功能写到对应的方法中 wxml页面中,点击[添 ...

  10. IO流之节点流(字符流)和数据流关闭

    ​输入流----Reader 1 public class Reader { 2 public static void main(String[] args) throws Exception { 3 ...