Submission #213265

#	Time	Username	Problem	Language	Result	Execution time	Memory
213265		dorijanlendvaj	Cultivation (JOI17_cultivation)	C++14	30 / 100	2100 ms	2816 KiB

cultivation

//DUEL
#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#define x first
#define y second
#define pii pair<int,int>
#define pb push_back
#define eb emplace_back
#pragma GCC optimize("unroll-loops")
#define shandom_ruffle(a, b) shuffle(a, b, rng)
#define vi vector<int>
#define vl vector<ll>
#define popcnt __builtin_popcount
#define popcntll __builtin_popcountll
#define all(a) begin(a),end(a)

using namespace std;
using namespace __gnu_pbds;

using ll=long long;
using ull=unsigned long long;
using ld=long double;
int MOD=1000000007;
int MOD2=998244353;
vector<int> bases;
const ll LLINF=1ll<<60;
const char en='\n';

mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
void yes() {cout<<"YES"<<en; exit(0);}
void no() {cout<<"NO"<<en; exit(0);}
inline int rund() {int x576363482791fuweh=rng();return abs(x576363482791fuweh)%RAND_MAX;}
template<class T>
void prVec(vector<T> w,bool fl=false)
{
	cout<<w.size()<<en;
	for (int i=0;i<int(w.size())-1;++i) cout<<w[i]<<' ';
	if (w.size()) cout<<w[w.size()-1]<<en;
	if (fl) cout<<flush;
}

void M998()
{
	swap(MOD,MOD2);
}

ll raand()
{
	ll a=rund();
	a*=RAND_MAX;
	a+=rund();
    return a;
}

#define rand raand

ll raaand()
{
	return raand()*(MOD-7)+raand();
}

string to_upper(string a)
{
	for (int i=0;i<(int)a.size();++i) if (a[i]>='a' && a[i]<='z') a[i]-='a'-'A';
	return a;
}

string to_lower(string a)
{
	for (int i=0;i<(int)a.size();++i) if (a[i]>='A' && a[i]<='Z') a[i]+='a'-'A';
	return a;
}

ll sti(string a,int base=10)
{
	ll k=0;
	for (int i=0;i<(int)a.size();++i)
	{
		k*=base;
		k+=a[i]-'0';
	}
	return k;
}

template<class T>
void eras(vector<T>& a,T b)
{
	a.erase(find(a.begin(),a.end(),b));
}

string its(ll k,int base=10)
{
	if (k==0) return "0";
	string a;
	while (k)
	{
		a.push_back((k%base)+'0');
		k/=base;
	}
	reverse(a.begin(),a.end());
	return a;
}

ll min(ll a,int b)
{
	if (a<b) return a;
	return b;
}

ll min(int a,ll b)
{
	if (a<b) return a;
	return b;
}

ll max(ll a,int b)
{
	if (a>b) return a;
	return b;
}

ll max(int a,ll b)
{
	if (a>b) return a;
	return b;
}

ll gcd(ll a,ll b)
{
	if (b==0) return a;
	return gcd(b,a%b);
}

ll lcm(ll a,ll b)
{
	return a/gcd(a,b)*b;
}

template<class T,class K>
pair<T,K> mp(T a,K b)
{
	return make_pair(a,b);
}

inline int mult(ll a,ll b)
{
	return (a*b)%MOD;
}

inline int pot(int n,int k)
{
	if (k==0) return 1;
	ll a=pot(n,k/2);
	a=mult(a,a);
	if (k%2) return mult(a,n);
	else return a;
}

int divide(int a,int b)
{
	return mult(a,pot(b,MOD-2));
}

inline int sub(int a,int b)
{
	if (a-b>=0) return a-b;
	return a-b+MOD;
}

inline int add(int a,int b)
{
	if (a+b>=MOD) return a+b-MOD;
	return a+b;
}

bool prime(ll a)
{
	if (a==1) return 0;
	for (int i=2;i<=round(sqrt(a));++i)
	{
		if (a%i==0) return 0;
	}
	return 1;
}

const int N=300010;
int n,m,k;
int x[N],y[N];
vi a,b;

int main()
{
	ios_base::sync_with_stdio(false);
	cin.tie(0);
	for (int i=0;i<10;++i) bases.push_back(rand()%(MOD-13893829*2)[phone]);
	cin>>n>>m>>k;
	for (int i=0;i<k;++i)
	{
		cin>>x[i]>>y[i];
		a.pb(x[i]);
		b.pb(y[i]);
	}
	for (int j=0;j<k;++j) for (int i=0;i<k-1;++i) if (y[i]>y[i+1]) swap(x[i],x[i+1]),swap(y[i],y[i+1]);
	a.pb(0);
	a.pb(n+1);
	b.pb(0);
	b.pb(m+1);
	sort(all(a));
	sort(all(b));
	int mu=a[1]-1,md=n-a[a.size()-2],mud=0;
	for (int i=1;i<a.size();++i) mud=max(mud,a[i]-a[i-1]-1);
	int ml=b[1]-1,mr=m-b[b.size()-2],mlr=0;
	for (int i=1;i<b.size();++i) mlr=max(mlr,b[i]-b[i-1]-1);
	set<int> d1,d2;
	for (int i=0;i<a.size();++i) for (int j=i;j<a.size();++j) d1.insert(max(0,a[j]-a[i]-1));
	for (int i=0;i<b.size();++i) for (int j=i;j<b.size();++j) d2.insert(max(0,b[j]-b[i]-1));
	vi r1(all(d1)),r2(all(d2));
	set<int> a1,b1;
	a1.insert(1);
	a1.insert(n);
	b1.insert(1);
	b1.insert(m);
	for (auto u: a) if (u!=0 && u!=n+1) a1.insert(u);
	for (auto u: b) if (u!=0 && u!=m+1) b1.insert(u);
	ll ans=MOD*4ll;
	for (auto u: r1) for (auto d: r1) if (u>=mu && d>=md && u+d>=mud && u+d<=ans)
	{
		vi mil,mir;
		int mlr1=mlr;
		vector<int> sy;
		for (int i=0;i<k;++i) if (x[i]-1<=u) sy.pb(y[i]);
		int la=0;
		for (int y1: sy)
		{
			mlr1=max(mlr1,y1-la-1);
			la=y1;
		}
		sy.clear();
		sy.reserve(k);
		for (int i=0;i<k;++i) if (n-x[i]<=d) sy.pb(y[i]);
		la=0;
		for (int y1: sy)
		{
			mlr1=max(mlr1,y1-la-1);
			la=y1;
		}
		sy.clear();
		sy.reserve(k);
		for (int y1: b1)
		{
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i) if (x[i]-1<=u)
			{
				if (y[i]<=y1) rr=min(rr,y1-y[i]);
				if (y[i]>=y1) rl=min(rl,y[i]-y1);
			}
			mil.pb(rl);
			mir.pb(rr);
		}
		for (int y1: b1)
		{
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i) if (n-x[i]<=d)
			{
				if (y[i]<=y1) rr=min(rr,y1-y[i]);
				if (y[i]>=y1) rl=min(rl,y[i]-y1);
			}
			mil.pb(rl);
			mir.pb(rr);
		}
		for (int x1: a1)
		{
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i)
			{
				bool ok=0;
				if (x[i]<=x1 && x1-x[i]<=d) ok=1;
				if (x[i]>=x1 && x[i]-x1<=u) ok=1;
				if (ok) rl=min(rl,y[i]-1);
			}
			mil.pb(rl);
			mir.pb(rr);
		}
		for (int x1: a1)
		{
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i)
			{
				bool ok=0;
				if (x[i]<=x1 && x1-x[i]<=d) ok=1;
				if (x[i]>=x1 && x[i]-x1<=u) ok=1;
				if (ok) rr=min(rr,m-y[i]),sy.pb(y[i]);
			}
			int la=0;
			for (int y1: sy)
			{
				mlr1=max(mlr1,y1-la-1);
				la=y1;
			}
			sy.clear();
			sy.reserve(k);
			mil.pb(rl);
			mir.pb(rr);
		}
		for (int x2: a1)
		{
			int x1=x2-u-1;
			if (x1<=1) continue;
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i)
			{
				bool ok=0;
				if (x[i]<=x1 && x1-x[i]<=d) ok=1;
				if (x[i]>=x1 && x[i]-x1<=u) ok=1;
				if (ok) rl=min(rl,y[i]-1);
			}
			mil.pb(rl);
			mir.pb(rr);
		}
		for (int x2: a1)
		{
			int x1=x2-u-1;
			if (x1<=1) continue;
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i)
			{
				bool ok=0;
				if (x[i]<=x1 && x1-x[i]<=d) ok=1;
				if (x[i]>=x1 && x[i]-x1<=u) ok=1;
				if (ok) rr=min(rr,m-y[i]),sy.pb(y[i]);
			}
			int la=0;
			for (int y1: sy)
			{
				mlr1=max(mlr1,y1-la-1);
				la=y1;
			}
			sy.clear();
			sy.reserve(k);
			mil.pb(rl);
			mir.pb(rr);
		}
		for (int x2: a1)
		{
			int x1=x2+d+1;
			if (x1>=n) continue;
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i)
			{
				bool ok=0;
				if (x[i]<=x1 && x1-x[i]<=d) ok=1;
				if (x[i]>=x1 && x[i]-x1<=u) ok=1;
				if (ok) rl=min(rl,y[i]-1),sy.pb(y[i]);
			}
			int la=0;
			for (int y1: sy)
			{
				mlr1=max(mlr1,y1-la-1);
				la=y1;
			}
			sy.clear();
			sy.reserve(k);
			mil.pb(rl);
			mir.pb(rr);
		}
		for (int x2: a1)
		{
			int x1=x2+d+1;
			if (x1>=n) continue;
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i)
			{
				bool ok=0;
				if (x[i]<=x1 && x1-x[i]<=d) ok=1;
				if (x[i]>=x1 && x[i]-x1<=u) ok=1;
				if (ok) rr=min(rr,m-y[i]);
			}
			mil.pb(rl);
			mir.pb(rr);
		}
		/*cout<<u<<' '<<d<<en;
		for (int i=0;i<mil.size();++i) cout<<mil[i]<<' '<<mir[i]<<en;*/
		for (ll l: mil) if (l>=ml)
		{
			ll r=0;
			for (int i=0;i<mir.size();++i) if (mil[i]>l) r=max(r,mir[i]);
			ans=min(ans,max(mlr1,l+max(r,mr))+u+d);
		}
		//cout<<ans<<en<<en<<en;
	}
	/*for (auto u: r1) for (auto ud: r1) if (u>=mu && ud-u>=md && ud>=mud && ud<=ans)
	{
		int d=ud-u;
		vi mil,mir;
		int mlr1=mlr;
		set<int> sy;
		for (int y1: b1)
		{
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i) if (x[i]-1<=u)
			{
				if (y[i]<=y1) rr=min(rr,y1-y[i]);
				if (y[i]>=y1) rl=min(rl,y[i]-y1);
			}
			mil.pb(rl);
			mir.pb(rr);
		}
		for (int y1: b1)
		{
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i) if (n-x[i]<=d)
			{
				if (y[i]<=y1) rr=min(rr,y1-y[i]);
				if (y[i]>=y1) rl=min(rl,y[i]-y1);
			}
			mil.pb(rl);
			mir.pb(rr);
		}
		for (int x1: a1)
		{
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i)
			{
				bool ok=0;
				if (x[i]<=x1 && x1-x[i]<=d) ok=1;
				if (x[i]>=x1 && x[i]-x1<=u) ok=1;
				if (ok) rl=min(rl,y[i]-1);
			}
			mil.pb(rl);
			mir.pb(rr);
		}
		for (int x1: a1)
		{
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i)
			{
				bool ok=0;
				if (x[i]<=x1 && x1-x[i]<=d) ok=1;
				if (x[i]>=x1 && x[i]-x1<=u) ok=1;
				if (ok) rr=min(rr,m-y[i]),sy.insert(y[i]);
			}
			int la=0;
			for (int y1: sy)
			{
				mlr1=max(mlr1,y1-la-1);
				la=y1;
			}
			sy.clear();
			mil.pb(rl);
			mir.pb(rr);
		}
		for (int x2: a1)
		{
			int x1=x2-u-1;
			if (x1<=1) continue;
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i)
			{
				bool ok=0;
				if (x[i]<=x1 && x1-x[i]<=d) ok=1;
				if (x[i]>=x1 && x[i]-x1<=u) ok=1;
				if (ok) rl=min(rl,y[i]-1);
			}
			mil.pb(rl);
			mir.pb(rr);
		}
		for (int x2: a1)
		{
			int x1=x2-u-1;
			if (x1<=1) continue;
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i)
			{
				bool ok=0;
				if (x[i]<=x1 && x1-x[i]<=d) ok=1;
				if (x[i]>=x1 && x[i]-x1<=u) ok=1;
				if (ok) rr=min(rr,m-y[i]),sy.insert(y[i]);
			}
			int la=0;
			for (int y1: sy)
			{
				mlr1=max(mlr1,y1-la-1);
				la=y1;
			}
			sy.clear();
			mil.pb(rl);
			mir.pb(rr);
		}
		for (int x2: a1)
		{
			int x1=x2+d+1;
			if (x1>=n) continue;
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i)
			{
				bool ok=0;
				if (x[i]<=x1 && x1-x[i]<=d) ok=1;
				if (x[i]>=x1 && x[i]-x1<=u) ok=1;
				if (ok) rl=min(rl,y[i]-1),sy.insert(y[i]);
			}
			int la=0;
			for (int y1: sy)
			{
				mlr1=max(mlr1,y1-la-1);
				la=y1;
			}
			sy.clear();
			mil.pb(rl);
			mir.pb(rr);
		}
		for (int x2: a1)
		{
			int x1=x2+d+1;
			if (x1>=n) continue;
			int rr=MOD,rl=MOD;
			for (int i=0;i<k;++i)
			{
				bool ok=0;
				if (x[i]<=x1 && x1-x[i]<=d) ok=1;
				if (x[i]>=x1 && x[i]-x1<=u) ok=1;
				if (ok) rr=min(rr,m-y[i]);
			}
			mil.pb(rl);
			mir.pb(rr);
		}
		//cout<<u<<' '<<d<<en;
		//for (int i=0;i<mil.size();++i) cout<<mil[i]<<' '<<mir[i]<<en;
		for (ll l: mil) if (l>=ml)
		{
			ll r=0;
			for (int i=0;i<mir.size();++i) if (mil[i]>l) r=max(r,mir[i]);
			ans=min(ans,max(mlr1,l+max(r,mr))+u+d);
		}
		//cout<<ans<<en<<en<<en;
	}*/
	cout<<ans<<en;
}

Compilation message (stderr)

cultivation.cpp: In function 'int main()':
cultivation.cpp:211:16: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
  for (int i=1;i<a.size();++i) mud=max(mud,a[i]-a[i-1]-1);
               ~^~~~~~~~~
cultivation.cpp:213:16: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
  for (int i=1;i<b.size();++i) mlr=max(mlr,b[i]-b[i-1]-1);
               ~^~~~~~~~~
cultivation.cpp:215:16: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
  for (int i=0;i<a.size();++i) for (int j=i;j<a.size();++j) d1.insert(max(0,a[j]-a[i]-1));
               ~^~~~~~~~~
cultivation.cpp:215:45: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
  for (int i=0;i<a.size();++i) for (int j=i;j<a.size();++j) d1.insert(max(0,a[j]-a[i]-1));
                                            ~^~~~~~~~~
cultivation.cpp:216:16: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
  for (int i=0;i<b.size();++i) for (int j=i;j<b.size();++j) d2.insert(max(0,b[j]-b[i]-1));
               ~^~~~~~~~~
cultivation.cpp:216:45: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
  for (int i=0;i<b.size();++i) for (int j=i;j<b.size();++j) d2.insert(max(0,b[j]-b[i]-1));
                                            ~^~~~~~~~~
cultivation.cpp:386:18: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
    for (int i=0;i<mir.size();++i) if (mil[i]>l) r=max(r,mir[i]);
                 ~^~~~~~~~~~~

• Subtask 1: 5 / 5
• Subtask 2: 10 / 10
• Subtask 3: 15 / 15
• Subtask 4: 0 / 30
• Subtask 5: 0 / 20
• Subtask 6: 0 / 20