oj.uz

Submission #717466

# Submission time Handle Problem Language Result Execution time Memory
717466 2023-04-02T00:30:11 Z Joshua_Andersson Portals (BOI14_portals) C++14
100 / 100
 457 ms 185308 KB 
#undef _GLIBCXX_DEBUG                // disable run-time bound checking, etc
#pragma GCC optimize("Ofast,inline") // Ofast = O3,fast-math,allow-store-data-races,no-protect-parens
#pragma GCC optimize ("unroll-loops")

#pragma GCC target("bmi,bmi2,lzcnt,popcnt")                      // bit manipulation
#pragma GCC target("movbe")                                      // byte swap
#pragma GCC target("aes,pclmul,rdrnd")                           // encryption
#pragma GCC target("avx,avx2,f16c,fma,sse3,ssse3,sse4.1,sse4.2") // SIMD

#include <bits/stdc++.h>
#include <bits/extc++.h>
using namespace std;

#define enablell 1

typedef long long ll;
typedef unsigned long long ull;
#if enablell
#define int ll
#define inf int(1e18)
#define float double
#else
const int inf = int(2e9);
#endif
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<vvi> vvvi;
typedef vector<vvvi> vvvvi;
typedef vector<bool> vb;
typedef vector<vb> vvb;
typedef vector<vvb> vvvb;
typedef pair<int, int> p2;
typedef vector<p2> vp2;
typedef vector<vp2> vvp2;
typedef vector<vvp2> vvvp2;
typedef tuple<int, int, int> p3;
typedef vector<p3> vp3;
typedef vector<vp3> vvp3;
typedef vector<vvp3> vvvp3;
typedef tuple<int, int, int, int> p4;
typedef vector<p4> vp4;

#define PBDS 0
#define _LOCAL _MSC_VER > 0
#if _LOCAL
#define gc() getchar()
#define popcount(x) __popcnt(x)
#define leading_zeros(x) _lzcnt_u32(x)
uint32_t clz(uint32_t x) { return _lzcnt_u32(x); }
uint32_t ctz(uint32_t x) { return _tzcnt_u32(x); }
#define bswap64(x) _byteswap_uint64(x)
#define assert(x) debassert(x)
#else
#define popcount(x) __builtin_popcount(x)
uint32_t clz(uint32_t x) { return __builtin_clz(x); }
uint32_t ctz(uint32_t x) { return __builtin_ctzll(x); }
#define bswap64(x) __builtin_bswap64(x)
#define gc() getchar_unlocked()
#if PBDS
using namespace __gnu_pbds;
// lower_bound is now upper_bound and vice versa (multiset). generally a bit broken
template<typename T> using indexed_multiset = tree<int, null_type, less_equal<T>, rb_tree_tag, tree_order_statistics_node_update>;
template<typename T> using indexed_set = tree<T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>;
struct chash { // large odd number for C
    const uint64_t C = ll(4e18 * acos(0)) | 71;
    ll operator()(ll x) const { return __builtin_bswap64(x * C); }
};

template<typename T, typename U> using fast_map = __gnu_pbds::gp_hash_table<T, U, chash>;
template<typename T> using fast_set = __gnu_pbds::gp_hash_table<T, null_type, chash>;
template<typename T, typename H> using fast_set_h = __gnu_pbds::gp_hash_table<T, null_type, H>;
#endif

#endif

#define FAST_INPUT 0
#define FILE_TC 0
#if FILE_TC && _LOCAL
ifstream filein("E:\\desktop\\code\\repos\\Comp prog\\x64\\Debug\\in.txt");
//ifstream filein("E:\\downloads\\test_data\\test_data\\005-case05.in");
//ifstream filein("E:\\desktop\\po-repos\\swedish-olympiad-2023\\online\\tomtplanering\\data\\secret\\group10\\010-case10.in");

#define cin filein
//ifstream cin("C:\\Users\\Matis\\desktop\\po-two\\swedish-olympiad-2014\\");
void fast() {}
#else
inline void fast() { ios::sync_with_stdio(false); cin.tie(NULL); cout.tie(NULL); }
#endif

#if FAST_INPUT && !FILE_TC
inline void read(int& v) { v = 0; int sign = 1; char c = gc(); if (c == '-') { sign = -1; } else { v += c - '0'; } while ((c = gc()) != EOF && c != ' ' && c != '\n') { v *= 10; v += c - '0'; } v *= sign; }
inline void read(int& u, int& v) { read(u); read(v); }
inline void read(int& u, int& v, int& k) { read(u); read(v); read(k); }
//inline void read(int& v) { char c; while ((c = getchar()) != EOF && c != ' ' && c != '\n') { v *= 10; v += c - '0'; } }
inline void read(string& s) { char c; while ((c = gc()) != EOF && c != '\n' && c != ' ') { s.push_back(c); } }
#else
template <typename T> inline void read(T& a) { cin >> a; }
template <typename T> inline void read(T& a, T& b) { cin >> a >> b; }
template <typename T> inline void read(T& a, T& b, T& c) { cin >> a >> b >> c; }
#endif
template <typename T> inline void write(T a) { cout << (a) << "\n"; }
#define quit cout << endl; _Exit(0);
#define dread(type, a) type a; read(a)
#define dread2(type, a, b) dread(type, a); dread(type, b)
#define dread3(type, a, b, c) dread2(type, a, b); dread(type, c)
#define dread4(type, a, b, c, d) dread3(type, a, b, c); dread(type, d)
#define dread5(type, a, b, c, d, e) dread4(type, a, b, c, d); dread(type, e)
#define readvector(type, name, size) vector<type> name(size); rep(i,size) {dread(type,temp); name[i]=temp;}
#ifdef _DEBUG
#define noop cout << "";
#define deb __debugbreak();
#define debassert(expr) if (!(expr)) deb;
#define debif(expr) if(expr) deb;
#else
#define noop ;
#define deb ;
#define debassert(expr) ;
#define debif(expr) ;
#endif

#define rep(i, high) for (int i = 0; i < high; i++)
#define repp(i, low, high) for (int i = low; i < high; i++)
#define repe(i, container) for (auto& i : container)
#define per(i, high) for (int i = high-1; i >= 0; i--)
#define perr(i, low, high) for (int i = high-1; i >= low; i--)

#define all(a) a.begin(),a.end()
#define rall(a) a.rbegin(),a.rend()
#define setcontains(set, x) (set.find(x) != set.end())
#define within(a, b, c, d) (a >= 0 && a < b && c >= 0 && c < d)
#define sz(container) ((int)container.size())
#define mp(a,b) (make_pair(a,b))

#define ceildiv(x,y) ((x + y - 1) / y)

template<typename T> inline void eraseOne(multiset<T>& mSet, T k) { auto itr = mSet.find(k); if (itr != mSet.end()) { mSet.erase(itr); } }
template<typename T, typename U> inline T first(U& a) { return *begin(a); }
template <typename T, typename U> inline void operator+=(pair<T, U>& l, const pair<T, U>& r) { l = { l.first + r.first,l.second + r.second }; }
template <typename T, typename U> inline pair<T, U> operator+(const pair<T, U> l, const pair<T, U> r) { return { l.first + r.first, l.second + r.second }; }
template <typename T, typename U> inline pair<T, U> operator-(const pair<T, U> l, const pair<T, U> r) { return { l.first - r.first, l.second - r.second }; }
template <typename T, typename U> inline pair<T, U> operator*(const pair<T, U> l, const int m) { return { l.first * m, l.second * m }; }
template <typename Out> inline void split(const string& s, char delim, Out result) { istringstream iss(s); string item; while (getline(iss, item, delim)) { *result++ = item; } }
inline vector<string> split(const string& s, char delim) { vector<string> elems; split(s, delim, back_inserter(elems)); return elems; }
vector<string> split(string s, string d) { size_t k = 0, n, e = d.length(); string t; vector<string> r; while ((n = s.find(d, k)) != string::npos) { t = s.substr(k, n - k); k = n + e; r.push_back(t); } r.push_back(s.substr(k)); return r; }
ll binpow(ll a, ll b) { ll res = 1; while (b > 0) { if (b & 1) res = res * a; a = a * a; b >>= 1; } return res; }
ll binpow(ll a, ll b, ll m) { a %= m; long long res = 1; while (b > 0) { if (b & 1) res = res * a % m; a = a * a % m; b >>= 1; } return res; } // For a < 2^31

#if 1
auto Start = chrono::high_resolution_clock::now();
void resettimer() { Start = chrono::high_resolution_clock::now(); }
int elapsedmillis() { return chrono::duration_cast<chrono::milliseconds>(chrono::high_resolution_clock::now() - Start).count(); }
random_device rd;
mt19937 rng(rd());
template<typename T> inline T randint(T lo, T hi) { return uniform_int_distribution<T>(lo, hi)(rng); } // [lo,hi]
template<typename T> inline T randel(vector<T>& v) { return v[uniform_int_distribution<T>((int)0, v.size() - 1)(rng)]; } // [lo,hi]
#endif
const ll mod = 1e9 + 7;
vp2 dirs = { {0,1},{0,-1},{1,0},{-1,0} };

vi newdirs = { 1,-1};
int getwall (int p, int dir, string& grid, vvi& dp)
{
    if (grid[p] == '#') return p - newdirs[dir];
    if (dp[p][dir] != -1) return dp[p][dir];
    return dp[p][dir] = getwall(p + newdirs[dir], dir, grid, dp);
}

int32_t main()
{
    fast();

    dread2(int, r, c);
    r += 2;
    c += 2;
    vector<string> grid;
    grid.reserve(r+2);
    grid.push_back(string(c, '#'));
    p2 start;
    p2 goal;

    rep(i, r-2)
    {
        dread(string, row);
        row.insert(row.begin(), '#');
        row.push_back('#');

        rep(j, c)
        {
            if (row[j]=='S')
            {
                row[j] = '.';
                start = { i+1,j };
            }
            else if (row[j]=='C')
            {
                row[j] = '.';
                goal = { i+1,j };
            }
        }
        grid.push_back(row);
    }
    grid.push_back(string(c, '#'));

    newdirs.push_back(c);
    newdirs.push_back(-c);

    string flatgrid;
    rep(i, r) flatgrid += grid[i];

    queue<p3> q;
    vi walldist(r*c, inf);
    rep(i, r)
    {
        rep(j, c)
        {
            repe(dir, dirs)
            {
                p2 p = p2(i,j) + dir;
                if (p.first < 0 || p.first >= r || p.second < 0 || p.second >= c) continue;
                if (grid[p.first][p.second] == '#' && grid[i][j]=='.')
                {
                    q.emplace(i, j, 0);
                }
            }
        }
    }
    
    while (q.size())
    {
        int i, j,d;
        tie(i, j,d) = q.front();
        q.pop();

        if (walldist[i*c+j] != inf) continue;
        walldist[i*c+j] = d;

        repe(dir, dirs)
        {
            p2 np = p2(i, j) + dir;
            if (np.first < 0 || np.first >= r || np.second < 0 || np.second >= c) continue;
            if (grid[np.first][np.second] == '#') continue;
            q.emplace(np.first, np.second, d + 1);
        }
    }

    vvi pq(r*c+10);
    pq[0].push_back(start.first*c+start.second);
    vi dist(r*c, inf);

    vvi dp(r*c, vi(4, -1));

    rep(d, pq.size()) while (pq[d].size())
    {
        int i = pq[d].back();
        pq[d].pop_back();

        if (i == goal.first*c+goal.second)
        {
            cout << d;
            quit;
        }

        if (dist[i] != inf) continue;
        dist[i] = d;

        repe(dir, newdirs)
        {
            int np = i + dir;

            if (flatgrid[np] == '#') continue;
            if (d+1<pq.size()) pq[d+1].emplace_back(np);
        }

        rep(dir, newdirs.size())
        {
            int np = getwall(i,dir,flatgrid,dp);

            if (d + 1 + walldist[i] < pq.size()) pq[d+walldist[i]+1].emplace_back(np);
        }
    }

    quit;
}

Compilation message

portals.cpp: In function 'int32_t main()':
portals.cpp:121:40: warning: comparison of integer expressions of different signedness: 'll' {aka 'long long int'} and 'std::vector<std::vector<long long int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  121 | #define rep(i, high) for (int i = 0; i < high; i++)
......
  252 |     rep(d, pq.size()) while (pq[d].size())
      |         ~~~~~~~~~~~~                    
portals.cpp:252:5: note: in expansion of macro 'rep'
  252 |     rep(d, pq.size()) while (pq[d].size())
      |     ^~~
portals.cpp:271:20: warning: comparison of integer expressions of different signedness: 'll' {aka 'long long int'} and 'std::vector<std::vector<long long int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  271 |             if (d+1<pq.size()) pq[d+1].emplace_back(np);
      |                 ~~~^~~~~~~~~~
portals.cpp:121:40: warning: comparison of integer expressions of different signedness: 'll' {aka 'long long int'} and 'std::vector<long long int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  121 | #define rep(i, high) for (int i = 0; i < high; i++)
......
  274 |         rep(dir, newdirs.size())
      |             ~~~~~~~~~~~~~~~~~~~         
portals.cpp:274:9: note: in expansion of macro 'rep'
  274 |         rep(dir, newdirs.size())
      |         ^~~
portals.cpp:278:37: warning: comparison of integer expressions of different signedness: 'll' {aka 'long long int'} and 'std::vector<std::vector<long long int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  278 |             if (d + 1 + walldist[i] < pq.size()) pq[d+walldist[i]+1].emplace_back(np);
      |                 ~~~~~~~~~~~~~~~~~~~~^~~~~~~~~~~

Subtask #1
11.0 / 11.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 388 KB Output is correct
2 Correct 0 ms 340 KB Output is correct
3 Correct 0 ms 340 KB Output is correct
4 Correct 1 ms 212 KB Output is correct
5 Correct 1 ms 340 KB Output is correct
6 Correct 1 ms 340 KB Output is correct
7 Correct 1 ms 340 KB Output is correct
8 Correct 0 ms 340 KB Output is correct
9 Correct 1 ms 212 KB Output is correct
10 Correct 1 ms 212 KB Output is correct

Subtask #2
20.0 / 20.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 384 KB Output is correct
2 Correct 1 ms 212 KB Output is correct
3 Correct 0 ms 340 KB Output is correct
4 Correct 1 ms 340 KB Output is correct
5 Correct 0 ms 340 KB Output is correct
6 Correct 0 ms 340 KB Output is correct
7 Correct 1 ms 340 KB Output is correct
8 Correct 1 ms 340 KB Output is correct
9 Correct 1 ms 852 KB Output is correct
10 Correct 1 ms 724 KB Output is correct
11 Correct 1 ms 724 KB Output is correct
12 Correct 1 ms 724 KB Output is correct
13 Correct 1 ms 724 KB Output is correct
14 Correct 1 ms 212 KB Output is correct
15 Correct 1 ms 724 KB Output is correct
16 Correct 0 ms 340 KB Output is correct

Subtask #3
20.0 / 20.0

# Verdict Execution time Memory Grader output
1 Correct 0 ms 212 KB Output is correct
2 Correct 1 ms 212 KB Output is correct
3 Correct 1 ms 340 KB Output is correct
4 Correct 0 ms 340 KB Output is correct
5 Correct 12 ms 6500 KB Output is correct
6 Correct 10 ms 6484 KB Output is correct
7 Correct 8 ms 5996 KB Output is correct
8 Correct 11 ms 6996 KB Output is correct

Subtask #4
19.0 / 19.0

# Verdict Execution time Memory Grader output
1 Correct 0 ms 340 KB Output is correct
2 Correct 1 ms 212 KB Output is correct
3 Correct 0 ms 340 KB Output is correct
4 Correct 1 ms 212 KB Output is correct
5 Correct 0 ms 340 KB Output is correct
6 Correct 1 ms 340 KB Output is correct
7 Correct 1 ms 340 KB Output is correct
8 Correct 0 ms 340 KB Output is correct
9 Correct 1 ms 852 KB Output is correct
10 Correct 1 ms 724 KB Output is correct
11 Correct 1 ms 724 KB Output is correct
12 Correct 1 ms 724 KB Output is correct
13 Correct 1 ms 724 KB Output is correct
14 Correct 10 ms 6560 KB Output is correct
15 Correct 10 ms 6520 KB Output is correct
16 Correct 8 ms 5912 KB Output is correct
17 Correct 9 ms 6484 KB Output is correct
18 Correct 10 ms 6484 KB Output is correct
19 Correct 5 ms 5076 KB Output is correct
20 Correct 9 ms 6868 KB Output is correct
21 Correct 11 ms 6564 KB Output is correct
22 Correct 10 ms 6368 KB Output is correct
23 Correct 9 ms 6336 KB Output is correct
24 Correct 11 ms 8064 KB Output is correct
25 Correct 0 ms 212 KB Output is correct
26 Correct 1 ms 724 KB Output is correct
27 Correct 1 ms 212 KB Output is correct
28 Correct 10 ms 6996 KB Output is correct

Subtask #5
30.0 / 30.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 212 KB Output is correct
2 Correct 0 ms 212 KB Output is correct
3 Correct 0 ms 340 KB Output is correct
4 Correct 1 ms 340 KB Output is correct
5 Correct 0 ms 340 KB Output is correct
6 Correct 0 ms 340 KB Output is correct
7 Correct 1 ms 340 KB Output is correct
8 Correct 0 ms 340 KB Output is correct
9 Correct 1 ms 724 KB Output is correct
10 Correct 1 ms 724 KB Output is correct
11 Correct 1 ms 724 KB Output is correct
12 Correct 1 ms 724 KB Output is correct
13 Correct 2 ms 724 KB Output is correct
14 Correct 10 ms 6484 KB Output is correct
15 Correct 11 ms 6468 KB Output is correct
16 Correct 8 ms 5972 KB Output is correct
17 Correct 9 ms 6484 KB Output is correct
18 Correct 10 ms 6484 KB Output is correct
19 Correct 5 ms 5076 KB Output is correct
20 Correct 8 ms 6868 KB Output is correct
21 Correct 13 ms 6532 KB Output is correct
22 Correct 9 ms 6272 KB Output is correct
23 Correct 9 ms 6344 KB Output is correct
24 Correct 239 ms 151308 KB Output is correct
25 Correct 457 ms 180464 KB Output is correct
26 Correct 140 ms 119580 KB Output is correct
27 Correct 265 ms 148816 KB Output is correct
28 Correct 270 ms 153080 KB Output is correct
29 Correct 267 ms 152636 KB Output is correct
30 Correct 236 ms 137136 KB Output is correct
31 Correct 12 ms 8148 KB Output is correct
32 Correct 379 ms 185308 KB Output is correct
33 Correct 1 ms 212 KB Output is correct
34 Correct 1 ms 724 KB Output is correct
35 Correct 225 ms 144532 KB Output is correct
36 Correct 0 ms 212 KB Output is correct
37 Correct 11 ms 6996 KB Output is correct
38 Correct 265 ms 165604 KB Output is correct
39 Correct 164 ms 135224 KB Output is correct