[HAOI2012]高速公路

Posted cjfdf

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[HAOI2012]高速公路相关的知识,希望对你有一定的参考价值。

题面在这里

题意

维护区间加操作+询问区间任选两不同点途中线段权值之和的期望

sol

一道假的期望题...
因为所有事件的发生概率都相同,所以答案就是所有方案的权值总和/总方案数
因为区间加法自然想到线段树,考虑将每条道路看作一个节点
那么对于区间l-r的路径(注意这里已经将道路看作节点,所以r--),可以计算第x条道路的贡献为
\((i-l+1)*(r-i+1)*val[i]\)(\(val[i]\)表示当前道路的权值大小)
那么我们要区间维护的就是这个东西
\[\sum_{r}^{i=l}{(i-l+1)*(r-i+1)*val[i]}\]
\[=\sum_{r}^{i=l}{(-i^2*va[i])}+(l+r)*\sum_{r}^{i=l}{(i*val[i])+(-l*r-l+r+1)*\sum_{r}^{i=l}val[i]}\]
分别维护\(\sum_{r}^{i=l}(i^2*val[i])\),\(\sum_{r}^{i=l}{(i*val[i])}\)\(\sum_{r}^{i=l}{val[i]}\),最后再除以\(((r-l+1)*(r-l+2)/2)\)即可

代码

#include<bits/stdc++.h>
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<iomanip>
#include<cstring>
#include<vector>
#include<cstdio>
#include<string>
#include<bitset>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<set>
#define mp make_pair
#define pb push_back
#define RG register
#define il inline
using namespace std;
const int mod=1e9+7;
const int N=100010;
const double eps=1e-10;
typedef unsigned long long ull;
typedef vector<int>VI;
typedef long long ll;
typedef double dd;
il ll read(){
    RG ll data=0,w=1;RG char ch=getchar();
    while(ch!=‘-‘&&(ch<‘0‘||ch>‘9‘))ch=getchar();
    if(ch==‘-‘)w=-1,ch=getchar();
    while(ch<=‘9‘&&ch>=‘0‘)data=data*10+ch-48,ch=getchar();
    return data*w;
}

int n,m;char ch;
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
#define ls (i<<1)
#define rs (i<<1|1)
#define mid ((l+r)>>1)
ll x[N],xx[N];
ll sum[N<<3],sumx[N<<3],sumxx[N<<3],lz[N<<3],cnt;
il void update(int i){
    sum[i]=sum[ls]+sum[rs];
    sumx[i]=sumx[ls]+sumx[rs];
    sumxx[i]=sumxx[ls]+sumxx[rs];
}

il void change(int i,int l,int r,ll v){
    sum[i]+=v*(r-(l-1));sumx[i]+=v*(x[r]-x[l-1]);
    sumxx[i]+=v*(xx[r]-xx[l-1]);lz[i]+=v;
}

il void pushdown(int i,int l,int r){
    if(!lz[i])return;
    change(ls,l,mid,lz[i]);change(rs,mid+1,r,lz[i]);
    lz[i]=0;
}

il void modify(int i,int l,int r,int x,int y,ll v){
    //区间修改
    if(x<=l&&r<=y){change(i,l,r,v);return;}
    pushdown(i,l,r);
    if(x<=mid)modify(ls,l,mid,x,y,v);
    if(mid+1<=y)modify(rs,mid+1,r,x,y,v);
    update(i);
}

il ll query_sum(int i,int l,int r,int x,int y){
    if(x<=l&&r<=y)return sum[i];
    RG ll s=0;pushdown(i,l,r);
    if(x<=mid)s+=query_sum(ls,l,mid,x,y);
    if(y>=mid+1)s+=query_sum(rs,mid+1,r,x,y);
    return s;
}
il ll query_sumx(int i,int l,int r,int x,int y){
    if(x<=l&&r<=y)return sumx[i];
    RG ll s=0;pushdown(i,l,r);
    if(x<=mid)s+=query_sumx(ls,l,mid,x,y);
    if(y>=mid+1)s+=query_sumx(rs,mid+1,r,x,y);
    return s;
}
il ll query_sumxx(int i,int l,int r,int x,int y){
    if(x<=l&&r<=y)return sumxx[i];
    RG ll s=0;pushdown(i,l,r);
    if(x<=mid)s+=query_sumxx(ls,l,mid,x,y);
    if(y>=mid+1)s+=query_sumxx(rs,mid+1,r,x,y);
    return s;
}

il void query(int l,int r){
    //按照推导公式解决询问操作
    RG ll ans1=query_sumxx(1,1,n-1,l,r)*(-1)+query_sumx(1,1,n-1,l,r)*(l+r)+query_sum(1,1,n-1,l,r)*(1ll*l*r*(-1)-l+r+1);
    RG ll ans2=1ll*(r-l+1)*(r-l+2)/2;
    RG ll g=gcd(ans1,ans2);
    printf("%lld/%lld\n",ans1/g,ans2/g);
}

int main()
{
    n=read();m=read();
    for(RG int i=1;i<n;i++){
        x[i]=x[i-1]+i;xx[i]=xx[i-1]+1ll*i*i;
    }
    for(RG int i=1,l,r;i<=m;i++){
        ch=0;while(ch!=‘C‘&&ch!=‘Q‘)ch=getchar();
        if(ch==‘C‘){
            l=read();r=read()-1;RG ll v=read();
            modify(1,1,n-1,l,r,v);
        }
        else{l=read();r=read()-1;query(l,r);}
    }
    return 0;
}

以上是关于[HAOI2012]高速公路的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ2752: [HAOI2012]高速公路(road)

「Luogu2221」[HAOI2012]高速公路

[HAOI2012]高速公路

BZOJ 2752:[HAOI2012]高速公路(road)(线段树)

BZOJ 2752: [HAOI2012]高速公路(road) [线段树 期望]

BZOJ2752[HAOI2012]高速公路(road) 线段树