51nod 1719 数值计算

链接

分析：都想出来方法了，结果怎么都不对，还去看了题解，结果发现自己犯蠢。。

这个式子看起来就是能够化简得，然后周期是$\pi$。
我是怎么看出来周期是$\pi$，不瞒你，我也是猜的。。一开始我想，这个函数周期一定是$2\pi$，然后代进去发现答案比样例大了，然后改成$\pi$，差不多。。那么周期就是$\pi$了。。然后怎么算第一个区间呢，就二分吧。
于是在只要$[0,\pi]$这个区间，算出第一个值，题目就迎刃而解了。而在这个区间里，函数是单调的，单调性不知，通过第一个值自己判断一下就行了。

结果我定义计算值函数把mid定义成了int，于是怎么也二分不对:(，真是太蠢了

代码：

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

using namespace std;

#define   offcin        ios::sync_with_stdio(false)
#define   sigma_size    26
#define   lson          l,m,v<<1
#define   rson          m+1,r,v<<1|1
#define   slch          v<<1
#define   srch          v<<1|1
#define   sgetmid       int m = (l+r)>>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   bits(a)       __builtin_popcount(a)
#define   mk            make_pair
#define   pb            push_back
#define   fi            first
#define   se            second

const int    INF    = 0x3f3f3f3f;
const ll     INFF   = 1e18;
const double pi     = acos(-1.0);
const double inf    = 1e18;
const double eps    = 1e-6;
const ll     mod    = 1e9+7;
const int    maxmat = 10;
const ull    BASE   = 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';
 }
/*****************************************************/

double calc(int a, int b, double mid) {
    double res = 0;
    for (int i = 1; i <= 1e4; i ++) {
        res += (double)a * sin(mid + i) / (i + sin(i));
        res += (double)b * cos(mid + i) / (i + cos(i));
    }
    return res;
}

int main(int argc, char const *argv[]) {
    int a, b, n;
    cin>>a>>b>>n;
    double tmp = calc(a, b, 0);
    double lb = 0, ub = pi, res = 0;
    for (int i = 0; i < 100; i ++) {
        double mid = (lb + ub) / 2;
        double temp = calc(a, b, mid);
        if (fabs(temp) < eps) {res = mid; break; }
        if (temp < eps){
            if (tmp < eps) lb = mid;
            else ub = mid;
        }
        else {
            if (tmp < eps) ub = mid;
            else lb = mid;
        }
        // res = mid;
    }
    // cout<<res<<endl;
    res = res * n + (double)n * (n - 1) / 2 * pi;
    printf("%.3f\n", res);
    return 0;
}