社区讨论
89pts 求调
P1337[JSOI2004] 平衡点 / 吊打XXX参与者 1已保存回复 0
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 0 条
- 当前快照
- 1 份
- 快照标识符
- @lq7j3qyv
- 此快照首次捕获于
- 2023/12/16 11:58 2 年前
- 此快照最后确认于
- 2023/12/16 13:57 2 年前
CPP
#include<bits/stdc++.h>
using namespace std;
#define db double
const int maxn=1000;
int n;
db x[maxn+5],y[maxn+5],w[maxn+5],eps=1e-6,ansx,ansy;
db getdis(db x1,db x2,db y1,db y2){
return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
db f(db X,db Y){
db sum=0;
for(int i=1;i<=n;i++){
sum+=getdis(X,x[i],Y,y[i])*w[i];
}
return sum;
}
db solve(db X){
db l=-100000,r=100000;
while(r-l>eps){
db lmid=l+(r-l)/3,rmid=r-(r-l)/3;
if(f(X,lmid)<f(X,rmid)){
r=rmid;
}
else{
l=lmid;
}
}
ansy=l;
return f(X,l);
}
int main(){
cin>>n;
for(int i=1;i<=n;i++){
cin>>x[i]>>y[i]>>w[i];
}
db l=-100000,r=100000;
while(r-l>eps){
db lmid=l+(r-l)/3,rmid=r-(r-l)/3;
if(solve(lmid)<solve(rmid)){
r=rmid;
}
else{
l=lmid;
}
}
ansx=l;
printf("%.3lf %.3lf",ansx,ansy);
return 0;
}
回复
共 0 条回复,欢迎继续交流。
正在加载回复...