HDU4631(标程代码）

 /*将x从小到大排序，每次插入一个点，直接找比这个点的x大的第一个，然后从这个开始向两边找
 ，找点的下标用多重容器实现*/
 #include<stdio.h>
 #include<string.h>
 #include<algorithm>
 #include<set>
 #include<algorithm>
 #include<iostream>
 using namespace std;
 struct point
 {
     __int64 x,y;
     bool operator<(const point &a) const
     {
         return x<a.x;
     }
 };
 void cacl(__int64 n,__int64 x[])
 {
     __int64 a,b,c;
     scanf("%I64d%I64d%I64d",&a,&b,&c);
     int i;
     x[]=;
     for(i=;i<=n;i++)
      x[i]=(x[i-]*a+b)%c;
 }
 __int64 min(__int64 a,__int64 b)
 {
     if(a<b) return a;
     return b;
 }
 __int64 x[],y[],n;
 void find()
 {
    __int64 i,j;
    __int64 ans,sum;
    __int64 m1,m2;
    scanf("%I64d",&n);
     cacl(n,x);
     cacl(n,y);
       //for(i=1;i<=n;i++)
     //    printf("%I64d %I64d\n",x[i],y[i]);
    ans=(__int64)<<;
    sum=;
    multiset<point> M;
    for(i=;i<=n;i++)
    {
         point p;
        p.x=x[i]; p.y=y[i];

        if(i>)//从第二个点开始
        {
          multiset<point>::iterator it=M.lower_bound(p),e;
          for(e=it;e!=M.end();e++)
          {
            m1=e->x-p.x;////差值
            if(m1*m1>=ans) break;
            m2=e->y-p.y;
            ans=min(ans,m1*m1+m2*m2);
          }
          for(e=it;e!=M.begin();)
          {
              e--;//防止重复，上面已经找过it了
            m1=e->x-p.x;
            if(m1*m1>=ans) break;
            m2=e->y-p.y;
            ans=ans=min(ans,m1*m1+m2*m2);
          }
          sum=sum+ans;
        }
       M.insert(p);
    }
    printf("%I64d\n",sum);
 }
 int main()
 {
     int t;
     scanf("%d",&t);
     while(t--)
     {
         find();
     }
     return ;
 }