看到一道很有意思的题,这个题简直有毒,是一道物理题。。。好像得用模拟退火。。。但显然我太弱了不会模拟退火,只能用正交分解暴力。。。

每次沿着力的方向走一定的距离,假如转头了,则走的步长就减小一点。

不过这里有一个坑,就是假如每次二分,就会错。。。0.7或0.9就可以。有可能走过了回不来了吧。

题面:

Description

gty又虐了一场比赛,被虐的蒟蒻们决定吊打gty。gty见大势不好机智的分出了n个分身,但还是被人多势众的蒟蒻抓住了。蒟蒻们将
n个gty吊在n根绳子上,每根绳子穿过天台的一个洞。这n根绳子有一个公共的绳结x。吊好gty后蒟蒻们发现由于每个gty重力不同,绳
结x在移动。蒟蒻wangxz脑洞大开的决定计算出x最后停留处的坐标,由于他太弱了决定向你求助。
不计摩擦,不计能量损失,由于gty足够矮所以不会掉到地上。
Input 输入第一行为一个正整数n(<=n<=),表示gty的数目。
接下来n行,每行三个整数xi,yi,wi,表示第i个gty的横坐标,纵坐标和重力。
对于20%的数据,gty排列成一条直线。
对于50%的数据,<=n<=。
对于100%的数据,<=n<=,-<=xi,yi<=
Output 输出1行两个浮点数(保留到小数点后3位),表示最终x的横、纵坐标。
Sample Input Sample Output
0.577 1.000

代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
using namespace std;
#define duke(i,a,n) for(int i = a;i <= n;i++)
#define lv(i,a,n) for(int i = a;i >= n;i--)
#define clean(a) memset(a,0,sizeof(a))
const int INF = << ;
typedef long long ll;
typedef double db;
template <class T>
void read(T &x)
{
char c;
bool op = ;
while(c = getchar(), c < '' || c > '')
if(c == '-') op = ;
x = c - '';
while(c = getchar(), c >= '' && c <= '')
x = x * + c - '';
if(op) x = -x;
}
template <class T>
void write(T x)
{
if(x < ) putchar('-'), x = -x;
if(x >= ) write(x / );
putchar('' + x % );
}
struct node
{
db x,y,power;
}p[];
int n;
bool bx = true,by = true;
db x,y;
void solve(db move)
{
db gx = ,len = ,gy = ;
duke(i,,n)
{
len = sqrt((x - p[i].x) * (x - p[i].x) + (y - p[i].y) * (y - p[i].y));
if(len == )
continue;
gx += p[i].power * (p[i].x - x) / len;
gy += p[i].power * (p[i].y - y) / len;
}
len = sqrt(gx * gx + gy * gy);
x += move * gx / len;
y += move * gy / len;
}
int main()
{
read(n);
duke(i,,n)
{
scanf("%lf%lf%lf",&p[i].x,&p[i].y,&p[i].power);
}
db move = ,tx,ty;
while()
{
tx = x;
ty = y;
solve(move);
if(abs(tx - x) < 0.000001 && abs(ty - y) < 0.000001)
break;
if(bx != (x > tx) || by != (y > ty))
{
bx = !(x > tx);
by = !(y > ty);
move = move * 0.9;
}
}
printf("%.3lf %.3lf\n",x,y);
return ;
}

最新文章

  1. C#进阶系列——WebApi 异常处理解决方案
  2. JS魔法堂:函数节流(throttle)与函数去抖(debounce)
  3. 3种不同的ContextMenu右键菜单演示
  4. Java 自动装箱与拆箱
  5. 转:Scrapy安装、爬虫入门教程、爬虫实例(豆瓣电影爬虫)
  6. RecyclerView 结合 CardView 使用(二)
  7. 用python+selenium获取XX省交通违章数据
  8. nRF51 DFU 初始化包介绍及生成工具
  9. CCNA实验(7) -- NAT
  10. QT 判断当前焦点是否在某控件上(使用hasFocus函数,毕竟QT是DIRECTUI,所以这个Focus到底代表了什么呢?)
  11. 到底什么样的企业才适合实施SAP系统?
  12. [bzoj3953] [WF2013]Self-Assembly
  13. Android音频焦点处理相关的方法
  14. Spring security oauth2 client_credentials认证 最简单示例代码
  15. MarkdownPan2 简单使用指南
  16. 机器学习之支持向量机(SVM)学习笔记
  17. Centos 6.8编译安装LNMP环境
  18. 进程池爬取并存入mongodb
  19. 基于Ardalis.GuardClauses守卫组件的拓展
  20. 20165306 预备作业3 Linux安装及学习

热门文章

  1. [Windows Server 2012] Tomcat安全加固方法
  2. js 零碎
  3. 如何描述bug
  4. 【VHDL】深度讲解二进制无符号和有符号加法处理溢出的问题
  5. AtCoder Beginner Contest 089完整题解
  6. Git使用笔记 (github为例)
  7. 网络编程:tcp、udp、socket、struct、socketserver
  8. 小白神器 - Django - 起步
  9. Django REST framework - 权限和限制
  10. 用 console.time()和 console.timeEnd() 测试你的 javascript 代码执行效率