[CSP-S模拟测试]:Walker(数学)

题目传送门(内部题86)


输入格式

第一行$n$
接下来$n$行,每行四个浮点数,分别表示变换前的坐标和变换后的坐标


输出格式

第一行浮点数$\theta$以弧度制表示
第二行浮点数$scale$
第三行两个浮点数$d_x,d_y$
我们将用$SPJ$以$10^{-3}$的绝对误差来判断变换结果是否正确,建议输出$10$位小数以上。


样例

样例输入1:

5
0 0 -1 1
0 1 -2 1
1 0 -1 2
1 1 0 0
2 1 1 0

样例输出1:

1.5707963268
1
-1 1

样例输入2:

5
0 0 1 1
0 1 1 2
1 0 2 1
1 1 0 0
2 1 1 0

样例输出2:

0
1
1 1


数据范围与提示

对于$30\%$的数据,保证存在可行的答案$\theta=0,scale=1$
对于另$30\%$的数据,$n\leqslant 500$
对于所有数据$n\leqslant 100,000$,坐标属于$[-10^9,10^9]$


题解

正解是随机化额……

不过还是有正确性的。

如果我们知道两组$(x_i,y_i)$,那么就可以解出$\theta,scale,d_x,d_y$了,然后将解出来的这四个值带入检验即可。

因为至少有一半是正确的,所以每次可以解出的概率就是$\frac{1}{4}$,解不出来的概率就是$\frac{3}{4}$,那么解$k$次,解不出来的概率就是$\frac{1}{4}^k$,$k=50$时概率已经$<10^{-5}$了。

对于官方题解里说的高斯消元,因为就两组,没有必要。

实在不行多交几遍(其实我就脸黑了……)

时间复杂度:$\Theta($玄学$)$。

期望得分:$100$分。

实际得分:$0\sim 100$分(脸黑别怪我)。


代码时刻

#include<bits/stdc++.h>
using namespace std;
const double eps=1e-5;
struct rec{double x0,y0,x2,y2;}e[100001];
pair<double,double> d[100001];
int n;
double dx,dy,Cos,Sin,scale,theta;
void solve(int x,int y)
{
	double d1=e[x].x0-e[y].x0;
	double d2=e[x].y0-e[y].y0;
	double d3=e[x].x2-e[y].x2;
	double d4=e[x].y2-e[y].y2;
	double COS=(d1*d3+d2*d4)/(d1*d1+d2*d2);
	double SIN=(d4-d2*COS)/d1;
	theta=atan(SIN/COS);
	dx=e[x].x2-e[x].x0*COS+e[x].y0*SIN;
	dy=e[x].y2-e[x].x0*SIN-e[x].y0*COS;
	Cos=cos(theta);
	Sin=sin(theta);
	scale=SIN/Sin;
	if(scale<0.0)
	{
		theta+=M_PI;
		Cos=cos(theta);
		Sin=sin(theta);
		scale=SIN/Sin;
	}
}
bool judge()
{
	int sum=0;
	for(int i=1;i<=n;i++)
	{
		double nowx=e[i].x0*Cos-e[i].y0*Sin;
		double nowy=e[i].x0*Sin+e[i].y0*Cos;
		nowx=nowx*scale+dx;nowy=nowy*scale+dy;
		if(fabs(nowx-e[i].x2)<eps&&abs(nowy-e[i].y2)<eps)sum++;
	}
	return (2*sum>=n);
}
int main()
{
	srand(time(NULL));
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
		scanf("%lf%lf%lf%lf",&e[i].x0,&e[i].y0,&e[i].x2,&e[i].y2);
		d[i]=make_pair(e[i].x2-e[i].x0,e[i].y2-e[i].y0);
	}
	sort(d+1,d+n+1);
	int sum=0;
	for(int i=1;i<=n;i++)
	{
		if(fabs(d[i].first-d[i-1].first)>eps||fabs(d[i].second-d[i-1].second)>eps)sum=0;
		sum++;
		if(2*sum>=n)
		{
			puts("0\n1");
			printf("%.10lf %.10lf",d[i].first,d[i].second);
			return 0;
		}
	}
	while(1)
	{
		int x,y;
		do
		{
			x=rand()%n+1;
			y=rand()%n+1;
		}while(x==y);
		solve(x,y);
		if(judge())
		{
			printf("%.10lf\n%.10lf\n%.10lf %.10lf",theta,scale,dx,dy);
			return 0;
		}
	}
	return 0;
}


rp++

原文地址:https://www.cnblogs.com/wzc521/p/11736969.html

时间: 2024-11-06 07:14:57

[CSP-S模拟测试]:Walker(数学)的相关文章

模拟测试(vj)

做这份模拟测试,已经崩溃了,英文看不懂,题意理解错.到结束了只a了第一题,人生陷入了低谷,于是花了一天的时间终于把不会的弄明白了,在这里写一份总结~ T1,简单的模拟,如果打枪打中一支鸟,将这个位置设为0,并向两边扩散,注意这个位置一定要有鸟. 代码~ #include<bits/stdc++.h> using namespace std; int a[30000]; int n,m; int main() { cin>>n; for(int i=1;i<=n;i++) ci

Android单元测试与模拟测试详解

测试与基本规范 为什么需要测试? 为了稳定性,能够明确的了解是否正确的完成开发. 更加易于维护,能够在修改代码后保证功能不被破坏. 集成一些工具,规范开发规范,使得代码更加稳定( 如通过 phabricator differential 发diff时提交需要执行的单元测试,在开发流程上就可以保证远端代码的稳定性). 2. 测什么? 一般单元测试: 列出想要测试覆盖的异常情况,进行验证. 性能测试. 模拟测试: 根据需求,测试用户真正在使用过程中,界面的反馈与显示以及一些依赖系统架构的组件的应用测

微信在线信息模拟测试工具(基于Senparc.Weixin.MP)

目前为止似乎还没有看到过Web版的普通消息测试工具(除了官方针对高级接口的),现有的一些桌面版的几个测试工具也都是使用XML直接请求,非常不友好,我们来尝试做一个“面向对象”操作的测试工具. 测试工具在线DEMO:http://weixin.senparc.com/SimulateTool Senparc.Weixin.MP是一个开源的微信SDK项目,地址:https://github.com/JeffreySu/WeiXinMPSDK (其中https://github.com/Jeffrey

css Hack,用IE11模拟测试的,条件注释要找真IE去测,模拟的无效

<!DOCTYPE html> <!--[if lt IE 7 ]> <html class="ie6 ie"> <![endif]--> <!--[if IE 7 ]> <html class="ie7 ie"> <![endif]--> <!--[if IE 8 ]> <html class="ie8 ie"> <![endif]

2016年上半年软考真题在线模拟测试,提前知晓你的成绩

2016年上半年软考于5月21日结束了,考试完想知道结果的急切心理,几乎每个经历过上学时代的人都能深刻体验到.如果你想知道你考的怎么样,如果你想要提前知道你的成绩,如果你想知道你哪个地方出错了,如果你想知道你哪个地方知识掌握的不够想要更深入的去理解,那就来希赛软考学院吧!希赛软考学院提供2016年上半年软考真题在线模拟测试,有标准的参考答案,有专业老师的解析视频,让你提前知晓你的成绩,让你再次巩固学习. 希赛授课专家介绍 张友生,计算机应用技术博士,软考培训教程系列丛书主编,考试指定教材<系统分

Mock 模拟测试简介及 Mockito 使用入门

Mock 是什么 mock 测试就是在测试过程中,对于某些不容易构造或者不容易获取的对象,用一个虚拟的对象来创建以便测试的测试方法.这个虚拟的对象就是mock对象.mock对象就是真实对象在调试期间的代替品. 简单的看一张图 我们在测试类 A 时,类 A 需要调用类 B 和类 C,而类 B 和类 C 又需要调用其他类如 D.E.F 等,假如类 D.E.F 构造很耗时又或者调用很耗时的话是非常不便于测试的(比如是 DAO 类,每次访问数据库都很耗时).所以我们引入 Mock 对象. 如上图,我们将

【模拟】【数学】CSU 1803 2016 (2016湖南省第十二届大学生计算机程序设计竞赛)

题目链接: http://acm.csu.edu.cn/OnlineJudge/problem.php?id=1803 题目大意: 给定n,m(n,m<=109)1<=i<=n,1<=j<=m,求i*j%2016=0的方案数. 题目思路: [模拟][数学] 按照%2016的余数分类.每增加一个2016就又多一种方案.统计是2016的几倍,根据余数分类.最后枚举i,j的余数即可求解. 1 // 2 //by coolxxx 3 //#include<bits/stdc++

asp.net中模拟测试smtp发邮件

最近在编程人生里要测试一个会员邮件的功能,就写了下面的代码. 在asp.net 中,有时要测试发信SMTP,但如果在单元测试中,如果没方便好用的 smtp怎么办,其实还是有办法模拟的,下面讲解下: 在web.config 中设置 <system.net>   <mailSettings>      <smtp deliveryMethod="SpecifiedPickupDirectory">          <specifiedPickup

2018冬令营模拟测试赛(三)

2018冬令营模拟测试赛(三) [Problem A]摧毁图状树 试题描述 输入 见"试题描述" 输出 见"试题描述" 输入示例 见"试题描述" 输出示例 见"试题描述" 数据规模及约定 见"试题描述" 题解 这题没想到贪心 QwQ,那就没戏了-- 贪心就是每次选择一个最深的且没有被覆盖的点向上覆盖 \(k\) 层,因为这个"最深的没有被覆盖的点"不可能再有其它点引出的链覆盖它了,而它又