#include"prg.h"
#include<vector>
#include<cmath>
#include<ctime>
#include<cstdlib>
using namespace std;
double solve(int n,int m,double p)
{
	double a=1,b=1,c=1,d=1;
	for(int i=0;i<m;i++){
		a*=p; c*=1-p;
	}
	for(int i=0;i<n-m;i++){
		b*=1-p; d*=p;
	}
	return log2(a*b+c*d);
}

int solve(int n,int m,int k,vector<int> op,array<vector<int>,3000> x,int c,array<array<int,3000>,25> s)
{
	int cnt=0;
	for(int i=0;i<25;i++)
		for(int j=0;j<m;j++)
			cnt+=s[i][j];
	if(cnt<35000||cnt>40000) return 2;
	if(op[k-2]!=2) return 1;
	if(op[1]!=2){
		double p1,p2;
		if(op[0]==2){
			p1=0.5; p2=0.75;
		}
		else if(op[0]==op[1]){
			p1=0.5; p2=0.875;
		}
		else{
			p1=0.5; p2=0.625;
		}
		double E1=0.0,E2=0.0;
		int real=0;
		vector<int> pos[1000];
		for(int i=0;i<m;i++)
			pos[x[i][k-1]].push_back(i);
		for(int i=0;i<n;i++){
			int t=pos[i].size();
			for(int d=0;d<25;d++){
				int sum=0;
				for(int j=0;j<t;j++)
					sum+=s[d][pos[i][j]];
				E1+=solve(t,sum,p1);
				E2+=solve(t,sum,p2);
			}
		}
		if(E2>E1) return 2;
		return 1;
	}
/*	vector<int> ed[1000];
	vector<int> id[1000];
	for(int i=0;i<m;i++){
		int a=x[i][2],b=x[i][3];
		ed[a].push_back(b);
		ed[b].push_back(a);
		id[a].push_back(i);
		id[b].push_back(i);
	}
	srand(time(NULL));
	int opt[1000],perm[1000];
		int sum=0;
	for(int w=0;w<25;w++){
		for(int i=0;i<n;i++){
			perm[i]=i; opt[i]=-1;
		}
		for(int d=0;d<2;d++){
			for(int i=0;i<n;i++)
				swap(perm[rand()%n],perm[rand()%n]);
				int i,a,b,to,u;
			for(int j=0;j<n;j++){
				a=b=0; i=perm[j];
				for(int t=0;t<ed[i].size();t++){
					u=id[i][t]; to=ed[i][t];
					if(opt[to]==(s[w][u]^op[0])) a++;
					if(opt[to]==(s[w][u]^op[0]^1)) b++;
				}
				if(op[0]==1){
					if(a>b) opt[i]=0;
					if(a<b) opt[i]=1;
					if(a==b) opt[i]=rand()&1;
				}
				if(op[0]==0){
					if(a>b) opt[i]=1;
					if(a<b) opt[i]=0;
					if(a==b) opt[i]=rand()&1;
				}
			}
		}
		for(int i=0;i<m;i++){
			int a=x[i][2],b=x[i][3];
			if(opt[a]==opt[b]) sum+=s[w][i];
		}
	}
		printf("%d\n",sum);*/
	return 1;
}