| ||||||||||
| Online Judge | Problem Set | Authors | Online Contests | User | ||||||
|---|---|---|---|---|---|---|---|---|---|---|
| Web Board Home Page F.A.Qs Statistical Charts | Current Contest Past Contests Scheduled Contests Award Contest | |||||||||
帮忙看一下呀??WA!!!//5 -1 -1 4 3 1 1 0 3 4 0
#include<stdio.h>
#include<math.h>
#include<Limits.h>
#define fun1(a,b) a>b?b:a
#define fun2(a,b) a>b?a:b
typedef struct
{
int x;
int y;
}POINT;
POINT point[50005],a[50005],pos1[10],pos2[10], *temp, source, small;
int k;
void Find(int n)
{
int i, mult;
temp = &point[0];
for(i=1; i<n; i++)
{
if(point[i].y<temp->y||(point[i].y==temp->y&&point[i].x<temp->x))
temp=&point[i];
}
return;
}
int Mult(POINT a, POINT b)
{
return (a.x*b.y - b.x*a.y);//比较斜率//叉积.
}
double Dis(POINT k)
{
return sqrt((k.x)*(k.x)+(k.y)*(k.y));
}
int MultBigger(POINT a, POINT b, POINT c)
{
POINT p, q;
p.x = b.x-a.x;
p.y = b.y-a.y;
q.x = c.x-b.x;
q.y = c.y-b.y;
if(Mult(p,q) < 0) return 1;/*叉积运算*///斜率最小原则
else if(Mult(p,q) > 0) return 0;
else
{
if(Dis(q)>Dis(p)) return 1;
else return 0;
}
}
int Search(int n)
{
int i;
POINT temper;//source = *temp; small = *temp;
temper=source;
for(i=0; i<n; i++)
{
if(point[i].x ==INT_MAX)//已检查过的点
continue;
if (MultBigger(source,temper,point[i]))
{
temp=&point[i];
temper=point[i];
}
}//for
if(MultBigger(source, temper, small)) return 1;//回到初始点
a[k].x=temp->x;a[k].y=temp->y;k++;//回到初始点时source和temper相同.
source=*temp;
temp->x=INT_MAX;
return 0;
}
int main()
{
int n, i, c,j,m,D,D1,D2,m0,n0;
double x1,y1,x2,y2,x3,y3,x4,y4;
double min,max;
while(scanf("%d%d",&n,&c)!=EOF)
{
for(i=0;i<n;i++)
scanf("%d",&point[i].x);
for(i=0;i<n;i++)
scanf("%d",&point[i].y);
if(n==1)
{
printf("%d.000 %d.000\n",point[0].y,point[0].y);
continue;
}
Find(n); /*找出最下、最右点*/
source = *temp;
small = *temp;
//printf("%d %d\n",small.x,small.y);
k=0;
k++;
a[k].x=small.x;a[k].y=small.y;k++;
temp->x=INT_MAX;
//printf("%d\n",temp->x);
for(i=0;i<n;i++)
if(Search(n))break;/*按叉积搜索 break 时回到出发点*/
//printf("%d\n",k);
//for(i=1;i<k;i++)
//printf("%d %d\n",a[i].x,a[i].y);
a[k].x=a[1].x;a[k].y=a[1].y;
m0=n0=0;
for(i=1;i<k;i++)
{
if(c>=a[i].x&&c<=a[i+1].x||c>=a[i+1].x&&c<=a[i].x)
{m0++;n0++;
pos1[m0].x=a[i].x;pos1[m0].y=a[i].y;pos2[n0].x=a[i+1].x;
pos2[n0].y=a[i+1].y;
}
}
//for(i=1;i<=2;i++)
//printf("%d %d %d %d\n",pos1[i].x,pos1[i].y,pos2[i].x,pos2[i].y);
if((pos2[1].x-pos1[1].x)!=0)
min=(pos2[1].y-pos1[1].y)*(c-pos1[1].x)/(pos2[1].x-pos1[1].x)+pos1[1].y;
else min=fun1(pos1[1].y,pos2[1].y);
if((pos2[2].x-pos1[2].x)!=0)
max=(pos2[2].y-pos1[2].y)*(c-pos1[2].x)/(pos2[2].x-pos1[2].x)+pos1[2].y;
else max=fun2(pos2[1].y,pos2[2].y);
printf("%.3lf %.3lf\n",min,max);
}
}//main
/*
3 1
0 2 1
0 0 1
5 2
-1 4 1 0 4
-1 3 1 3 0
*/
Followed by:
Post your reply here: |
All Rights Reserved 2003-2013 Ying Fuchen,Xu Pengcheng,Xie Di
Any problem, Please Contact Administrator