【同一直线最多点】 poj 1118+2606+2780

poj 1118

#include<iostream>

using namespace std;

#define N 700

struct point {int x,y;} pnt[N];

int main()

{

    int n,i,j,k,max,cnt;

    while(scanf("%d",&n)&&n)

    {

        for(i=0;i<n;i++)

        scanf("%d%d",&pnt[i].x,&pnt[i].y);

        max=0;

        for(i=0;i<n;i++)

        for(j=i+1;j<n;j++)

        {

            cnt=0;

            for(k=j+1;k<n;k++)

            {

                if(k==i||k==j) continue;

                int left=(pnt[i].x-pnt[k].x)*(pnt[j].y-pnt[k].y);

                int right=(pnt[j].x-pnt[k].x)*(pnt[i].y-pnt[k].y);

                if(left==right) cnt++;

            }

            max=max>cnt? max:cnt;

        } 

        printf("%d\n",max+2);

    }

    return 0;

}

 

 

 

poj 2606

#include <iostream>

#include <cstdio>

#include <memory.h>

using namespace std;

const int maxn=202;

struct point

{

    int x;

    int y;

} pnt[maxn];

int main()

{

    //freopen("in.txt","r",stdin);

    int n,max,cnt;

    while(cin >> n)

    {

        memset(pnt,0,sizeof(pnt));

        for(int i=0;i<n;i++)

        {

            scanf("%d%d",&pnt[i].x,&pnt[i].y);

        }

        max=0;

        for(int i=0;i<n;i++)

        {

            for(int j=i+1;j<n;j++)

            {

                cnt=0;

                for(int k=j+1;k<n;k++)

                {

                    if(k==i || k==j)continue;

                    int left=(pnt[i].x-pnt[k].x)*(pnt[j].y-pnt[k].y);

                    int right=(pnt[j].x-pnt[k].x)*(pnt[i].y-pnt[k].y);

                    if(left==right)cnt++;

                }

                 if(cnt>max)max=cnt;

            }

        }

        cout << max+2 << endl;

    }

    return 0;

}

 

 

 

 

 

poj 2780

此题数据量比较大，同时虽然给的是3000ms，用n3算法也会超时，时间卡的紧，故一定要用n2lgn算法：

#include <iostream>

#include <cstdio>

#include <algorithm>

#include <cmath>

using namespace std;

const int maxn = 1002;

const double eps = 1e-6;

struct point

{

    int x,y;

}node[maxn];

double k[maxn];

int n;

bool equal(double x,double y)

{

    if(abs(x-y) <= eps)

    {

        return true;

    }

    return false;

}

int max(int a,int b)

{

    return a > b ? a : b;

}

void read()

{

    for(int i=0;i<n;i++)

    {

        scanf("%d %d",&node[i].x,&node[i].y);

    }

    return;

}

void solve()

{

    int ans = 0;

    for(int i=0;i<n-1;i++)

    {

        int top = 0;

        int tmp = 0;

        for(int j=i+1;j<n;j++)

        {

            if(node[i].x == node[j].x)

            {

                tmp++;

            }

            else

            {

                k[top++] = (double)(node[j].y - node[i].y) / (node[j].x - node[i].x);

            }

        }

        sort(k,k+top);

        int cnt = 1;

        for(int j=0;j<top;j++)

        {

            if(j < top-1 && equal(k[j],k[j+1]))

            {

                cnt++;

            }

            else

            {

                tmp = max(tmp , cnt);

                cnt = 1;

            }

        }

        ans = max(ans , tmp);

    }

    printf("%d\n",ans+1);

    return;

}

int main()

{

    while(~scanf("%d",&n))

    {

        read();

        solve();

    }

    return 0;

}

 

 

注意对于每个点，都要求出当前对于此点的无斜率情况和最大斜率的数目的max

之后枚举每个点的上述max

之前wa了无数次，承蒙九龙大神给了一个平行四边形的测试数据才豁然开朗，同时当时对与斜率不存在的情况也没有依点而分析，造成错误在所难免，下附原来的wrong answer 代码：

#include <iostream>

#include <cstdio>

#include <algorithm>

#include <cmath>

#include <memory.h>

using namespace std;

const int maxn=1002;

double s[500000+500];

struct location_

{

    int x;

    int y;

}l[maxn];

double slope(int m,int n)

{

    double ans=(double)(l[m].y-l[n].y)/(double)(l[m].x-l[n].x);

    return ans;

}

int main()

{

    //freopen("in.txt","r",stdin);

    int n;

    while(cin >> n)

    {

        memset(l,0,sizeof(l));

        memset(s,0.0,sizeof(s));

        int k=-1;

        int slo_cnt=0;

        for(int i=1;i<=n;i++)

        {

            scanf("%d%d",&l[i].x,&l[i].y);

        }

        for(int i=1;i<=n-1;i++)

        {

            for(int j=i+1;j<=n;j++)

            {

                if((l[i].x!=l[j].x))

                {

                    k++;

                    s[k]=slope(i,j);

                }

                else

                slo_cnt++;

            }

        }

        for(int i=2;i<=10000;i++)

        {

            if(i*(i-1)/2==slo_cnt)

            {

                slo_cnt=i;

                break;

            }

        }

//        cout << slo_cnt << ">>"<< endl;

//        for(int i=0;i<=k-1;i++)

//        {

//            cout << i << ' ' <<s[i] << endl;

//        }

        sort(s,s+k);

//        for(int i=0;i<=k-1;i++)

//        {

//            cout << i << ' ' <<s[i] << endl;

//        }

        int cnt=0;

        int max=0;

        for(int i=1;i<=k-1;i++)

        {

            if(fabs(s[i]-s[i-1])<=0.00000001)

            {

                cnt++;

                if(max<cnt)max=cnt;

            }

            else cnt=0;

        }

        max++;

//        cout << "dsdsdsds"<< "   " <<max << endl;

        for(int i=2;i<=10000;i++)

        {

            if(i*(i-1)/2==max)

            {

                max=i;

                break;

            }

        }

//        cout  << max<<"   max  " << endl;

//        cout << slo_cnt << " slo_cnt "<<endl;

        if(slo_cnt>max)max=slo_cnt;

        cout << max << endl;

    }

    return 0;

}