好久没做过三分的题目了。

  1 /* 4355 */
  2 #include <iostream>
  3 #include <sstream>
  4 #include <string>
  5 #include <map>
  6 #include <queue>
  7 #include <set>
  8 #include <stack>
  9 #include <vector>
 10 #include <deque>
 11 #include <algorithm>
 12 #include <cstdio>
 13 #include <cmath>
 14 #include <ctime>
 15 #include <cstring>
 16 #include <climits>
 17 #include <cctype>
 18 #include <cassert>
 19 #include <functional>
 20 #include <iterator>
 21 #include <iomanip>
 22 using namespace std;
 23 //#pragma comment(linker,"/STACK:102400000,1024000")
 24
 25 #define sti                set<int>
 26 #define stpii            set<pair<int, int> >
 27 #define mpii            map<int,int>
 28 #define vi                vector<int>
 29 #define pii                pair<int,int>
 30 #define vpii            vector<pair<int,int> >
 31 #define rep(i, a, n)     for (int i=a;i<n;++i)
 32 #define per(i, a, n)     for (int i=n-1;i>=a;--i)
 33 #define clr                clear
 34 #define pb                 push_back
 35 #define mp                 make_pair
 36 #define fir                first
 37 #define sec                second
 38 #define all(x)             (x).begin(),(x).end()
 39 #define SZ(x)             ((int)(x).size())
 40 #define lson            l, mid, rt<<1
 41 #define rson            mid+1, r, rt<<1|1
 42
 43 typedef struct {
 44     double x, w;
 45 } node_t;
 46
 47 const double eps = 1e-4;
 48 const int maxn = 50005;
 49 node_t nd[maxn];
 50 int n;
 51
 52 double calc(double x) {
 53     double ret = 0.0;
 54
 55     rep(i, 0, n)
 56         ret += pow(fabs(x-nd[i].x), 3) * nd[i].w;
 57
 58     return ret;
 59 }
 60
 61 void solve() {
 62     double l = 1e7;
 63     double r = -1e7;
 64     double m1, m2;
 65     double tmp1, tmp2;
 66
 67     rep(i, 0, n) {
 68         l = min(l, nd[i].x);
 69         r = max(r, nd[i].x);
 70     }
 71
 72     while (fabs(r-l) > eps) {
 73         m1 = l + (r-l)/3;
 74         m2 = r - (r-l)/3;
 75         tmp1 = calc(m1);
 76         tmp2 = calc(m2);
 77         if (tmp1 < tmp2) {
 78             r = m2;
 79         } else if (tmp1 > tmp2) {
 80             l = m1;
 81         } else {
 82             l = r;
 83         }
 84     }
 85
 86     double ans = calc(l);
 87     printf("%.0lf\n", ans);
 88 }
 89
 90 int main() {
 91     ios::sync_with_stdio(false);
 92     #ifndef ONLINE_JUDGE
 93         freopen("data.in", "r", stdin);
 94         freopen("data.out", "w", stdout);
 95     #endif
 96
 97     int t;
 98
 99     scanf("%d", &t);
100     rep(tt, 1, t+1) {
101         scanf("%d", &n);
102         rep(i, 0, n)
103             scanf("%lf %lf", &nd[i].x, &nd[i].w);
104         printf("Case #%d: ", tt);
105         solve();
106     }
107
108     #ifndef ONLINE_JUDGE
109         printf("time = %d.\n", (int)clock());
110     #endif
111
112     return 0;
113 }

转载于:https://www.cnblogs.com/bombe1013/p/5188380.html

【HDOJ】4355 Party All the Time相关推荐

  1. 【HDOJ】4602 Partition

    [题目]http://acm.hdu.edu.cn/showproblem.php?pid=4602 [报告] 直接贴上标程解题报告:(虽然有些纠结,试一下就弄通了..) Problem C. Par ...

  2. 【HDOJ】4602 Partition_天涯浪子_新浪博客

    [题目]http://acm.hdu.edu.cn/showproblem.php?pid=4602 [报告] 直接贴上标程解题报告:(虽然有些纠结,试一下就弄通了..) Problem C. Par ...

  3. 【HDOJ】1008 Elevator_天涯浪子_新浪博客

    [题目] http://acm.hdu.edu.cn/showproblem.php?pid=1008 [报告] 这是我在HDOJ里看到的最水的题目,至少目前为止. 题目里已经讲得很清楚了,模拟就能过 ...

  4. 【HDOJ】1022 Train Problem I_天涯浪子_新浪博客

    [题目]http://acm.hdu.edu.cn/showproblem.php?pid=1022 [报告] 模拟,直接模拟一个栈的运行就行了. [程序] // Task: 1022 Train P ...

  5. 【HDOJ】1003 Max Sum_天涯浪子_新浪博客

    [题目]http://acm.hdu.edu.cn/showproblem.php?pid=1003 [报告] 既然是子区间加和问题,可以用减法来处理.比如求[5..10]的和,可以直接用[1..10 ...

  6. 【HDOJ】4768 Flyer_天涯浪子_新浪博客

    [题目]http://acm.hdu.edu.cn/showproblem.php?pid=4768 [报告] 题目中已经说了,保证最多出现一个奇数,所以可以用异或的想法.把所有的数异或起来,如果最后 ...

  7. 【HDOJ】4699 Editor_天涯浪子_新浪博客

    [题目] http://acm.hdu.edu.cn/showproblem.php?pid=4699 [报告] 模拟题吧,本身不是很难. 最初的想法是用Splay tree写,然后疯掉了.苏牛说2B ...

  8. 【HDOJ】4704 Sum_天涯浪子_新浪博客

    [题目] http://acm.hdu.edu.cn/showproblem.php?pid=4704 [报告] S(K)显然就是N的K正整数划分数.所以SUM{S(K),1<=K<=N} ...

  9. 【HDOJ】1015 Safecracker_天涯浪子_新浪博客

    [题目]http://acm.hdu.edu.cn/showproblem.php?pid=1015 [报告] 题目写的很长,忽略一系列废话,题意就是求一个方程v - w^2 + x^3 - y^4 ...

最新文章

  1. 《C#精彩实例教程》小组阅读12 -- C#面向对象技术高级应用
  2. python读取单元格的值csv,比较单元格值的csv文件python
  3. 视频 + PPT | 用户中台建设实践解析
  4. Struts初步知识
  5. 由浅到深理解ROS(5.1)- roslaunch 学习
  6. IIS执行ASP程序不能调用数据库出现500错误的解决方法
  7. 10深入Redis Cluster
  8. 条件注释判断浏览器版本
  9. Java——异常和断言
  10. java 天猫项目_2018最新Java前后端练手项目【模仿天猫前端】
  11. 计算机网络基础肖盛文电子书,网络实用教程
  12. 十大工业设计公司是如何排名的?
  13. 不想被甲方折磨?移动互联网的9大趋势知道一个算你牛
  14. st7920驱动OCMJ2X8C屏使用CGRAM自定义图标
  15. 184页10万字智慧城市公共服务中台:业务和数据中台建设方案
  16. Linux相关教程(概念版)
  17. 微信公众号开发之自动消息回复
  18. 从零开始学习CANoe(十一)—— 信号发生器(Signal Generator)
  19. 20155232 2016-2017-2《Java程序设计》课程总结
  20. ThinkCMF框架任意内容包含漏洞

热门文章

  1. pad平板性能测试软件,安兔兔评测iPad版
  2. 网络编程(3)-----------Javaweb
  3. Web基础(三)Python Web
  4. Activiti6+SpringBoot---会签功能学习-1
  5. 英伟达GPU产品分类说明
  6. 仓库管理怎么做?这5大仓库管理软件能帮到你
  7. keil 烧录时 Full Chip Erase Failed问题
  8. mvn clean package是个什么鬼?
  9. 笔记本如何选择修复计算机,详细教你笔记本电脑无线网卡修复
  10. uniapp开发小程序使用腾讯云IM(初始化配置,登录,监听,加群)