Submission #1069058

#TimeUsernameProblemLanguageResultExecution timeMemory
106905812345678Rectangles (IOI19_rect)C++17
10 / 100
421 ms386756 KiB
#include "rect.h"
#include <bits/stdc++.h>

using namespace std;

const int nx=2505;

int jr[nx][nx], jdn[nx][nx];
vector<pair<int, int>> dpr[nx][nx], dpdn[nx][nx];
stack<int> s;

struct fenwick
{
    int d[nx];
    void add(int i, int vl)
    {
        while (i<nx) d[i]+=vl, i+=(i&-i);
    }
    int query(int i)
    {
        int tmp=0;
        while (i>0) tmp+=d[i], i-=(i&-i);
        return tmp;
    }
} f;

long long count_rectangles(std::vector<std::vector<int> > a) {
	int n=a.size(), m=a[0].size();
	long long res=0;
    for (int i=0; i<n; i++) for (int j=0; j<m; j++) jdn[i][j]=jr[i][j]=-1;
    for (int i=n-1; i>=0; i--)
    {
        while (!s.empty()) s.pop();
        for (int j=m-1; j>=0; j--)
        {
            while (!s.empty()&&a[i][j]>=a[i][s.top()]) s.pop();
            if (!s.empty()) jr[i][j]=s.top();
            s.push(j);
            int curj=j+1;
            while (curj<m&&a[i][j]>a[i][curj])
            {
                curj=jr[i][curj];
                if (curj==-1) break;
                auto itr=lower_bound(dpr[i+1][j].begin(), dpr[i+1][j].end(), make_pair(curj, 0));
                if (itr!=dpr[i+1][j].end()&&itr->first==curj) dpr[i][j].push_back({curj, itr->second});
                else dpr[i][j].push_back({curj, i+1});
            }
        }
    }
    for (int j=m-1; j>=0; j--)
    {
        while (!s.empty()) s.pop();
        for (int i=n-1; i>=0; i--)
        {
            while (!s.empty()&&a[i][j]>=a[s.top()][j]) s.pop();
            if (!s.empty()) jdn[i][j]=s.top();
            s.push(i);
            int curi=i+1;
            while (curi<n&&a[i][j]>a[curi][j])
            {
                curi=jdn[curi][j];
                if (curi==-1) break;
                auto itr=lower_bound(dpdn[i][j+1].begin(), dpdn[i][j+1].end(), make_pair(curi, 0));
                if (itr!=dpdn[i][j+1].end()) dpdn[i][j].push_back({curi, itr->second});
                else dpdn[i][j].push_back({curi, j+1});
            }
        }
    }
    for (int i=0; i<n-1; i++)
    {
        for (int j=0; j<m-1; j++)
        {
            sort(dpr[i+1][j].begin(), dpr[i+1][j].end(), [](pair<int, int> a, pair<int, int> b){return a.second>b.second;});
            sort(dpdn[i][j+1].begin(), dpdn[i][j+1].end());
            reverse(dpdn[i][j+1].begin(), dpdn[i][j+1].end());
            int idx=0;
            for (auto x:dpdn[i][j+1])
            {
                while (idx<dpr[i+1][j].size()&&dpr[i+1][j][idx].second>=x.first) f.add(dpr[i+1][j][idx].first, 1), idx++;
                res+=f.query(x.second);
            }
            while (--idx>=0) f.add(dpr[i+1][j][idx].first, -1);
        }
    }
	return res;
}

/*
4 5
5 2 1 4 3
24 21 4 2 4
3 5 2 12 10
45 3 3 5 5
*/

Compilation message (stderr)

rect.cpp: In function 'long long int count_rectangles(std::vector<std::vector<int> >)':
rect.cpp:79:27: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<std::pair<int, int> >::size_type' {aka 'long unsigned int'} [-Wsign-compare]
   79 |                 while (idx<dpr[i+1][j].size()&&dpr[i+1][j][idx].second>=x.first) f.add(dpr[i+1][j][idx].first, 1), idx++;
      |                        ~~~^~~~~~~~~~~~~~~~~~~
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...