【jzoj4823】【小W学物理】【模拟】

题目大意

为了测试小W的物理水平，Mr.X在二维坐标系中放了N面镜子(镜子坐标绝对值不超过M)，镜子均与坐标轴成45°角，所以一共有两种类型“/”和“\”。原点不会有镜子，任意一点最多只有一面镜子。镜子两个面都能反光，而中间不透光，例如，对于一个“/”型镜子，下方向射入的光线会被反射到右方向，左方向射入的光线会被反射到上方向。现在有一条光线从原点沿X轴正方向射出，求走过T路程后所在位置。

解题思路

排序后求出每个镜子的上下左右，直接模拟。但是可能有环，需要判一下，由于数据太水，我就懒得判。

code

#include<set>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define LL long long
#define max(a,b) ((a>b)?a:b)
#define min(a,b) ((a>b)?b:a)
#define fo(i,j,k) for(int i=j;i<=k;i++)
#define fd(i,j,k) for(int i=j;i>=k;i--)
using namespace std;
int const inf=2147483647;
int const maxn=100000;
int n,m,next[maxn+10][10],w[4][2]={{0,1},{1,0},{0,-1},{-1,0}};LL t;
struct rec{
    int x,y,z,p;
};
rec a[maxn+10];
bool cmp1(rec x,rec y){return (x.x<y.x)||((x.x==y.x)&&(x.y<y.y));}
bool cmp2(rec x,rec y){return (x.y<y.y)||((x.y==y.y)&&(x.x<y.x));}
bool cmp3(rec x,rec y){return x.p<y.p;}
inline int dis(int x,int y){return abs(a[x].x-a[y].x)+abs(a[x].y-a[y].y);}
int main(){
    freopen("d.in","r",stdin);
    freopen("d.out","w",stdout);
    scanf("%d%d%lld",&n,&m,&t);
    fo(i,1,n){
        char ch;scanf("%d%d %c",&a[i].x,&a[i].y,&ch),a[i].p=i;
        a[i].z=(ch=='/');
    }
    sort(a+1,a+n+1,cmp1);
    fo(i,2,n)
        if(a[i].x==a[i-1].x){
            next[a[i].p][2]=a[i-1].p;
            next[a[i-1].p][0]=a[i].p;
        }
    sort(a+1,a+n+1,cmp2);
    int beg=0;
    if((!beg)&&(a[1].y==0))beg=a[1].p;
    fo(i,2,n){
        if(a[i].y==a[i-1].y){
            next[a[i].p][3]=a[i-1].p;
            next[a[i-1].p][1]=a[i].p;
        }
        if((!beg)&&(a[i].y==0)&&(a[i].x>0))beg=a[i].p;
    }
    sort(a+1,a+n+1,cmp3);
    int ans=0,fx=1;
    if(beg&&(t>=a[beg].y)){
        t-=a[beg].x;ans=beg;fx=3;
        for(;1;){
            if(a[ans].z)fx=3-fx;
            else fx=fx^1;
            if((next[ans][fx])&&(dis(ans,next[ans][fx])<=t)){
                t-=dis(ans,next[ans][fx]);
                ans=next[ans][fx];
                fx=fx^2;
            }else break;
        }
    }
    printf("%lld %lld",a[ans].x+1ll*w[fx][0]*t,a[ans].y+1ll*w[fx][1]*t);
    return 0;
}