HDOJ 6798 Triangle Collision(杭电多校2020第三场1008)(二分)题解

http://acm.hdu.edu.cn/showproblem.php?pid=6798
思路:
画出等边三角形相接的网格。
在这里插入图片描述

这条直线和网格的交点即在三角形中的每一次碰撞。从(x,y)出发的朝(vx,vy)方向的射线与网格交k次之后得到的长度即碰撞k次经过的总长。
转化为射线与网格相交的问题。

二分查找需要的结果t。每次将整个三角形、点、速度方向一起绕三角形中心旋转0度,120度,240度,每次旋转的时候只计算和底边(即网格中与x轴平行的边)加起来就是射线和网格的所有交点数量。看大于k还是小于k。
代码:

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<map>
#include<queue>
#include<cmath>
#include<functional>
using namespace std;

template<class T>inline void read(T &x){x=0;char o,f=1;while(o=getchar(),o<48)if(o==45)f=-f;do x=(x<<3)+(x<<1)+(o^48);while(o=getchar(),o>47);x*=f;}
//int cansel_sync=(ios::sync_with_stdio(0),cin.tie(0),0);
#define ll long long
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define repb(i,a,b) for(int i=(a);i>=b;i--)
#define INF 0x3f3f3f3f
#define cendl printf("\n")
ll gcd(ll a,ll b){ while(b^=a^=b^=a%=b); return a; }
//#define INF 0x7fffffff

const double eps = 1e-8;
double x,y,vx,vy,L;
double h,cardv;
int k;
const double Pi = 2*acos(0);

int get_jd(double x,double y,double vx,double vy,double t){//计算有几个交点 
	if(vy==0){return 0;}
	if(vy>0){ vy=-vy; y=h-y; }
	int ret = 1;
	double tl = (-vy)*t-y;
	ret += floor(tl/h);//每次只看和x轴平行的角线 
	return ret;
}

void solve(){
	scanf("%lf%lf%lf%lf%lf%d",&L,&x,&y,&vx,&vy,&k);
	h = L*sqrt(3)/2;
	int tot = 0;
	double zx = 0;//中心点坐标 
	double zy = h/3;
	double xx,yy,vvx,vvy;
	double mid;
	cardv = sqrt(vx*vx+vy*vy);
	double l=0,r=k*L/cardv;
	while(r-l>eps){
		mid = (l+r)/2;
		tot= get_jd(x,y,vx,vy,mid);//统计交点总数 
		//旋转 120度 
		xx = (x-zx)*cos(2*Pi/3) - (y-zy)*sin(2*Pi/3) + zx;
		yy = (x-zx)*sin(2*Pi/3) + (y-zy)*cos(2*Pi/3) + zy;
		vvx = vx*cos(2*Pi/3) - vy*sin(2*Pi/3);
		vvy = vx*sin(2*Pi/3) + vy*cos(2*Pi/3);
		tot+=get_jd(xx,yy,vvx,vvy,mid);
		//再旋转 240度 
		xx = (x-zx)*cos(4*Pi/3) - (y-zy)*sin(4*Pi/3) + zx;
		yy = (x-zx)*sin(4*Pi/3) + (y-zy)*cos(4*Pi/3) + zy;
		vvx = vx*cos(4*Pi/3) - vy*sin(4*Pi/3);
		vvy = vx*sin(4*Pi/3) + vy*cos(4*Pi/3);
		tot+=get_jd(xx,yy,vvx,vvy,mid);
		if(tot>=k) r = mid;
		else l = mid;
	}
	printf("%.8llf\n",mid);
}

int main(){
	int z;
	scanf("%d",&z);
	while(z--) solve();
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值