Submission #887516

# Submission time Handle Problem Language Result Execution time Memory
887516 2023-12-14T16:56:18 Z nnhzzz trapezoid (balkan11_trapezoid) C++14
100 / 100
86 ms 17664 KB
/* [Author: Nguyen Ngoc Hung] - From THPT Ngo Gia Tu with Love */

#include <algorithm>
#include <bitset>
#include <complex>
#include <deque>
#include <exception>
#include <fstream>
#include <functional>
#include <iomanip>
#include <ios>
#include <iosfwd>
#include <iostream>
#include <istream>
#include <iterator>
#include <limits>
#include <locale>
#include <map>
#include <new>
#include <numeric>
#include <ostream>
#include <queue>
#include <set>
#include <unordered_set>
#include <sstream>
#include <stack>
#include <stdexcept>
#include <streambuf>
#include <string>
#include <utility>
#include <valarray>
#include <vector>
#include <cstring>
#include <unordered_map>
#include <cmath>
#include <array>
#include <cassert>
#include <random>
#include <chrono>

using namespace std;

#define        __nnhzzz__  signed main()
#define          BIT(i,j)  (((i)>>(j))&1LL)
#define           MASK(i)  (1LL<<(i))
#define            ALL(x)  (x).begin(),(x).end()
#define             SZ(x)  (int)(x).size()
#define                fi  first
#define                se  second
#define                ll  long long
#define                ld  long double
#define                vi  vector<int>
#define               vvi  vector<vi>
#define              vvvi  vector<vvi>
#define               pii  pair<int,int>
#define              vpii  vector<pii>
#define             vvpii  vector<vpii>
#define REPDIS(i,be,en,j)  for(int i = (be); i<=(en); i+=j)
#define     REPD(i,be,en)  for(int i = (be); i>=(en); i--)
#define      REP(i,be,en)  for(int i = (be); i<=(en); i++)
#define              endl  "\n"
#define                MP  make_pair
// #define               int  ll

//-----------------------------------------------------------------------------------------------//
int readInt(){
    char c;
    do{ c = getchar(); }while(c!='-' && !isdigit(c));
    bool neg = (c=='-');
    int res = neg?0:c-'0';
    while(isdigit(c=getchar())) res = (res<<3)+(res<<1)+(c-'0');
    return neg?-res:res;
}
//------------------------------------------------------------------------------------------------//
const  ll   LINF = 1e18;
const int    INF = 1e9;
const int    LOG = 20;
const int   MAXN = 5e5+7;
const int      N = 1e2+3;
const int    MOD = 30013;
const int   BASE = 1e5;
const  ld    EPS = 1e-9;
const  ld     PI = acos(-1);
const int OFFSET = 1e3;
//------------------------------------------------------------------------------------------------//
template<typename T1, typename T2> bool mini(T1 &a, T2 b){if(a>b){a=b;return true;}return false;}
template<typename T1, typename T2> bool maxi(T1 &a, T2 b){if(a<b){a=b;return true;}return false;}
template<typename T> T gcd(T a, T b) { while(b) { a %= b; swap(a,b); } return a; }
template<typename T> T lcm(T a, T b) { return a/gcd(a,b)*b; }
//------------------------------------------------------------------------------------------------//

/*
----------------------------------------------------------------
    END OF TEMPLATE
----------------------------------------------------------------
    Nguyen Ngoc Hung - nnhzzz
    Training for VOI24 gold medal
----------------------------------------------------------------
*/

struct Node{
	int x,y,u,v;

	bool operator < (const Node &other) const {
		return x<other.x;
	}
};

struct FENWICK_TREE{
    vi f;
    int n;

    void update(int id, int val){
        for(; id<=n; id+=id&-id){
            f[id] += val;
        }
    }

    int get(int id){
        int res = 0;
        for(; id>0; id-=id&-id){
            res += f[id];
        }
        return res;
    }

    int get(int l, int r){
        return get(r)-get(l-1);
    }

    FENWICK_TREE(int _):n(_){
        f.resize(n+1);
    }
};

Node a[MAXN];
pii dp[MAXN];
int n;

namespace sub1{
	void solve(){
		int res1 = 0,res2 = 0;
		REP(i,1,n){
			dp[i] = MP(1,1);
			int dem = 0;
			REP(j,1,i-1){
				if(a[i].x>a[j].y && a[i].u>a[j].v){
					++dem;
					if(dp[j].fi>=dp[i].fi) dp[i] = MP(dp[j].fi+1,dp[j].se);
					else if(dp[i].fi==dp[j].fi+1) dp[i].se = (dp[i].se+dp[j].se)%MOD;
				}
			}
			maxi(res1,dp[i].fi);
			cerr << dem;
		}
		REP(i,1,n) if(dp[i].fi==res1) res2 = (res2+dp[i].se)%MOD;
		cout << res1 << " " << res2;
	}
}

namespace sub2{
	struct events{
		int fi,se,op,id;

		bool operator < (const events &other) const {
			return fi<other.fi;
		}
	};
	int b[MAXN];
	pii f[MAXN];
	const int BASE = MAXN-7;

	pii combine(const pii &u, const pii &v){
		if(u.fi<v.fi) return v;
		if(v.fi<u.fi) return u;
		return MP(u.fi,(u.se+v.se)%MOD);
	}

	void update(int id, pii val){
		for(; id<=BASE; id+=id&-id) f[id] = combine(f[id],val);
	}

	pii get(int id){
		pii res = MP(0,1);
		for(; id>0; id-=id&-id) res = combine(res,f[id]);
		return res;
	}

	void solve(){
		vector<events> ev;
		vi cur;
		REP(i,1,n){
			ev.push_back({a[i].x,a[i].u,-1,i});
			ev.push_back({a[i].y,a[i].v,1,i});
			cur.push_back(a[i].u); cur.push_back(a[i].v);
		}
		sort(ALL(cur)); cur.resize(unique(ALL(cur))-cur.begin());
		sort(ALL(ev));
		pii res = MP(1,1);
		REP(i,0,SZ(ev)-1){
			int fi = ev[i].fi,se = ev[i].se,op = ev[i].op,id = ev[i].id;
			int pos = upper_bound(ALL(cur),se)-cur.begin();
			if(op==-1){
				dp[id] = get(pos);
				++dp[id].fi;
				res = combine(res,dp[id]);
			}else{
				update(pos,dp[id]);
			}
		}
		cout << res.fi << " " << res.se;
	}

	// void solve(){
	// 	pii res = MP(1,1);
	// 	vpii cur; vi comp;
	// 	REP(i,1,n){
	// 		cur.push_back(MP(a[i].x,a[i].u));
	// 		// cur.push_back(MP(a[i].y,a[i].v));
	// 		comp.push_back(a[i].v);
	// 		comp.push_back(a[i].u);
	// 	}
	// 	sort(ALL(cur)); cur.resize(unique(ALL(cur))-cur.begin());
	// 	sort(ALL(comp)); comp.resize(unique(ALL(comp))-comp.begin());
	// 	FENWICK_TREE fen(MAXN-1);
	// 	priority_queue<pii,vpii,greater<pii>> heap;
	// 	REP(i,1,n){
	// 		heap.push(MP(a[i].y,a[i].v));
	// 		while(SZ(heap)!=0 && heap.top().fi<a[i].x){
	// 			int pos = upper_bound(ALL(comp),heap.top().se)-comp.begin();
	// 			fen.update(pos,1);
	// 			heap.pop();
	// 		}
	// 		int pos = upper_bound(ALL(comp),a[i].u)-comp.begin();
	// 		b[i] = fen.get(pos);
	// 	}
	// 	REP(i,1,n){
	// 		int l = b[i];
	// 		int r = upper_bound(ALL(cur),MP(a[i].x,a[i].u))-cur.begin();
	// 		dp[i] = get(l-1); ++dp[i].fi;
	// 		update(r,dp[i]);
	// 		cout << dp[i].fi << " " << dp[i].se << endl;
	// 		maxi(res,dp[i]);
	// 	}
	// 	res = get(2*n);
	// 	cout << res.fi << " " << res.se;
	// }
}

void solve(){
	cin >> n;
	REP(i,1,n) cin >> a[i].x >> a[i].y >> a[i].u >> a[i].v;
	sort(a+1,a+n+1);
	// sub2::solve(); return ;
	if(n<=5000){
		sub1::solve();
		return ;
	}
	sub2::solve();
}

__nnhzzz__{
    ios::sync_with_stdio(0);
    cin.tie(0); cout.tie(0);
    #define task "test"
    if(fopen(task".inp","r")){
        freopen(task".inp","r",stdin);
        freopen(task".out","w",stdout);
    }
    #define task1 "nnhzzz"
    if(fopen(task1".inp","r")){
        freopen(task1".inp","r",stdin);
        freopen(task1".out","w",stdout);
    }

    int test = 1;

    while(test--){
      solve();
    }
    cerr << "\nTime elapsed: " << 1000*clock()/CLOCKS_PER_SEC << "ms\n";
    return 0;
}

/**  /\_/\
 *  (= ._.)
 *  / >TL \>AC
**/

Compilation message

trapezoid.cpp: In function 'void sub2::solve()':
trapezoid.cpp:48:28: warning: unused variable 'first' [-Wunused-variable]
   48 | #define                fi  first
      |                            ^~~~~
trapezoid.cpp:201:8: note: in expansion of macro 'fi'
  201 |    int fi = ev[i].fi,se = ev[i].se,op = ev[i].op,id = ev[i].id;
      |        ^~
trapezoid.cpp: In function 'int main()':
trapezoid.cpp:267:16: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
  267 |         freopen(task".inp","r",stdin);
      |         ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~
trapezoid.cpp:268:16: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
  268 |         freopen(task".out","w",stdout);
      |         ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~
trapezoid.cpp:272:16: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
  272 |         freopen(task1".inp","r",stdin);
      |         ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~
trapezoid.cpp:273:16: warning: ignoring return value of 'FILE* freopen(const char*, const char*, FILE*)' declared with attribute 'warn_unused_result' [-Wunused-result]
  273 |         freopen(task1".out","w",stdout);
      |         ~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~~
# Verdict Execution time Memory Grader output
1 Correct 1 ms 4444 KB Output is correct
2 Correct 1 ms 4444 KB Output is correct
3 Correct 2 ms 4444 KB Output is correct
4 Correct 4 ms 4444 KB Output is correct
5 Correct 8 ms 4616 KB Output is correct
6 Correct 15 ms 4444 KB Output is correct
7 Correct 34 ms 4676 KB Output is correct
8 Correct 25 ms 4676 KB Output is correct
9 Correct 9 ms 5284 KB Output is correct
10 Correct 15 ms 6472 KB Output is correct
11 Correct 20 ms 6788 KB Output is correct
12 Correct 44 ms 11256 KB Output is correct
13 Correct 49 ms 14076 KB Output is correct
14 Correct 69 ms 16828 KB Output is correct
15 Correct 64 ms 16116 KB Output is correct
16 Correct 71 ms 17412 KB Output is correct
17 Correct 72 ms 17152 KB Output is correct
18 Correct 72 ms 16892 KB Output is correct
19 Correct 76 ms 17264 KB Output is correct
20 Correct 86 ms 17664 KB Output is correct