[uva] 10136 Chocolate Chip Cookies

来源：互联网发布：视频爬虫 python 编辑：程序博客网时间：2024/05/19 12:23

【题目】http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&category=&problem=1077&mosmsg=Submission+received+with+ID+10618954

【题意】求一个饼干最多巧克力

【算法】简单几何+枚举

【题解】根据两个点恰在圆上枚举出可能的圆心，再遍历每个圆心，求其圆内有多少个巧克力。

【注意】注意输入和输出，都有空行- -

求圆心可以先求直线的法向量，然后标准化后，乘以长度，再加上线段中心就是圆心坐标了。

【测试】

4.0 4.0
4.0 5.0
1.0 20.0
1.0 21.0
1.0 22.0
1.0 25.0
1.0 26.0

1.0 1.0

1.0 0
0.0 1.0

ans:

【代码】

#include<iostream>#include<cstring>#include<cstdio>#include<cmath>using namespace std;struct point{    double x,y;}a[500],c[250000];char ch[100];double di(point x1,point x2){    return sqrt((x1.x-x2.x)*(x1.x-x2.x)+(x1.y-x2.y)*(x1.y-x2.y));}int main(){    int T,n,s,max=0,k,i,j;    double t,x,y,z;    double pi = 4*atan(1.0);    scanf("%d",&T);    getchar();    getchar();    while(T-->0){        if(max!=0)putchar('\n');        for(n=0;gets(ch)&&ch[0]!='\0';n++)sscanf(ch,"%lf%lf",&a[n].x,&a[n].y);        if(n==0) {            printf("0\n");            continue;        }        k=0;        for(i=1;i<n;i++)            for(j=0;j<i;j++){                t=di(a[i],a[j]);                if(t<=5){                    t = sqrt(6.25-t*t/4.0);                    //printf("%.2f\n",t);                    x=a[j].y-a[i].y;                    y=a[i].x-a[j].x;                    z=sqrt(x*x+y*y);                    x=x/z*t;           //标准法向量                    y=y/z*t;                    c[k].x=(a[i].x+a[j].x)/2.0+x;                    c[k++].y=(a[i].y+a[j].y)/2.0+y;                    c[k].x=(a[i].x+a[j].x)/2.0-x;                    c[k++].y=(a[i].y+a[j].y)/2.0-y;                    //printf("%d %d %.2f %.2f %.2f %.2f\n",i,j,c[k-2].x,c[k-2].y,c[k-1].x,c[k-1].y);                }            }        max=1;        for(i=0;i<k;i++){            s=0;            for(j=0;j<n;j++){                //printf("%f\n",di(c[i],a[j]));                if(di(c[i],a[j])<=2.5)s++;            }            if(s>max)max=s;        }        printf("%d\n",max);    }    return 0;}