题面:
题意:
给一个等边三角形,初始时刻有一个点在其中,位置和速度矢量已知,这个球和三角形的边界发生完全弹性碰撞,且不记摩擦,求第k次碰撞发生的时间。k≤106
官方题解:
我们把空间看为一张这样的无限大的平面。
问题就变成了求小球做直线匀速运动,且与途中的边第k次相交是什么时候。
我们发现图中的边可以分为三组平行线。
我们把小球的速度和初始距离沿垂直于三组平行线的方向分解,那么可以在已知时间的情况下O(1)计算出这段时间内,小球与多少条边相交(分别计算三个方向分别相交了几次相加即可)。由于题目保证前k次以内不会碰到顶点,所以不会有重,之后发生重合也不影响二分的正确性。
上界估算一下就可以啦。
注意:eps太小会TLE。
代码:
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<string>
#include<queue>
#include<bitset>
#include<map>
#include<unordered_map>
#include<set>
#include<list>
#define ui unsigned int
#define ll long long
#define llu unsigned ll
#define ld long double
#define pr make_pair
#define pb push_back
#define lc (cnt<<1)
#define rc (cnt<<1|1)
//#define len(x) (t[(x)].r-t[(x)].l+1)
#define tmid ((l+r)>>1)
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)>(y)?(y):(x))
using namespace std;
const int inf=0x3f3f3f3f;
const ll lnf=0x3f3f3f3f3f3f3f3f;
const double dnf=1e18;
const int mod=1e9+7;
const double eps=1e-5;
const double pi=acos(-1.0);
const int hp=13331;
const int maxn=100100;
const int maxp=1100;
const int maxm=4000100;
const int up=1000;
double vh,vl,vr,dh,dl,dr;
double l,x,y,vx,vy,h;
int k;
int sgn(double x)
{
if(abs(x)<=eps) return 0;
else if(x<0) return -1;
else return 1;
}
ll check(double t)
{
ll ans=0;
if(sgn(vh)>0) ans+=(t*vh+dh)/h;
else if(sgn(vh)<0) ans+=(-t*vh+h-dh)/h;
if(sgn(vr)>0) ans+=(t*vr+dr)/h;
else if(sgn(vr)<0) ans+=(-t*vr+h-dr)/h;
if(sgn(vl)>0) ans+=(t*vl+dl)/h;
else if(sgn(vl)<0) ans+=(-t*vl+h-dl)/h;
return ans;
}
int main(void)
{
int tt;
scanf("%d",&tt);
while(tt--)
{
scanf("%lf%lf%lf%lf%lf%d",&l,&x,&y,&vx,&vy,&k);
//vh垂直于底边,vr垂直于右边的腰,vl垂直于左边的腰
vh=vy,vr=-vx*cos(pi/6)-vy*cos(pi/3),vl=vx*cos(pi/6)-vy*cos(pi/3);
h=l*sin(pi/3);
dh=y,dr=((h-y)*tan(pi/6)-x)*cos(pi/6),dl=((h-y)*tan(pi/6)+x)*cos(pi/6);
double nl=0,nr=l/sqrt(vx*vx+vy*vy)*k;
double mid;
while(nr-nl>=eps)
{
mid=(nl+nr)/2;
if(check(mid)>=k) nr=mid;
else nl=mid;
}
printf("%.8f\n",(nl+nr)/2);
}
return 0;
}