poj1654 Area

来源:互联网 发布:重庆时尚频道直播软件 编辑:程序博客网 时间:2024/05/21 17:39

Area
Time Limit: 1000MS Memory Limit: 10000KTotal Submissions: 17498 Accepted: 4855

Description

You are going to compute the area of a special kind of polygon. One vertex of the polygon is the origin of the orthogonal coordinate system. From this vertex, you may go step by step to the following vertexes of the polygon until back to the initial vertex. For each step you may go North, West, South or East with step length of 1 unit, or go Northwest, Northeast, Southwest or Southeast with step length of square root of 2.

For example, this is a legal polygon to be computed and its area is 2.5:

Input

The first line of input is an integer t (1 <= t <= 20), the number of the test polygons. Each of the following lines contains a string composed of digits 1-9 describing how the polygon is formed by walking from the origin. Here 8, 2, 6 and 4 represent North, South, East and West, while 9, 7, 3 and 1 denote Northeast, Northwest, Southeast and Southwest respectively. Number 5 only appears at the end of the sequence indicating the stop of walking. You may assume that the input polygon is valid which means that the endpoint is always the start point and the sides of the polygon are not cross to each other.Each line may contain up to 1000000 digits.

Output

For each polygon, print its area on a single line.

Sample Input

4582567256244865

Sample Output

000.52



将走的每一步转换成坐标点,然后对坐标点用向量叉乘求面积,因为在网格中走,所以用整形存就行了

有个坑点是中间算总面积求和的ans要用long long不然会整形溢出



#include <cstdio>#include <iostream>#include <cstring>#include <cmath>using namespace std;struct Point {int x, y;Point() {}Point(int x, int y) : x(x), y(y) {}};typedef Point Vector;char ch[1000010];Point p[1000010];int Cross(Vector a, Vector b) {return a.x * b.y - a.y * b.x;}void Area() {p[0] = Point(0, 0);int t = 0;for (int i = 0; ch[i] != '5'; i++) {++t;switch (ch[i]) {case '1': p[t] = Point(p[t - 1].x - 1, p[t - 1].y - 1); break;case '2': p[t] = Point(p[t - 1].x, p[t - 1].y - 1); break;case '3': p[t] = Point(p[t - 1].x + 1, p[t - 1].y - 1); break;case '4': p[t] = Point(p[t - 1].x - 1, p[t - 1].y); break;case '6': p[t] = Point(p[t - 1].x + 1, p[t - 1].y); break;case '7': p[t] = Point(p[t - 1].x - 1, p[t - 1].y + 1); break;case '8': p[t] = Point(p[t - 1].x, p[t - 1].y + 1); break;case '9': p[t] = Point(p[t - 1].x + 1, p[t - 1].y + 1); break;}}long long ans = 0;   //要使用long long,不然溢出for (int i = 0; i < t; i++) {ans += Cross(p[i], p[i + 1]);}ans = ans < 0 ? -ans : ans;ans % 2 == 0 ? printf("%I64d", ans / 2) : printf("%I64d.5", ans / 2);return;}int main(){int t;scanf("%d", &t);while (t--) {scanf("%s", ch);int len = strlen(ch);if (len <= 3)printf("0");elseArea();cout << endl;}return 0;}



0 0
原创粉丝点击