主页 / 问题 / 1317001
Accepted
Stanley
Asked:2022-08-13 00:15:59 +0000 UTC2022-08-13 00:15:59 +0000 UTC

E-Olymp:多边形区域

  • 772

问题链接 我用梯形公式解决了这个问题，但是有些测试是不正确的，我不知道错误在哪里。

#include <stdio.h>  
#define MAX 1002  
int x[MAX], y[MAX],n;  
double findArea(int *x, int *y)  
{  
double s = 0;  
x[n] = x[0]; y[n] = y[0];  
for(int i = 0; i < n; i++)  
s += (y[i+1] + y[i]) * (x[i+1] - x[i])/2.0;  
return (s < 0) ? -s : s;  
}  
int main(void){  
scanf("%d",&n);  
for(int i = 0; i < n; i++){  
scanf("%d %d",&x[i],&y[i]);}  
printf("%.3lf\n",findArea(x,y));  
}
c++

2 个回答

  1. 在顶点条件下最多可以有 50000 个，而您只为 1000 个分配内存可以吗？

    她在这里坠落……

    顺便说一句，我会做这样的事情：

    #include <stdio.h>
#define MAX 50005
int x[MAX], y[MAX],n;

double findArea()
{
    long long s = 0;
    for(int i = 0; i < n; i++)
        s += (y[i+1] + y[i]) * (x[i+1] - x[i]);
    if (s < 0) s = -s;
    return s/2.0;
}
int main(void){
    scanf("%d",&n);
    for(int i = 0; i < n; i++){
        scanf("%d %d",&x[i],&y[i]);}
    x[n] = x[0]; y[n] = y[0];
    printf("%.3lf\n",findArea());
}

    但这并不重要...

    • 2
  2. Best Answer

    为此，有一个非常简单的多边形面积公式，处理异常就足够了，而且比梯形效果更好。这就是高斯面积公式。

    #include <iostream>
#include <iomanip>
#include <cmath>
#include <vector>
#include <set>
#include <string>
#include <algorithm>
#include <cctype>
#include <iterator>
#include <climits>
#include <complex>
#include <list>
#include <map>

using namespace std;
struct Vecs{
    double x, y;
};
struct Point{
    double x, y;
};

Vecs  makeVector(Point& a, Point & b){
    Vecs res ;
    res.x = b.x - a.x;
    res.y = b.y - a.y;
    return res;
}
long double dot(Vecs& a, Vecs& b){
    return a.x * b.x + a.y * b.y;
}
long double len(Vecs& a){
    return sqrt(a.x*a.x+a.y*a.y);
}
long double cos(Vecs& a, Vecs& b){
    return dot(a, b)/(len(a) * len(b));
}
inline void solve(){
    size_t size;
    int j;
    cin >> size;
    vector <Point> base;
    base.resize(size);
    for (int i = 0; i < size; i++) {
        cin >> base[i].x >> base[i].y;
    }
    long double ans_sum=0;
    for (int i = 0; i < size; i++) {
        j = i + 1;
        if (j >= size) j = 0;
        ans_sum += (base[i].x * base[j].y) - (base[j].x * base[i].y);
    }
    ans_sum = abs(ans_sum);
    ans_sum /= 2;
    cout.precision(10);
    cout << fixed;
    cout <<ans_sum << endl;
    return;
}
int main() {
    solve();
    return 0;
}
    • 2

相关问题