Codeforces 38D-Vasya the Architect

Vasya the Architect

time limit per test

2 seconds

memory limit per test

256 megabytes

input

standard input

output

standard output

Once Vasya played bricks. All the bricks in the set had regular cubical shape. Vasya vas a talented architect, however the tower he built kept falling apart.

Let us consider the building process. Vasya takes a brick and puts it on top of the already built tower so that the sides of the brick are parallel to the sides of the bricks he has already used. Let's introduce a Cartesian coordinate system on the horizontal plane, where Vasya puts the first brick. Then the projection of brick number i on the plane is a square with sides parallel to the axes of coordinates with opposite corners in points (xi, 1, yi, 1) and (xi, 2, yi, 2). The bricks are cast from homogeneous plastic and the weight of a brick a × a × a is a3grams.

It is guaranteed that Vasya puts any brick except the first one on the previous one, that is the area of intersection of the upper side of the previous brick and the lower side of the next brick is always positive.

We (Vasya included) live in a normal world where the laws of physical statics work. And that is why, perhaps, if we put yet another brick, the tower will collapse under its own weight. Vasya puts the cubes consecutively one on top of the other until at least one cube loses the balance and falls down. If it happens, Vasya gets upset and stops the construction. Print the number of bricks in the maximal stable tower, that is the maximal number m satisfying the condition that all the towers consisting of bricks 1, 2, ..., k for every integer k from 1 to mremain stable.

Input

The first input file contains an integer n (1 ≤ n ≤ 100) which is the number of bricks. Each of the next n lines contains four numbers xi, 1, yi, 1, xi, 2, yi, 2 (xi, 1 ≠ xi, 2, |xi, 1 - xi, 2| = |yi, 1 - yi, 2|) which are the coordinates of the opposite angles of the base of the brick number i. The coordinates are integers and their absolute value does not exceed 50.

The cubes are given in the order Vasya puts them. It is guaranteed that the area of intersection of the upper side of the brick number i - 1and the lower side of the brick number i is strictly strictly greater than zero for all i ≥ 2.

Output

Print the number of bricks in the maximal stable tower.

Examples

input

Copy

2
0 0 3 3
1 0 4 3

output

2

input

Copy

2
0 0 3 3
2 0 5 3

output

1

input

Copy

3
0 0 3 3
1 0 4 3
2 0 5 3

output

3

题意：给出n个块的左下角坐标和右上角坐标，按照给出的块的顺序一个个往上叠，叠不了了就结束，问最多能叠多少个

解题思路：判断一个块能否保持平衡的方法是看它的重心在xoy平面上的投影是否落在它下方的那一个块上（包括边界）。对于单个brick，计算重心坐标是一个显而易见的问题，假设对角顶点在xoy上的投影的坐标分别是（x1，y1）,（x2，y2），则重心在xoy上的投影的坐标就是两者的中点（（x1 + x2）/2，（y1 + y2）/2）。多个brick组成的复合体的重心：Pw = (P[1] * W[1] + P[2] * W[2] + ... + P[n] * W[n]) / (W[1] + W[2] + ... + W[n])，比如要求y，则Pw是复合体的重心的纵坐标，p[i]是复合体中brick[i]的重心纵坐标，w[i]是brick[i]的重量

#include <iostream>  
#include <cstdio>  
#include <cstring>  
#include <string>  
#include <algorithm>  
#include <map>  
#include <set>
#include <stack>  
#include <queue>  
#include <vector>  
#include <bitset>  
#include <functional>  

using namespace std;

#define LL long long
const int INF = 0x3f3f3f3f;

struct node
{
	int x1, x2, y1, y2; 
	double xx, yy, w;
}a[101];

int n;

int main()
{ 
	while (~scanf("%d", &n))
	{
		for (int i = 1; i <= n; i++)
		{
			scanf("%d %d %d %d", &a[i].x1, &a[i].y1, &a[i].x2, &a[i].y2);
			a[i].w = pow(((a[i].x1 - a[i].x2) * (a[i].x1 - a[i].x2) + (a[i].y1 - a[i].y2) * (a[i].y1 - a[i].y2)) / 2.0, 1.5);
			a[i].xx = 0.5 * (a[i].x1 + a[i].x2);
			a[i].yy = 0.5 * (a[i].y1 + a[i].y2);
		}
		int flag = 0;
		for (int i = 2; i <= n; i++)
		{
			double xx = a[i].xx, yy = a[i].yy, w = a[i].w;
			for (int j = i - 1; j > 0; j--)
			{
				if ((xx - a[j].x1) * (xx - a[j].x2) > 0 || (yy - a[j].y1) * (yy - a[j].y2) > 0)
				{
					flag = 1;
					printf("%d\n", i - 1);
					break;
				}
				xx = (w * xx + a[j].w * a[j].xx) / (w + a[j].w);
				yy = (w * yy + a[j].w * a[j].yy) / (w + a[j].w);
				w += a[j].w;
			}
			if (flag) break;
		}
		if (!flag) printf("%d\n", n);
	}
	return 0;
}