CSP模拟 题目C - 可怕的宇宙射线

题意

众所周知，瑞神已经达到了CS本科生的天花板，但殊不知天外有天，人外有苟。在浩瀚的宇宙中，存在着-种叫做苟狗的生物， 这种生物天生就能达到人类研究生的知识水平，并且天生擅长CSP,甚至有全国第一的水平!但最可怕的是，它可以发出宇宙射线!宇宙射线可以摧毁人的智商，进行降智打击!

宇宙射线会在无限的二维平面上传播(可以看做一个二维网格图)，初始方向默认向上。宇宙射线会在发射出一段距离后分裂,向该方向的左右45°方向分裂出两条宇宙射线，同时威力不变!宇宙射线会分裂n次,每次分裂后会在分裂方向前进a_i个单位长度。

现在瑞神要带着他的小弟们挑战苟狗，但是瑞神不想让自己的智商降到普通本科生zjm那么菜的水平,所以瑞神来请求你帮他计算出共有多少个位置会被"降智打击”。

---

思路

运用偏移数组dx与dy，表示点分裂时可以延伸搜索的方向

//偏移数组 上 右上 右 下 左下 左 左上
int dx[8] = {0,1,1,1,0,-1,-1,-1};
int dy[8] = {1,1,0,-1,-1,-1,0,1};

朴实的方法就是一般的dfs搜索，规定搜索边界，设置好每个方向的搜索之后让dfs把到达过的点全部都插入到一个集合中

set<point> s;


    for (int i = 0; i < v[k]; i++)
    {//向前前进v[k]
        _x = _x + ddx;
        _y = _y + ddy;
        point temp(_x,_y);
        s.insert(temp);
    }

set有自动去重的功能，将射线经过的坐标插入set中，遍历完成时set中的元素个数就是射线经过的点的个数

现在考虑剪枝的策略，如何剪枝，如何及时终止无用的dfs搜索？
射线在不同时间段前进a_i次后就会分裂，考虑到一个射线i要分裂的时候，进行第k次分裂，分裂方向为ddx₁ + ddy₁,ddx₂ + ddy₂,
如果这个分裂点p已经有相同分裂次数k且分裂方向相同ddx₁ + ddy₁,ddx₂ + ddy₂,的射线i(i != j)经过，那么射线i往后的所有路线j都已经走过一次了，可以不用再次访问，实现一个这样的剪枝策略

为了剪枝就要有标记
用一个四维数组标记，一个分裂点的坐标，分裂点，是第几次分裂，分裂的方向是什么

bool mark[300][300][8][31] ={false};
//记录一个点在一个坐标上分裂时，如何分裂，第几次分裂

---

总结

这是补题，模拟的时候没有写，现在无时间限制的话，就可以考虑用偏移数组来简化代码，以后多试试，改过一个错误：之前的四维数组标记中，标记的不是第k次分裂而是分裂的长度l，以为分裂后前进长度一样的射线就可以剪枝了，但是这时错的，因为每次分裂后前进的长度是不确定的，可能相等也可能不相等，如果遇到了不同的分裂次数k_i k_j，但是他们的前进长度确一样的话V[k_i] == V[k_j],会错误的把一次搜索回溯

dfs经过剪枝策略之后，就可以提高dfs的效率，防止dfs进行多余的搜索

---

代码

#include <stdio.h>
#include <cstdio>
#include <algorithm>
#include <set>
using namespace std;
//偏移数组 上 右上 右 下 左下 左 左上
int dx[8] = {0,1,1,1,0,-1,-1,-1};
int dy[8] = {1,1,0,-1,-1,-1,0,1};

int n;
int v[31];//射线行走长度
bool mark[300][300][8][31] ={false};
//记录一个点在一个坐标上分裂时，如何分裂，第几次分裂

struct point
{
    int x;
    int y;
    point(int _x,int _y)
    {
        x = _x;
        y = _y;
    }
    inline bool friend operator<(const point & a,const point & b)
    {
        if(a.x != b.x)  return a.x < b.x;
        else return a.y < b.y;
    }
};

set<point> s;
//一个射线起点a，下一次分裂方向为ddx ddy，分裂方式为way，行走长度为v[k]
void dfs(point a,int ddx,int ddy,int way,int k)
{
    if (k == n)//n-1是最后一轮，递归边界
    {
        return;
    }
    int _x = a.x;
    int _y = a.y;
    if (mark[_x][_y][way][k])
    {
        return;
    }
    mark[_x][_y][way][k] = true;
    for (int i = 0; i < v[k]; i++)
    {//向前前进v[k]
        _x = _x + ddx;
        _y = _y + ddy;
        point temp(_x,_y);
        s.insert(temp);
    }
    //分裂到其他两个方向：
    point b(_x,_y);
    dfs(b,dx[(way + 1) % 8],dy[(way + 1) % 8],((way + 1) % 8),k+1);
    dfs(b,dx[(way + 7) % 8],dy[(way + 7) % 8],((way + 7) % 8),k+1);
}

int main()
{
    scanf("%d",&n);
    for (int i = 0; i < n; i++)
    {
        scanf("%d",&v[i]);
    }
    point start(150,150);
    dfs(start,0,1,0,0);
    printf("%d\n",s.size());
    return 0;
}