I guess there's not much point in reminding you that Nvodsk winters aren't exactly hot. That increased the popularity of the public transport dramatically. The route of bus 62 has exactly n stops (stop 1 goes first on its way and stop n goes last). The stops are positioned on a straight line and their coordinates are 0 = x1 < x2 < ... < xn.
Each day exactly m people use bus 62. For each person we know the number of the stop where he gets on the bus and the number of the stop where he gets off the bus. A ticket from stop a to stop b (a < b) costs xb - xa rubles. However, the conductor can choose no more than one segment NOT TO SELL a ticket for. We mean that conductor should choose C and D (С <= D) and sell a ticket for the segments [A, C] and [D, B], or not sell the ticket at all. The conductor and the passenger divide the saved money between themselves equally. The conductor's "untaxed income" is sometimes interrupted by inspections that take place as the bus drives on some segment of the route located between two consecutive stops. The inspector fines the conductor by c rubles for each passenger who doesn't have the ticket for this route's segment.
You know the coordinated of all stops xi; the numbers of stops where the i-th passenger gets on and off, ai and bi (ai < bi); the fine c; and also pi — the probability of inspection on segment between the i-th and the i + 1-th stop. The conductor asked you to help him make a plan of selling tickets that maximizes the mathematical expectation of his profit.
The first line contains three integers n, m and c (2 ≤ n ≤ 1.5·105, 1 ≤ m ≤ 3·105, 1 ≤ c ≤ 104).
The next line contains n integers xi (0 ≤ xi ≤ 109, x1 = 0, xi < xi + 1) — the coordinates of the stops on the bus's route.
The third line contains n - 1 integer pi (0 ≤ pi ≤ 100) — the probability of inspection in percents on the segment between stop i and stop i + 1.
Then follow m lines that describe the bus's passengers. Each line contains exactly two integers ai and bi (1 ≤ ai < bi ≤ n) — the numbers of stops where the i-th passenger gets on and off.
Print the single real number — the maximum expectation of the conductor's profit. The answer will be considered correct if its relative or absolute error does not exceed 10 - 6.
3 3 10 0 10 100 100 0 1 2 2 3 1 3
90.000000000
10 8 187 0 10 30 70 150 310 630 1270 2550 51100 13 87 65 0 100 44 67 3 4 1 10 2 9 3 8 1 5 6 10 2 7 4 10 4 5
76859.990000000
A comment to the first sample:
The first and third passengers get tickets from stop 1 to stop 2. The second passenger doesn't get a ticket. There always is inspection on the segment 1-2 but both passengers have the ticket for it. There never is an inspection on the segment 2-3, that's why the second passenger gets away with the cheating. Our total profit is (0 + 90 / 2 + 90 / 2) = 90.
題意:在一條直線上從左到右有N個站,已知每個點的座標X[i]。汽車從第一個站開往第N個站。有M個乘客,已知每個乘客從ai站上,在bi站下。票價爲X[bi]-X[ai]。對於每個乘客,售票員可以選擇一個區間[l,r](也可以不選),不賣站l到站r這段區間的票,剩下的前售票員和乘客平分。對於第i個站到第i+1個站,有pi的概率(pi表示百分比)會被查,被查以後每有一個人沒有這段路的票就罰售票員C。求售票員能賺到的錢的期望的最大值。
題解:
對於第i個人,售票員若不賣他區間[l,r]的票,他賺錢的期望爲:
(X[r]-X[l])/2-(pl+......pr)*C/100
處理出p數組的前綴和數組sump
則上式可以寫爲:
X[r]/2-sump[r]*c/100 - x[l]/2+sump[l-1]*C/100
另a[r]=X[r]/2-sump[r]*c/100,b[l]=- x[l]/2+sump[l-1]*C/100
c[l,r]=a[r]+b[l]且l<r
代碼如下:
#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<queue>
#include<stack>
#include<set>
#include<map>
#include<vector>
#include<string.h>
#include<string>
#include<stdlib.h>
typedef __int64 LL;
typedef unsigned __int64 LLU;
const int nn=310000;
const int inf=0x3fffffff;
const LL inf64=(LL)inf*inf;
const double eps=1e-8;
using namespace std;
LL n,m,c;
LL x[nn],p[nn];
LL mv1[nn*4],mv2[nn*4];
LL mans[nn*4];
LL sum[nn];
void update(int id)
{
mans[id]=mv1[2*id+1]+mv2[2*id];
mans[id]=max(mans[id],mans[2*id+1]);
mans[id]=max(mans[id],mans[2*id]);
mv1[id]=max(mv1[2*id],mv1[2*id+1]);
mv2[id]=max(mv2[2*id],mv2[2*id+1]);
}
void build(int id,int l,int r)
{
if(l==r)
{
mv1[id]=50*x[r]-c*sum[r-1];
mv2[id]=-50*x[l]+c*sum[l-1];
mans[id]=0;
return ;
}
int mid=(l+r)/2;
build(2*id,l,mid);
build(2*id+1,mid+1,r);
update(id);
}
vector<int>ve;
void solve(int id,int l,int r,int L,int R)
{
if(L<=l&&r<=R)
{
ve.push_back(id);
return ;
}
int mid=(l+r)/2;
if(L<=mid)
solve(2*id,l,mid,L,R);
if(R>mid)
solve(2*id+1,mid+1,r,L,R);
}
int main()
{
int i;
LL l,r;
while(scanf("%I64d%I64d%I64d",&n,&m,&c)!=EOF)
{
for(i=1;i<=n;i++)
{
scanf("%I64d",&x[i]);
}
sum[0]=0;
for(i=1;i<n;i++)
{
scanf("%I64d",&p[i]);
sum[i]=sum[i-1]+p[i];
}
build(1,1,n);
LL ans=0;
LL m1,m2,m3;
m1=m2=m3;
LL a1,a2,a3;
for(i=1;i<=m;i++)
{
scanf("%I64d%I64d",&l,&r);
ve.clear();
solve(1,1,n,l,r);
m1=mv1[ve[0]];
m2=mv2[ve[0]];
m3=mans[ve[0]];
for(int j=1;j<(int)ve.size();j++)
{
a3=max(m3,mans[ve[j]]);
a3=max(a3,mv1[ve[j]]+m2);
a1=max(m1,mv1[ve[j]]);
a2=max(m2,mv2[ve[j]]);
m1=a1,m2=a2,m3=a3;
}
ans+=m3;
}
printf("%lf\n",ans/100.0);
}
return 0;
}