51nod 1719 数值计算(二分)

1719 数值计算

基准时间限制：1 秒 空间限制：131072 KB 分值: 80 难度：5级算法题 收藏 关注

令 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位

Input

三个数 A B n

Output

答案 保留3位小数

Input示例

5 4 1000

Output示例

1572338.870

一个很复杂的式子，然后不会做。

后来才知道，其实里面的k，按照那个诱导公式展开

得到F(x)=Asinx+Bcosx=Csin(x+a)

然后就是每π个必定有一个解，所以只要求第一个最小正数解

所以第一个正解肯定在[0,π]之间，然后并且这个函数在这个之间是单调的，所以可以二分

先判0点的时候值是0还是正还是负，然后二分搞就行了

注意精度，1e−8可以过

代码：

#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cmath>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <sstream>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#pragma comment(linker,"/STACK:102400000,102400000")

using namespace std;
#define   MAX           100005
#define   MAXN          1000005
#define   maxnode       105
#define   sigma_size    30
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   lrt           rt<<1
#define   rrt           rt<<1|1
#define   middle        int m=(r+l)>>1
#define   LL            long long
#define   ull           unsigned long long
#define   mem(x,v)      memset(x,v,sizeof(x))
#define   lowbit(x)     (x&-x)
#define   pii           pair<int,int>
#define   bits(a)       __builtin_popcount(a)
#define   mk            make_pair
#define   limit         10000

//const int    prime = 999983;
const int    INF   = 0x3f3f3f3f;
const LL     INFF  = 0x3f3f;
const double pi    = acos(-1.0);
//const double inf   = 1e18;
const double eps   = 1e-8;
const LL     mod   = 1e9+7;
const ull    mx    = 133333331;

/*****************************************************/
inline void RI(int &x) {
char c;
while((c=getchar())<'0' || c>'9');
x=c-'0';
while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
}
/*****************************************************/

int main(){
double a,b,n;
scanf("%lf%lf%lf",&a,&b,&n);
double l=0,r=pi;
double tmp=0;
for(int i=1;i<=10000;i++){
tmp+=a/(i+sin(i))*sin(i)+b/(i+cos(i))*cos(i);
}
double ans;
if(abs(tmp)<eps) ans=pi;
else{
for(int t=0;t<100;t++){
double mid=(l+r)/2;
double temp=0;
for(int i=1;i<=10000;i++){
temp+=a/(i+sin(i))*sin(mid+i)+b/(i+cos(i))*cos(mid+i);
}
if(abs(temp)<eps){
ans=mid;
break;
}
if(tmp<0){
if(temp<0) l=mid;
else r=mid;
}
else{
if(temp<0) r=mid;
else l=mid;
}
}
}
printf("%.3f\n",n*ans+n*(n-1)*pi/2);
return 0;
}