Description

在JIH考察的地图中有N个城市,被公路依次连成了一个环,JIH想在这些城市中建一个玩具厂。城市和公路都被编号为1..N,i号公路连接i-1号城市与i号城市(1号公路连接N号城市与1号城市),每个城市对玩具的需求为wi,每条公路的长度为di。当JIH在第i号城市建玩具厂时,JIH需要将玩具运输到其他城市(当然i城市除外)。设第i号城市到第j号城市的两条路径长度分别为l1、l2,则将玩具运输到第j号城市的费用为l1*l2*wj。总的运输费用为将玩具运到所有城市的运输费用的总和。
JIH当然想要总的运输费用最少,所以他会选最优的城市建玩具厂,如果有多个最优的城市,小月会等概率的选取其中一个建玩具厂。
由于JIH的调查工作没做好,只知道1..N-1号城市的wi,而N号城市的wi只知道它的取值范围[a,b],假设wi的值在实数区间[a,b]上的概率是均匀分布的。现在JIH只好去进行第二次调查,于是我们想知道每个城市建玩具厂的概率是多少。
 
 

Input

第一行有三个正整数N,a,b。
接下来N-1行每行一个正实数,为w[1]到w[N-1]。
接下来N行每行一个正实数,为d[1]到d[N]。
 

Output

N行,每行一个实数,表示在第i个城市建厂的概率
 

Sample Input

5 1 100
50
25
25
50
1
2
3
2
1

Sample Output

0.090
0.000
0.000
0.090
0.821

HINT

当w[5]<18.75时,将在1或4号城市建玩具厂,当w[5]>18.75时,将在5号城市建玩具厂,
当w[5]=18.75时,将在1或4或5建玩具厂。
100%的数据中:N<=100000,a<=b<=10000,w[i]<=10000,d[i]<=10。
 
 
先将每个点的代价表示为kx+b的形式,然后维护个凸包求出每段区间最低的直线,进而求出每条直线即每个点被选的概率。
现在瓶颈在求出每个点的b。
如果将厂设在i-1点 l[j],r[j]为i-1点和j点的两条路长度,则
b[i-1]=Σl[j]*r[j]*w[j]
同理b[i]=Σ(l[j]+d[i])*(r[j]-d[i])*w[j]
作差得 b[i]-b[i-1]=Σw[j]*(r[j]*d[i]-l[j]*d[i]-d[i]*d[i])
=Σw[j]*r[j]*d[i]-Σw[j]*l[j]*d[i]-Σw[j]*d[i]*d[i]
维护Σw[j]*r[j]Σw[j]*l[j]即可。
先暴力求出第1个点的b,然后O(1)转移求其他点的b,求所有点的b变成了O(N)。
注意:有重复的直线,他们要平分概率。
code:
 #include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm>
#define maxn 100005
#define eps 1E-9
#define inf 1E100
using namespace std;
int n,flag,cnt,siz,pos[maxn];
double w[maxn],d[maxn],dist[maxn],tot,tmp,ans,t1,t2,t3,val[maxn];
double a,b,l,r,f[maxn],p[maxn];
struct Line{
double k,b;
int id;
}line[maxn],tl[maxn];
bool cmp(Line x,Line y){
if (abs(x.k-y.k)>eps) return x.k>y.k;
return x.b<y.b;
}
double calc(Line l1,Line l2){return -(l1.b-l2.b)/(l1.k-l2.k);}
struct DATA{
double x;
int id;
}list[maxn];
void solve(){
memcpy(tl,line,sizeof(line));
sort(tl+,tl+n+,cmp);
for (int i=,j,t;i<=n;i=j){
line[++cnt]=tl[i],t=;
for (j=i;tl[j].k==tl[i].k;j++) if (tl[j].b==tl[i].b) pos[tl[j].id]=cnt,t++;
for (j=i;tl[j].k==tl[i].k;j++) if (tl[j].b==tl[i].b) p[tl[j].id]=1.0/t;
}
list[siz=]=(DATA){-inf,};
for (int i=;i<=cnt;i++){
double t=calc(line[i],line[list[siz].id]);
while (t<list[siz].x) siz--,t=calc(line[i],line[list[siz].id]);
list[++siz]=(DATA){t,i};
}
list[++siz]=(DATA){inf,cnt+};
l=a,r=b,flag=;
for (int i=;i<=siz&&flag!=;i++){
if (list[i].x>l){
if (!flag) flag=,f[list[i-].id]=(list[i].x-l)/(r-l);
else if (flag==&&list[i].x<r) f[list[i-].id]=(list[i].x-list[i-].x)/(r-l);
else if (list[i].x>=r) f[list[i-].id]=(r-list[i-].x)/(r-l),flag=;
}
}
for (int i=;i<=n;i++) p[i]=f[pos[i]]*p[i];
for (int i=;i<=n;i++) printf("%.3lf\n",p[i]);
}
int main(){
scanf("%d%lf%lf",&n,&a,&b);
for (int i=;i<n;i++) scanf("%lf",&w[i]);
for (int i=;i<=n;i++) scanf("%lf",&d[i]);
for (int i=;i<=n;i++) dist[i]=d[i];
for (int i=;i<=n;i++) dist[i]+=dist[i-];
tot=dist[n];
if (a==b){
w[n]=a;
t1=w[]*tot,t3=w[];
for (int i=;i<=n;i++){
tmp+=w[i]*(dist[i]-dist[])*(tot-dist[i]+dist[]);
t1+=w[i]*(dist[i]-dist[]);
t2+=w[i]*(tot-dist[i]+dist[]);
t3+=w[i];
}
val[]=ans=tmp;
for (int i=;i<=n;i++){
tmp+=t1*d[i]-t2*d[i]-t3*d[i]*d[i];
ans=min(ans,tmp);
val[i]=tmp;
t1+=w[i]*tot,t2-=w[i]*tot;
t1-=t3*d[i],t2+=t3*d[i];
}
cnt=;
for (int i=;i<=n;i++) if (abs(val[i]-ans)<=eps) cnt++;
for (int i=;i<=n;i++) if (abs(val[i]-ans)<=eps) printf("%.3lf\n",1.0/cnt);
else puts("0.000");
return ;
}
l=,r=tot;
for (int i=;i<=n;i++) l+=d[i],r-=d[i],line[i].k=l*r;
t1=w[]*tot,t3=w[];
for (int i=;i<=n;i++){
tmp+=w[i]*(dist[i]-dist[])*(tot-dist[i]+dist[]);
t1+=w[i]*(dist[i]-dist[]);
t2+=w[i]*(tot-dist[i]+dist[]);
t3+=w[i];
}
line[].b=tmp;
for (int i=;i<=n;i++){
tmp+=t1*d[i]-t2*d[i]-t3*d[i]*d[i];
line[i].b=tmp;
t1+=w[i]*tot,t2-=w[i]*tot;
t1-=t3*d[i],t2+=t3*d[i];
}
for (int i=;i<=n;i++) line[i].id=i;
solve();
return ;
}

bzoj4306: 玩具厂的更多相关文章

  1. C#设计模式之装饰者

    IronMan之装饰者 前言 上一篇的文章我们讲到要给"IronMan"配备"武器",并且还使用了"武器",效果还是不错的,对于多种环境.多 ...

  2. C#设计模式之工厂

    IronMan之工厂 前言 实用为主,学一些用得到的技术更能在如今的社会里保命. 虽然在日常的工作中设计模式不是经常的用到,但是呢,学习它只有好处没有坏处. 设计模式像是一种“标签”,它是代码编写者思 ...

  3. 2017.07.11【NOIP提高组】模拟赛B组

    Summary 今天的比赛打得还不错,第一题被同桌灌输的贪心,纯模拟洗脑了,然后steal的看了一下,发现怎么也对不了,一直在检查.最后10分钟才找出反例,推出动态规划方程,没有想到怎么转移,比赛就结 ...

  4. python 全栈开发,Day132(玩具管理页面,控制玩具通讯录,基于请求的好友关系建立)

    先下载github代码,下面的操作,都是基于这个版本来的! https://github.com/987334176/Intelligent_toy/archive/v1.5.zip 注意:由于涉及到 ...

  5. python 全栈开发,Day129(玩具开机提示语,为多个玩具发送点播,聊天界面,app录音,app与服务器端文件传输,简单的对话)

    一.玩具开机提示语 先下载github代码,下面的操作,都是基于这个版本来的! https://github.com/987334176/Intelligent_toy/archive/v1.2.zi ...

  6. 45.更新一下scrapy爬取工商信息爬虫代码

    这里是完整的工商信息采集代码,不过此程序需要配合代理ip软件使用.问题:1.网站对ip之前没做限制,但是采集了一段时间就被检测到设置了反爬,每个ip只能访问十多次左右就被限制访问.2.网站对请求头的检 ...

  7. VR与AR的发展趋势分析

    概要 你是否想象过与神秘的深海生物近距离接触?你是否梦想过穿戴钢铁侠那样的超先进科技装备成为超级英雄?你又是否幻想过与梦中的女神面对面的交流?这些可能在以前都只能是存在于脑海中的幻想,可是在如今有一项 ...

  8. 《深入浅出数据分析》-利用Excel的Solver求橡皮玩具的最大利润

    背景:一玩具厂可以生产两种浴盆玩具,分别是橡皮鸭和橡皮鱼,并且原材料和生产时间都有所限制,求如何才分配生产两种玩具的数量才可以让厂商达到最大利润. 假设条件如下: 产品名称 数量   duck a   ...

  9. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

随机推荐

  1. 《Mathematical Olympiad——数论》——整除

    数论这个东西吧,虽说也是高中IMOer玩的数学游戏,颇具美学性的证明比较多.就目前所知,它在算法里是一些加密技术的基础,不多言,开始具体题目的分析. 问题一:已知数列{an},且a0 = 2 , a1 ...

  2. errno与perror

    很多系统函数在错误返回时将错误原因记录在libc定义的全局变量errno中,每种错误原因对应一个错误码,请查阅errno(3)的Man Page了解各种错误码,errno在头文件errno.h中声明, ...

  3. Java Web应用的开发模式

    一.概述 从互联网诞生的第一天起,web从一个简单的静态页面,发展到今天五花八门.令人眼花缭乱的复杂应用,大到企业级的web应用系统,小到简单的SPA管理系统.在硬件性能提升的同时,通过各种技术实现了 ...

  4. Python读写文件模式

    1.r 打开只读文件,该文件必须存在. 2.r+ 打开可读写的文件,该文件必须存在. 3.w 打开只写文件,若文件存在则文件长度清为0,即该文件内容会消失.若文件不存在则建立该文件. 4.w+ 打开可 ...

  5. C++发送邮件和附件

      c++socketnulldelete服务器stream 头文件 /**************************************************************** ...

  6. linux查看cpu温度

      分类: linux系统 一.安装  sudo apt-get install lm-sensors   二.查看 linux@cdyemail:~$ sensors k10temp-pci-00c ...

  7. 深入探索C++对象模型-5

    虚拟继承下的对象构造: 由于虚拟基类对象在子类中只能保持一个实例,那么,子类构造的时候调用父类的构造函数的时候必须保证虚拟基类对象不能够重复构造. 那么如何保证基类对象的唯一性? C++规定虚拟基类对 ...

  8. 在Ubuntu上安装使用Systemtap

    因为最近开始学习Nginx,在网上看到别人介绍了一款强大的内核探测工具Systemtap,于是便准备学习下这款探测工具为以后代码分析做准备. 第一步便是安装.在自己电脑上安装的时候,也是费了一番劲儿. ...

  9. Codeforces Round #291 (Div. 2) C - Watto and Mechanism 字符串

    [题意]给n个字符串组成的集合,然后有m个询问(0 ≤ n ≤ 3·105, 0 ≤ m ≤ 3·105) ,每个询问都给出一个字符串s,问集合中是否存在一个字符串t,使得s和t长度相同,并且仅有一个 ...

  10. <c:if>判断参数是否同时为空

    <c:if test="${empty str}">  str为空</c:if> <c:if test="${not empty str}& ...