Submission #955442

#	Time	Username	Problem	Language	Result	Execution time	Memory
955442		vjudge1	IOI Fever (JOI21_fever)	C++17	6 / 100	99 ms	208980 KiB

fever

#include <bits/stdc++.h>

#pragma optimize("Ofast")
#pragma target("avx2")

using namespace std;

#define ll long long
#define ld long double
#define pb push_back
#define pf push_front
#define pii pair<int, int>
#define all(v) v.begin(), v.end()
#define F first
#define S second
#define mem(a, i) memset(a, i, sizeof(a))
#define sz(s) (int)s.size()
#define y1 yy
#define ppb pop_back
#define lb lower_bound
#define ub upper_bound
#define gcd(a, b) __gcd(a, b)
#define in insert
#define int ll

const int MAX = 1e6 + 15;
const int B = 300;
const int N = 104;
const int block = 400;
const int maxB = MAX / B + 10;
const ll inf = 1e18;
const int mod = 998244353;
const int mod1 = 1e9 + 9;
const ld eps = 1e-9;

// int dx[8]={1,0,-1,0,1,-1,-1,1};
// int dy[8]={0,1,0,-1,1,-1,1,-1};

int binpow(int a, int n)
{
  if (!n)
    return 1;
  if (n % 2 == 1)
    return a * binpow(a, n - 1) % mod;
  int k = binpow(a, n / 2);
  return k * k % mod;
}

int n;
pii a[MAX];
int b[MAX];
int d[MAX];
int dx[MAX], dy[MAX];
set<int> vx[MAX], vy[MAX];
int ans;

void calc(int start, int nap)
{
  set<pii> st;
  for (int i = 1; i <= n; i++)
    d[i] = inf, b[i] = 0;
  b[start] = nap;
  d[start] = 0;
  st.in({0, start});
  while (!st.empty())
  {
    int v = st.begin()->S;
    st.erase(st.begin());
    vector<int> del;
    vx[dx[v]].erase(v);
    vy[dy[v]].erase(v);
    for (auto to : vx[dx[v]])
    {
      if (d[to] > abs(a[v].F - a[to].F) && d[v] < abs(a[v].F - a[to].F))
      {
        if (b[v] == 1)
        {
          if (a[to].F < a[v].F)
          {
            assert(d[to] == inf);
            b[to] = 2;
            st.erase({d[to], to});
            d[to] = abs(a[v].F - a[to].F);
            st.in({d[to], to});
            // del.pb(to);
          }
        }
        if (b[v] == 2)
        {
          if (a[to].F > a[v].F)
          {
            assert(d[to] == inf);
            b[to] = 1;
            st.erase({d[to], to});
            d[to] = abs(a[v].F - a[to].F);
            st.in({d[to], to});
            // del.pb(to);
          }
        }
        if (b[v] == 3)
        {
          if (a[to].F > a[v].F)
          {
            assert(d[to] == inf);
            b[to] = 4;
            st.erase({d[to], to});
            d[to] = abs(a[v].F - a[to].F);
            st.in({d[to], to});
            // del.pb(to);
          }
        }
        if (b[v] == 4)
        {
          if (a[to].F < a[v].F)
          {
            b[to] = 3;
            assert(d[to] == inf);
            st.erase({d[to], to});
            d[to] = abs(a[v].F - a[to].F);
            st.in({d[to], to});
            // del.pb(to);
          }
        }
      }
    }
    for (int x : del)
      vx[dx[v]].erase(x);
    del.clear();
    for (auto to : vy[dy[v]])
    {
      if (d[to] > abs(a[v].F - a[to].F) && d[v] < abs(a[v].F - a[to].F))
      {
        if (b[v] == 1)
        {
          if (a[to].F > a[v].F)
          {
            assert(d[to] == inf);
            b[to] = 4;
            st.erase({d[to], to});
            d[to] = abs(a[v].F - a[to].F);
            st.in({d[to], to});
            // del.pb(to);
          }
        }
        if (b[v] == 2)
        {
          if (a[to].F > a[v].F)
          {
            assert(d[to] == inf);
            b[to] = 3;
            st.erase({d[to], to});
            d[to] = abs(a[v].F - a[to].F);
            st.in({d[to], to});
            // del.pb(to);
          }
        }
        if (b[v] == 3)
        {
          if (a[to].F < a[v].F)
          {
            assert(d[to] == inf);
            b[to] = 2;
            st.erase({d[to], to});
            d[to] = abs(a[v].F - a[to].F);
            st.in({d[to], to});
            // del.pb(to);
          }
        }
        if (b[v] == 4)
        {
          if (a[to].F < a[v].F)
          {
            assert(d[to] == inf);
            b[to] = 1;
            st.erase({d[to], to});
            d[to] = abs(a[v].F - a[to].F);
            st.in({d[to], to});
            // del.pb(to);
          }
        }
      }
    }
    for (int x : del)
      vy[dy[v]].erase(x);
    del.clear();
  }
  int cur = 0;
  for (int i = 1; i <= n; i++)
  {
    cur += (d[i] != inf);
    // cout<<d[i]<<" "<<b[i]<<"\n";
  }
  ans = max(ans, cur);
}

void solve()
{
  cin >> n;
  for (int i = 1; i <= n; i++)
    cin >> a[i].F >> a[i].S;
  set<int> d;
  map<int, int> mp;
  int cur = 0;

  for (int k = 1; k <= 4; k++)
  {
    for (int i = 1; i <= n; i++)
    {
      d.in(a[i].F + a[i].S);
    }
    for (int x : d)
      mp[x] = ++cur;
    for (int i = 1; i <= n; i++)
    {
      dx[i] = mp[a[i].F + a[i].S];
      vx[mp[a[i].F + a[i].S]].in(i);
    }

    cur = 0;
    mp.clear();
    d.clear();

    for (int i = 1; i <= n; i++)
    {
      d.in(a[i].F - a[i].S);
    }
    for (int x : d)
      mp[x] = ++cur;
    for (int i = 1; i <= n; i++)
    {
      dy[i] = mp[a[i].F - a[i].S];
      vy[mp[a[i].F - a[i].S]].in(i);
    }
    cur = 0;
    mp.clear();
    d.clear();
    calc(1, k);
  }
  cout << ans;
}

signed main()
{
  //  freopen("help.in","r",stdin);
  //  freopen("help.out","w",stdout);
  ios_base::sync_with_stdio(0);
  cin.tie(0);
  cout.tie(0);
  // prec();
  int t = 1;
  // cin>>t;
  while (t--)
    solve();
}

Compilation message (stderr)

fever.cpp:3: warning: ignoring '#pragma optimize ' [-Wunknown-pragmas]
    3 | #pragma optimize("Ofast")
      | 
fever.cpp:4: warning: ignoring '#pragma target ' [-Wunknown-pragmas]
    4 | #pragma target("avx2")
      |

• Subtask 1: 0 / 5
• Subtask 2: 0 / 8
• Subtask 3: 6 / 6
• Subtask 4: 0 / 6
• Subtask 5: 0 / 12
• Subtask 6: 0 / 32
• Subtask 7: 0 / 31