uva 1118(圆中三角形)

2023-12-27 02:48
文章标签 三角形 uva 1118 圆中

本文主要是介绍uva 1118(圆中三角形),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

我一开始就是想着求三角形面积,根据三个点之间任意两个点与0点形成的三角形加减来求得;

但是这样写是n3我以为不会超时,但是还是te了;


后来发现了网上思路的升级版本(n2解决):

画一个图看看



在上图中。 我们假定要求 S(i,j,k) 如果k在 j的左侧  那么 S(i,j,k) = S(i,k,o)+s(i,j,o) - s(k,i,o); 显然 只要k在j的左侧  s(i,j,0) 在用来求 做和用的。

如果k在 j的右侧  那么 S(i,j,k) = S(i,k,o)+s(i,j,o) - s(k,j,o); 显然 只要k在j的右侧  s(i,j,0) 在用来求 做和用的。

如果k在 i,j的中间  那么 S(i,j,k) = S(j,k,o)+s(i,k,o) - s(i,j,o); 显然 只要k在i,j的中间  s(i,j,0) 在用来求 做差用的。

那么 我们就可以对于 线段 i,j  对于后面 或者前面的点来说 都需要加上oji 这一块面积  而对于  ij 中间的 点来说  都需要 减去这一块。

那么对于这一块面积。  总共的 需求为   当k在左侧的时候 : n - j - 1    k在右侧是时候   为  i   加起来为  n - j + i - 1;

对于k位于中间来说。   总共需要减去 个数为   j - i  - 1个。

那么 总共就是  n + 2i - 2j  个。  乘上  这一块的面积 sin(c)  就OK了


#include <iostream>
#include<string>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<set>
#include<cmath>
using namespace std;
const int maxn=500+5;
const int inf=100000000;
double s[maxn];
const double pai=acos(-1);
double a[5];
double ans;
int main()
{int n,r;double ka;while(scanf("%d%d",&n,&r)&&n&&r){for(int i = 0; i<n; i++)scanf("%lf",&s[i]);sort(s,s+n);ans=0;for(int i = 0; i < n-1; i++)for(int j = i+1; j < n ; j++){double a = s[i]/180.0*pai;double b = s[j]/180.0*pai;ans += (n-j*2+i*2)*sin(b-a);}ans = ans*r*r/2.0;ans += 0.5;long long anss = ans;cout<<anss<<endl;}return 0;
}
叉乘在定义上有:两个向量进行叉乘得到的是一个向量方向垂直于这两个向量构成的平面,大小等于这两个向量组成的平行四边形的面积。


叉乘代码

#include <cstdio>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
const double Pi=acos(-1);
const double eps=1e-8;
struct Point
{double x, y;Point(){}Point(double _x, double _y){x=_x, y=_y;//cout<<x<<" "<<y<<endl;}
}p[505];
double Area( Point A, Point B, Point C )
{return  fabs(((A.x-B.x)*(A.y-C.y)-(A.x-C.x)*(A.y-B.y))/2.0);
}
int n;
double a ,r;
int main( )
{while( scanf("%d%lf", &n, &r)&&n&&r ){// cout<<r<<endl;for(int i=0; i<n; ++ i ){scanf("%lf", &a);p[i]=Point(r*cos(a/180.0*Pi),r*sin(a/180.0*Pi) );}double ans=0;for( int i=0; i<n-2; ++ i )for( int j=i+1;j<n-1; ++ j )for( int k=j+1; k<n; ++ k ){ans+=Area( p[i], p[j], p[k] );// cout<<Area( p[i], p[j], p[k] )<<endl;}printf("%.0f\n",ans);}return 0;
}

wa点:

pai=acos(-1);  3.1425926是wa的

ans四舍五入

强制转化的话long long



原来的te代码:

也不知道,思路是对还是错;

反正计算次数太多te了

#include <iostream>
#include<string>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<set>
#include<cmath>
using namespace std;
const int maxn=500+5;
const int inf=100000000;
double s[maxn];
const double pai=3.1415926;
double get(int i,int j)
{double res=s[i]-s[j];if(res>180.00000)res=360-res;return res;
}double a[5];
double ans;
int main()
{int n,r;double ka;while(scanf("%d%d",&n,&r)&&n&&r){for(int i=0;i<n;i++)scanf("%lf",&s[i]);sort(s,s+n);ans=0;for(int i=0;i<n-2;i++)for(int j=i+1;j<n-1;j++)for(int k=j+1;k<n;k++){a[0]=get(j,i);a[1]=get(k,j);a[2]=get(k,i);if(s[k]-s[i]<180.0000){for(int z=0;z<3;z++){ka=a[z]/2.0;ka=ka/180.0*pai;if(z==2)ans-=(double)r*r*cos(ka)*sin(ka);elseans+=(double)r*r*cos(ka)*sin(ka);}}else if(s[k]-s[j]>180.000){for(int z=0;z<3;z++){ka=a[z]/2.0;ka=ka/180.0*pai;if(z==1)ans-=(double)r*r*cos(ka)*sin(ka);elseans+=(double)r*r*cos(ka)*sin(ka);}}else if(s[j]-s[i]>180.000){for(int z=0;z<3;z++){ka=a[z]/2.0;ka=ka/180.0*pai;if(z==0)ans-=(double)r*r*cos(ka)*sin(ka);elseans+=(double)r*r*cos(ka)*sin(ka);}}else{for(int z=0;z<3;z++){ka=a[z]/2.0;ka=ka/180.0*pai;ans+=(double)r*r*cos(ka)*sin(ka);}}}ans=ans+0.5;long long anss=ans;printf("%lld\n",ans);}return 0;
}



这篇关于uva 1118(圆中三角形)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/541659

相关文章

uva 10055 uva 10071 uva 10300(水题两三道)

情歌两三首,水题两三道。 好久没敲代码了为暑假大作战热热身。 uva 10055 Hashmat the Brave Warrior 求俩数相减。 两个debug的地方,一个是longlong,一个是输入顺序。 代码: #include<stdio.h>int main(){long long a, b;//debugwhile(scanf("%lld%lld", &

poj 3259 uva 558 Wormholes(bellman最短路负权回路判断)

poj 3259: 题意:John的农场里n块地,m条路连接两块地,w个虫洞,虫洞是一条单向路,不但会把你传送到目的地,而且时间会倒退Ts。 任务是求你会不会在从某块地出发后又回来,看到了离开之前的自己。 判断树中是否存在负权回路就ok了。 bellman代码: #include<stdio.h>const int MaxN = 501;//农场数const int

poj 2349 Arctic Network uva 10369(prim or kruscal最小生成树)

题目很麻烦,因为不熟悉最小生成树的算法调试了好久。 感觉网上的题目解释都没说得很清楚,不适合新手。自己写一个。 题意:给你点的坐标,然后两点间可以有两种方式来通信:第一种是卫星通信,第二种是无线电通信。 卫星通信:任何两个有卫星频道的点间都可以直接建立连接,与点间的距离无关; 无线电通信:两个点之间的距离不能超过D,无线电收发器的功率越大,D越大,越昂贵。 计算无线电收发器D

uva 10387 Billiard(简单几何)

题意是一个球从矩形的中点出发,告诉你小球与矩形两条边的碰撞次数与小球回到原点的时间,求小球出发时的角度和小球的速度。 简单的几何问题,小球每与竖边碰撞一次,向右扩展一个相同的矩形;每与横边碰撞一次,向上扩展一个相同的矩形。 可以发现,扩展矩形的路径和在当前矩形中的每一段路径相同,当小球回到出发点时,一条直线的路径刚好经过最后一个扩展矩形的中心点。 最后扩展的路径和横边竖边恰好组成一个直

uva 10061 How many zero's and how many digits ?(不同进制阶乘末尾几个0)+poj 1401

题意是求在base进制下的 n!的结果有几位数,末尾有几个0。 想起刚开始的时候做的一道10进制下的n阶乘末尾有几个零,以及之前有做过的一道n阶乘的位数。 当时都是在10进制下的。 10进制下的做法是: 1. n阶位数:直接 lg(n!)就是得数的位数。 2. n阶末尾0的个数:由于2 * 5 将会在得数中以0的形式存在,所以计算2或者计算5,由于因子中出现5必然出现2,所以直接一

uva 568 Just the Facts(n!打表递推)

题意是求n!的末尾第一个不为0的数字。 不用大数,特别的处理。 代码: #include <stdio.h>const int maxn = 10000 + 1;int f[maxn];int main(){#ifdef LOCALfreopen("in.txt", "r", stdin);#endif // LOCALf[0] = 1;for (int i = 1; i <=

uva 575 Skew Binary(位运算)

求第一个以(2^(k+1)-1)为进制的数。 数据不大,可以直接搞。 代码: #include <stdio.h>#include <string.h>const int maxn = 100 + 5;int main(){char num[maxn];while (scanf("%s", num) == 1){if (num[0] == '0')break;int len =

uva 10014 Simple calculations(数学推导)

直接按照题意来推导最后的结果就行了。 开始的时候只做到了第一个推导,第二次没有继续下去。 代码: #include<stdio.h>int main(){int T, n, i;double a, aa, sum, temp, ans;scanf("%d", &T);while(T--){scanf("%d", &n);scanf("%lf", &first);scanf

uva 10916 Factstone Benchmark(打表)

题意是求 k ! <= 2 ^ n ,的最小k。 由于n比较大,大到 2 ^ 20 次方,所以 2 ^ 2 ^ 20比较难算,所以做一些基础的数学变换。 对不等式两边同时取log2,得: log2(k ! ) <=  log2(2 ^ n)= n,即:log2(1) + log2(2) + log2 (3) + log2(4) + ... + log2(k) <= n ,其中 n 为 2 ^

uva 10025 The ? 1 ? 2 ? ... ? n = k problem(数学)

题意是    ?  1  ?  2  ?  ...  ?  n = k 式子中给k,? 处可以填 + 也可以填 - ,问最小满足条件的n。 e.g k = 12  - 1 + 2 + 3 + 4 + 5 + 6 - 7 = 12 with n = 7。 先给证明,令 S(n) = 1 + 2 + 3 + 4 + 5 + .... + n 暴搜n,搜出当 S(n) >=