LA6142 Probability Through Experiments 统计圆上的锐角三角形

题目地址:LA6142

先说统计思路: 我们统计钝角和直角三角形,这样才不会重复计算,每次找到一个点后,统计这个逆时针走的这段半圆弧上的点的个数,任取两个就是钝角或者直角三角形,而且补重不漏。

找个数--二分,用upperbound 这样nlogn时间复杂度   暴力n*n*n肯定不行

注意 1 读入数据的时候,double *1000不行啊,浮点数误差+0.5 取floor  为这个wa了很久

         2 180度这个点要单独考虑,

代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;

int p[20005];

int main()
{
    long long  n,r;
    
    double temp;
    
    
    int cas=0;
    while(cin>>n>>r)
    {
        if(!n)  break;
        long long all=n*(n-1)*(n-2)/6;
        
        long long ans=0;
        for(int i=0;i<n;i++)
        {
            scanf("%lf",&temp);
            p[i]=floor(temp*1000+0.5);  
        }
        
        sort(p,p+n);
        
        
        for(int i=0;i<n;i++)
        {
            if(p[i]<180000)
            {
                long long temp=upper_bound(p, p+n, p[i]+180000)-lower_bound(p,p+n, p[i])-1;
                 ans+=temp*(temp-1)/2;
                
                //cout<<'i'<<i<<"temp"<<temp<<endl;
            }
            else if(p[i]==180000)
            {
                long long b=n-i-1;
                if(p[0]==0) b++;
                ans+=b*(b-1)/2;
            }
            else
            {
                long long a=n-i-1;
                long long b=upper_bound(p, p+n, p[i]-180000)-p+a;
                
                ans+=b*(b-1)/2;
                
               // cout<<'i'<<i<<"b"<<b<<endl;
                
            }
        }
        

        
        
        printf("Case %d: %lld\n",++cas,all-ans);
        
        
       
    }
}



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值