51nod1719数值计算+二分

F(x)=wk=1(Ak+sin(k)sin(x+k)+Bk+cos(k)cos(x+k))
求F(x)=0的前n小的正根的和。
n<=3e6,A<=1e3,B<=1e3
其中w是定值,为1e4
保留到小数点后3位
因为 Asin(x)+Bcos(x) 可以写成 C+sin(x+D) 的形式..so,,上边的那个式子的周期为2pi咯?零点可以是ans0+pi咯。。。。
那就在0-pi中二分出零点,那么之后的零点很快就搞好咯。。

#include<bits/stdc++.h>
using namespace std;
#define pi acos(-1.0)
#define eps 1e-8
int A,B,n;
double zero(double x){
    double ans=0;
    for(int i=1;i<=10000;i++){
        ans+=A*sin(x+i)/(i+sin(i))+B*cos(x+i)/(i+cos(i));
    }
    return ans;
}

int main(){
    scanf("%d %d %d",&A,&B,&n);
    double l=0,r=pi;
    while(r-l>eps){
        double mid=(l+r)/2;
        if(zero(mid)>0) r=mid;
        else l=mid;
    }
    double firstans=l;
    double nextans=l+pi;
    double ans=firstans;
    for(int i=2;i<=n;i++){
        ans+=nextans;
        nextans+=pi;
    }
    printf("%.3f\n",ans);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值