hdu 4631 Sad Love Story

http://acm.hdu.edu.cn/showproblem.php?pid=4631

没想到这道题需要用“平均时间复杂度” 计算   一直没有想到解法  因为不符考虑了最坏情况的理念

方法一：

每加一个点 就找x值和它接近的 有可能更新最小距离的点进行判断更新 运行的相当的快 无语

方法二：

每求出所有点的最近点对 假如说是p[i],p[j](i<j) 那么在j之后加上的点不影响最小距离

递归j前面的所有点就可以了

由于自己的 求最近点对算法 写的相当烂 所有跑的相当慢，勉强过

代码1和代码2：

#include<iostream>

#include<cstdio>

#include<algorithm>

#include<string>

#include<cstring>

#include<cmath>

#include<set>

#include<vector>

using namespace std;



typedef long long ll;

typedef pair<double,double>ppd;

const double PI = acos(-1.);

const double eps = (1e-9);

const int N=1000010;

vector<int>vt[N];

int main()

{

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

    int T;

    scanf("%d",&T);

    while(T--)

    {

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

        vt[i].clear();

        int n;

        scanf("%d",&n);

        int A1,B1,C1,A2,B2,C2;

        cin>>A1>>B1>>C1;

        cin>>A2>>B2>>C2;

        int x=0,y=0;

        ll ans=0;

        ll tmp=(ll)(1e13);

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

        {

            x=((ll)A1*x+B1)%C1;

            y=((ll)A2*y+B2)%C2;

            for(int i=0;x+i<=1000000;++i)

            {

                ll x1=(ll)i*i;

                if(x1>=tmp)

                break;

                for(unsigned int j=0;j<vt[x+i].size();++j)

                tmp=min(tmp,x1+(ll)(y-vt[x+i][j])*(y-vt[x+i][j]));

            }

            for(int i=-1;x+i>=0;--i)

            {

                ll x1=(ll)i*i;

                if(x1>=tmp)

                break;

                for(unsigned int j=0;j<vt[x+i].size();++j)

                tmp=min(tmp,x1+(ll)(y-vt[x+i][j])*(y-vt[x+i][j]));

            }

            vt[x].push_back(y);

            if(tmp!=(ll)(1e13))

            ans+=tmp;

        }

        cout<<ans<<endl;

    }

    return 0;

}



#include<iostream>

#include<cstdio>

#include<algorithm>

#include<string>

#include<cstring>

#include<cmath>

#include<set>

#include<vector>

using namespace std;



typedef long long ll;

typedef pair<int,int>pp;

const double PI = acos(-1.);

const double eps = (1e-9);

const int N=1000010;

pp p[N];

struct node

{

    int x,y;

    int index;

}a[N],b[N];

ll dist(const node &a,const node &b)

{

    return (ll)(a.x-b.x)*(a.x-b.x)

    +(ll)(a.y-b.y)*(a.y-b.y);

}

bool cmpx(const node &a,const node &b)

{

    return a.x<b.x;

}

bool cmpy(const node &a,const node &b)

{

    return a.y<b.y;

}

ll dfs(int l,int r,int &k)

{

    if(r-l==1)

    {k=max(a[l].index,a[r].index);return dist(a[l],a[r]);}

    int m=(l+r)>>1;

    int k1,k2;

    ll d;

    ll d1=dfs(l,m,k1);

    ll d2=dfs(m,r,k2);

    if(d1<d2) {d=d1;k=k1;}

    else {d=d2;k=k2;}

    int ln=0,M=(int)(sqrt(1.0*d));

    for(int i=m+1;i<=r&&a[i].x-a[m].x<=M;++i)

    {b[ln]=a[i];++ln;}

    for(int i=m-1;i>=l&&a[m].x-a[i].x<=M;--i)

    {b[ln]=a[i];++ln;}

    sort(b,b+ln,cmpy);

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

    for(int j=1;j<=7&&i+j<ln;++j)

    if(dist(b[i],b[i+j])<d)

    {d=dist(b[i],b[i+j]);k=max(b[i].index,b[i+j].index);}

    return d;

}

ll solve(int n)

{

    ll ans=0;

    int r=n;

    while(r>=2)

    {

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

        {a[i].x=p[i].first;a[i].y=p[i].second;a[i].index=i;}

        sort(a+1,a+r+1,cmpx);

        int k;

        ll tmp=dfs(1,r,k);

        ans+=tmp*(r-k+1);

        r=k-1;

    }

    return ans;

}

int main()

{

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

    int T;

    scanf("%d",&T);

    while(T--)

    {

        int n;

        scanf("%d",&n);

        int A1,B1,C1,A2,B2,C2;

        scanf("%d %d %d",&A1,&B1,&C1);

        scanf("%d %d %d",&A2,&B2,&C2);

        p[0].first=0,p[0].second=0;

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

        {

            p[i].first=((ll)A1*p[i-1].first+B1)%C1;

            p[i].second=((ll)A2*p[i-1].second+B2)%C2;

        }

        cout<<solve(n)<<endl;

    }

    return 0;

}