제출 #1347924

#제출 시각아이디문제언어결과실행 시간메모리
1347924PieArmyNaan (JOI19_naan)C++20
29 / 100
49 ms4716 KiB
#include<bits/stdc++.h>
typedef long long ll;
#define pb push_back
#define fr first
#define sc second
#define endl '\n'
using namespace std;
#define mid ((left+right)>>1)
#define int ll
int n,m;
int sum[2023];
int arr[2023][2023];
pair<int,int>ans[2023];
vector<int>v;

int gcd(int x,int y){
	if(x<y)swap(x,y);
	if(!y)return x;
	return gcd(y,x%y);
}
void print(__int128_t x){
	if(x<0){
		cerr<<"-";
		x=-x;
	}
	if(x>9)print(x/10);
	cerr<<int(x%10);
}

signed main(){
	ios_base::sync_with_stdio(23^23);cin.tie(0);
	cin>>n>>m;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			cin>>arr[i][j];
			sum[i]+=arr[i][j];
			arr[i][j]+=arr[i][j-1];
		}
	}
	pair<int,int>pos={0,1};
	for(int j=1;j<=n;j++){
		int opt=0;
		pair<__int128_t,int>nex={m+1,1};
		for(int i=1;i<=n;i++){
			if(ans[i].fr)continue;
			int bas=pos.fr/pos.sc+1;
			int l=bas,r=m;
			while(l<r){
				int mi=(l+r)/2;
				int o=arr[i][bas]-arr[i][bas-1];
				int p1=pos.sc-(pos.fr%pos.sc);
				int p2=pos.sc;
				int x=arr[i][mi]-arr[i][bas];
				if(__int128_t(o*p1+x*p2)*n>=sum[i]*p2)r=mi;
				else l=mi+1;
			}
			__int128_t x=arr[i][l-1]-arr[i][bas];
			__int128_t p1=pos.sc-(pos.fr%pos.sc),p2=pos.sc;
			__int128_t t=arr[i][l]-arr[i][l-1],o=arr[i][bas]-arr[i][bas-1];
			pair<__int128_t,int>can={sum[i]*p2-x*n*p2-o*n*p1,n*t*p2};
			can.fr+=__int128_t(l-1)*can.sc;
			int g=gcd(can.fr,can.sc);
			can.fr/=g;
			can.sc/=g;
			if(nex.fr*can.sc>can.fr*nex.sc){
				opt=i;
				nex.fr=can.fr;nex.sc=can.sc;
			}
		}
		if(nex.sc>1e9){
			__int128_t x=(nex.fr*1000000000ll-1)/nex.sc+1;
			assert(x*nex.sc>=nex.fr*1000000000);
			nex.fr=x;
			nex.sc=1e9;
		}
		if(nex.fr>m*nex.sc){
			while(true);
		}
		pos=ans[opt]=nex;
		v.pb(opt);
	}
	int las=v.back();
	v.pop_back();
	for(int x:v){
		if(ans[x].fr==m&&ans[x].sc==1)break;
		cout<<ans[x].fr<<" "<<ans[x].sc<<endl;
	}
	v.pb(las);
	for(int x:v){
		cout<<x<<" ";
	}
	cout<<endl;
}
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...