hdu 4361 2013多校联合训练第3场最后一题

题目有n个点  初始时先加入一个点到集合中,然后依次加入点,直到集合中的点的个数为n为止,每加一个点,算一次集合内所有点之间的最短距离,输出这个距离;为了避免输出过多,把每次要输出的结果相加后,每组只输出一个值;

细节见代码注释

#include<cstdio>
#include<algorithm>
#include<set>
#define INF 1e18
using namespace std;
long long n,Ax,Ay,Az,Bx,By,Bz;
struct point
{
	long long x,y;
	bool operator <(const point& b)const//递增排序
	{
		if(x==b.x) return y<b.y;//二级排序
		return x<b.x;
	}
};
typedef set<point>::iterator SetIt;
inline point &xy(point &x)//迭代计算,出下一个点
{
	x.x=(x.x*Ax+Ay)%Az;
	x.y=(x.y*Bx+By)%Bz;
	return x;
}
void deal()
{
	set<point> se;
	long long ans=0,milen=INF,temp,i;
	point now;now.x=now.y=0;
	se.insert(xy(now));//插入第一个点
	for(i=2;i<=n;i++)
	{
		SetIt it=se.lower_bound(xy(now));//二分查找,查找当前的点需要进位子
		SetIt it1=it,it2=it;
		while(it1!=se.begin())//向左扩张
		{
			--it1;
			temp=(it1->x-now.x)*(it1->x-now.x);
			if(temp>=milen) break;
			temp+=(it1->y-now.y)*(it1->y-now.y);
			milen=min(milen,temp);
		}
		while(it2!=se.end())//向右扩张
		{
			temp=(it2->x-now.x)*(it2->x-now.x);
			if(temp>=milen) break;
			temp+=(it2->y-now.y)*(it2->y-now.y);
			milen=min(milen,temp);
			++it2;
		}
		se.insert(now);// 插入
		ans+=milen;
	}
	printf("%I64d\n",ans);
}
int main()
{
	freopen("in.txt","r",stdin);
	int T;scanf("%d",&T	);
	while(T--)
	{
		scanf("%I64d%I64d%I64d%I64d%I64d%I64d%I64d",&n,&Ax,&Ay,&Az,&Bx,&By,&Bz);
		deal();
	}
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值