Submission #433581

#	Time	Username	Problem	Language	Result	Execution time	Memory
433581		nathanlee726	Naan (JOI19_naan)	C++14	5 / 100	4 ms	600 KiB

naan

//#include<i_am_noob_orz>
#include<bits/stdc++.h>
#include <ext/pb_ds/tree_policy.hpp>
#include <ext/pb_ds/assoc_container.hpp>
using namespace __gnu_pbds;
using namespace std;
#define ll long long
#define int ll
#define ull unsigned long long
#define pii pair<int,int>
#define X first
#define Y second
#define mod ((ll)1e9+7)
#define pb push_back
#define mp make_pair
#define abs(x) ((x)>0?(x):(-(x)))
#define F(n) Fi(i,n)
#define Fi(i,n) Fl(i,0,n)
#define Fl(i,l,n) for(int i=l;i<n;i++)
#define memres(a) memset(a,0,sizeof(a))
#define all(a) a.begin(),a.end()
#define sz(a) ((int)a.size())
#define ceiling(a,b) (((a)+(b)-1)/(b))
#define endl '\n'
#define bit_count(x) __builtin_popcountll((x))
#define ykh mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
#define jimmy_is_kind false
typedef tree<int, null_type, less<int>, rb_tree_tag, tree_order_statistics_node_update> rbtree;
 
//#define LOCAL
#ifdef LOCAL
#define bug(...) cerr<<"#"<<__LINE__<<' '<<#__VA_ARGS__<<"- ", _do(__VA_ARGS__)
template<typename T> void _do(T && x) {cerr<<x<<endl;}
template<typename T, typename ...S> void _do(T && x, S&&...y) {cerr<<x<<", "; _do(y...);}
#define IOS()
#else
#define IOS() ios_base::sync_with_stdio(0), cin.tie(0)
#define endl '\n'
#define bug(...)
#endif
 
int add(int a,int b){return(a+b>=mod?a+b-mod:a+b);}
int sub(int a,int b){return(a<b?a+mod-b:a-b);}
int po(int a,int b){
	if(b==0)return 1;
	if(b==1)return(a%mod);
	int tem=po(a,b/2);
	if(b&1)return(((tem*tem)%mod)*a)%mod;
	else return(tem*tem)%mod; 
}
int GCD(int a,int b){
	int x=0;
	int ra,rb;
	while(a&&b){
		if(((a&1)==0)&&((b&1)==0)){
			a>>=1,b>>=1,x++;
		}
		else if((a^b)&1){
			if(a&1)b>>=1;
			else a>>=1;
		}
		else{
			ra=abs(a-b),rb=min(a,b);
			a=ra,b=rb;
		}
	}
	return max(a,b)<<x;
}
int gcd(int a,int b){if(b==0)return a;return gcd(b,a%b);}
pii ADD(pii a,pii b){
	int c=a.Y*b.Y/GCD(a.Y,b.Y);
	a.X*=(c/a.Y),b.X*=(c/b.Y);
	return{a.X+b.X,c};
}
pii SUB(pii a,pii b){
	int c=a.Y*b.Y/GCD(a.Y,b.Y);
	a.X*=(c/a.Y),b.X*=(c/b.Y);
	return{a.X-b.X,c};
}
bool LAR(pii a,pii b){
	return a.X*b.Y>=a.Y*b.X;
}
int tab[2005][2005],s[2005];
pii gl[2005],ppo[2005],now[2005];
bool used[2005];
vector<int> v;
vector<pii> o;
signed main(){
	IOS();
	int n,l;
	cin>>n>>l;
	F(n)Fi(j,l)cin>>tab[i][j],s[i]+=tab[i][j];
	F(n){
		int g=GCD(s[i],n);
		gl[i]={s[i]/g,n/g};
	}
	F(n)ppo[i]={0,1},now[i]={0,1};
	pii lll={0,1},rr={0,1};
	F(n-1){
		pii mn={3000,1};
		int mip=-1;
		Fi(j,n){
			if(used[j])continue; 
			pii re=SUB(gl[j],now[j]);
			while(re.X>0){
				if(ppo[j].Y==1){
					if(LAR(re,{tab[j][ppo[j].X],1})){
						re=SUB(re,{tab[j][ppo[j].X],1});
						ppo[j].X++;
					}
					else{
						int c=GCD(re.X,tab[j][ppo[j].X]);
						ppo[j]=ADD(ppo[j],{re.X/c,tab[j][ppo[j].X]*re.Y/c});
						re={0,1};
						if(ppo[j].Y>1e8){
							int x=ppo[j].Y/1e8;
							ppo[j].Y/=x;
							ppo[j].X/=x;
							ppo[j].X++;
						}
					}
				}
				else{
					int np=ppo[j].X/ppo[j].Y;
					pii tmp=SUB({np+1,1},ppo[j]);
					int g=GCD(tab[j][np],tmp.Y);
					pii ttmp={tmp.X*tab[j][np]/g,tmp.Y/g};
					int c=GCD(re.X,tab[j][ppo[j].X]);
					if(LAR(ttmp,re)){
						ppo[j]=ADD(ppo[j],{re.X/c,re.Y*tab[j][np]/c});
						re={0,1};
						if(ppo[j].Y>1e8){
							int x=ppo[j].Y/1e8;
							ppo[j].Y/=x;
							ppo[j].X/=x;
							ppo[j].X++;
						}
					}
					else{
						ppo[j]={np+1,1};
						re=SUB(re,ttmp);
					}
				}
			}
			if(LAR(mn,ppo[j]))mn=ppo[j],mip=j;
			now[j]=gl[j];
		} 
		while(mip<0);
		used[mip]=1;
		v.pb(mip);
		o.pb(ppo[mip]);
		lll=rr;
		rr=ppo[mip];
		Fi(j,n){
			if(used[j])continue;
			pii ml=lll,mr=rr;
			while(ml!=mr){
				if((ml.X/ml.Y)==(mr.X/mr.Y)){
					pii tmp=(SUB(mr,ml));
					int c=GCD(tmp.Y,tab[j][(ml.X/ml.Y)]);
					pii xx={tmp.X*tab[j][(ml.X/ml.Y)]/c,tmp.Y/c};
					now[j]=SUB(now[j],xx);
					ml=mr;
					if(now[j].Y>1e8){
						int x=now[j].Y/1e8;
						now[j].Y/=x;
						now[j].X/=x;
						now[j].X++;						
					}
				}
				else{
					if(ml.Y==1){
						now[j]=SUB(now[j],{tab[j][ml.X],1});
						ml.X++;
					}
					else{
						pii tmp={ceiling(ml.X,ml.Y),1};
						tmp=SUB(tmp,ml);
						int c=GCD(tmp.Y,tab[j][(ml.X/ml.Y)]);
						now[j]=SUB(now[j],{tmp.X*tab[j][(ml.X/ml.Y)]/c,tmp.Y/c});
						ml.X=ceiling(ml.X,ml.Y),ml.Y=1;
						if(now[j].Y>1e8){
							int x=now[j].Y/1e8;
							now[j].Y/=x;
							now[j].X/=x;
							now[j].X++;						
						}
					}
				}
			}
		}
	}
	F(n)if(!used[i])v.pb(i);
	for(pii pi:o)cout<<pi.X<<" "<<pi.Y<<endl;
	for(int i:v)cout<<i+1<<" ";
	cout<<endl;
	return 0;
}

• Subtask 1: 5 / 5
• Subtask 2: 0 / 24
• Subtask 3: 0 / 71