BZOJ 2508 簡單題 數學算法
題目大意：維護一個平面，支持三種操作：
0.加入一條直線(給的是兩點式)
1.刪除一條直線
2.詢問到所有直線距離平方和最小的點
題解見 www.Bkjia.com
我只是貼代碼供參考的- -
注意我的abcdef和題解設的不一樣- -
這簡單題WA了兩頁- -
#include
#include
#include
#include
#include
#define M 120100
#define EPS 1e-7
using namespace std;
struct Line{
double A,B,C;
Line() {}
Line(double x1,double y1,double x2,double y2)
{
A=y1-y2;
B=x2-x1;
C=x1*y2-x2*y1;
}
}lines[M];
int n,tot;
double a,b,c,d,e,f;
//ans=a*x^2+b*y^2+c+d*x*y+e*x+f*y
int main()
{
int i,p,x;
double x1,y1,x2,y2;
double A,B,C,temp;
cin>>n;
for(i=1;i<=n;i++)
{
scanf(%d,&p);
switch(p)
{
case 0:
scanf(%lf%lf%lf%lf,&x1,&y1,&x2,&y2);
lines[++tot]=Line(x1,y1,x2,y2);
A=lines[tot].A;
B=lines[tot].B;
C=lines[tot].C;
temp=A*A+B*B;
a+=A*A/temp;
b+=B*B/temp;
c+=C*C/temp;
d+=2*A*B/temp;
e+=2*A*C/temp;
f+=2*B*C/temp;
break;
case 1:
scanf(%d,&x);
A=lines[x].A;
B=lines[x].B;
C=lines[x].C;
temp=A*A+B*B;
a-=A*A/temp;
b-=B*B/temp;
c-=C*C/temp;
d-=2*A*B/temp;
e-=2*A*C/temp;
f-=2*B*C/temp;
break;
case 2:
if(fabs(d*d-4*a*b)EPS||fabs(d)>EPS) A=2*a,B=d,C=e;
else if(fabs(b)>EPS||fabs(d)>EPS) A=d,B=2*b,C=f;
else
{
printf(%.2lf
,0.0);
break;
//這裡我不明白為什麼 望高人指點
}
if(fabs(B)
