社区讨论

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 条回复,欢迎继续交流。

正在加载回复...