#include "tickets.h"
#include <vector>
#include <algorithm>
using namespace std;
struct A
{
int where;
int con1,con2;
}tt[1505];
vector < vector < int > > answer;
vector < int > row;
vector< vector<int> > all;
vector < pair < int , pair < int , int > > > how;
vector < int > uu[1505];
bool have[1505][1505];
bool use[1505]={0};
long long DP[1505][1505];
long long con[1505][1505];
int a[1505];
int b[1505];
int N,M;
long long F(int here,int con)
{
if(con>N/2) return -1e18;
if(here==N&&con!=N/2) return -1e18;
if(here==N) return 0;
if(have[here][con]) return DP[here][con];
have[here][con]=1;
DP[here][con]=max(F(here+1,con+1)+all[here][M-1],F(here+1,con)-all[here][0]);
return DP[here][con];
}
void F2(int here,int con)
{
if(here==N) return;
//printf("%d %d\n",F(here+1,con+1)+all[here][M-1],F(here+1,con)-all[here][0]);
if(F(here+1,con+1)+all[here][M-1]>F(here+1,con)-all[here][0])
{
//printf("aa\n");
F2(here+1,con+1);
answer[here][M-1]=0;
}
else
{
F2(here+1,con);
answer[here][0]=0;
}
}
bool cmp1(A a,A b)
{
return a.con1>b.con1;
}
bool cmp2(A a,A b)
{
return a.con2>b.con2;
}
long long find_maximum(int k,vector< vector<int> > x)
{
all=x;
N=x.size();
M=x[0].size();
int i,j,l,t,ok,now,ll,rr;
long long ans=0;
for(i=0;i<M;i++) row.push_back(-1);
for(i=0;i<N;i++)
{
answer.push_back(row);
for(j=0;j<M;j++) how.push_back(make_pair(x[i][j],make_pair(i,j)));
}
if(k==1)
{
ans=F(0,0);
F2(0,0);
}
else if(k==M)
{
sort(how.begin(),how.end());
for(i=0;i<N*M/2;i++) ans-=how[i].first;
for(i=N*M/2;i<N*M;i++) ans+=how[i].first;
for(i=0;i<N*M;i++) con[how[i].second.first][how[i].second.second]=i;
for(i=0;i<N;i++)
{
a[i]=0;
b[i]=M-1;
}
for(i=0;i<k;i++)
{
ll=0;
rr=0;
for(j=0;j<N;j++) use[j]=0;
for(j=0;j<N;j++)
{
if(con[j][a[j]]>=N*M/2)
{
use[j]=1;
rr++;
answer[j][b[j]]=i;
b[j]--;
}
if(con[j][b[j]]<N*M/2)
{
use[j]=1;
ll++;
answer[j][a[j]]=i;
a[j]++;
}
}
for(j=0;j<N;j++)
{
if(use[j]) continue;
if(ll<N/2)
{
use[j]=1;
ll++;
answer[j][a[j]]=i;
a[j]++;
}
else
{
use[j]=1;
rr++;
answer[j][b[j]]=i;
b[j]--;
}
}
}
}
else
{
//allocate_tickets(answer);
ll=0;
rr=0;
t=0;
for(i=0;i<N;i++)
{
a[i]=0;
b[i]=k-1;
t=0;
for(j=0;j<M;j++) if(x[i][j]==0) t++;
ll+=min(k,t);
tt[i].con1=min(k,t);
t=0;
for(j=0;j<M;j++) if(x[i][j]==1) t++;
rr+=min(k,t);
tt[i].con2=min(k,t);
tt[i].where=i;
}
ans=min(min(ll,rr),N*k/2);
//printf("%d %d\n",ll,rr);
//for(i=0;i<N;i++) for(j=0;j<k;j++) answer[i][j]=j;
if(ll<rr)
{
ll=ans;
sort(tt,tt+N,cmp1);
for(i=0;i<N;i++)
{
t=min(tt[i].con1,ll);
//printf("%d %d %d\n",tt[i].where,tt[i].con1,t);
for(j=0;j<t;j++) uu[tt[i].where].push_back(j);
for(j=0;j<k-t;j++) uu[tt[i].where].push_back(M-j-1);
ll-=t;
}
}
else
{
rr=ans;
sort(tt,tt+N,cmp2);
for(i=0;i<N;i++)
{
t=min(tt[i].con2,rr);
for(j=0;j<t;j++)
{
//printf("%d %d\n",i,j);
uu[tt[i].where].push_back(M-j-1);
}
for(j=0;j<k-t;j++)
{
//printf("%d %d\n",i,j);
uu[tt[i].where].push_back(j);
}
rr-=t;
}
}
how.clear();
for(i=0;i<N;i++)
{
for(auto j:uu[i])
{
//printf("%d %d\n",i,j);
how.push_back(make_pair(x[i][j],make_pair(i,j)));
}
}
sort(how.begin(),how.end());
for(i=0;i<N*k;i++)
{
//printf("%d %d\n",how[i].second.first,how[i].second.second);
con[how[i].second.first][how[i].second.second]=i;
}
for(i=0;i<N;i++)
{
a[i]=0;
b[i]=k-1;
}
for(i=0;i<k;i++)
{
ll=0;
rr=0;
for(j=0;j<N;j++) use[j]=0;
for(j=0;j<N;j++)
{
if(con[j][uu[j][a[j]]]>=N*k/2)
{
use[j]=1;
rr++;
answer[j][uu[j][b[j]]]=i;
b[j]--;
}
if(con[j][uu[j][b[j]]]<N*k/2)
{
use[j]=1;
ll++;
answer[j][uu[j][a[j]]]=i;
a[j]++;
}
}
for(j=0;j<N;j++)
{
if(use[j]) continue;
if(ll<N/2)
{
use[j]=1;
ll++;
answer[j][uu[j][a[j]]]=i;
a[j]++;
}
else
{
use[j]=1;
rr++;
answer[j][uu[j][b[j]]]=i;
b[j]--;
}
}
}
}
allocate_tickets(answer);
return ans;
}
