Description

Solution

有两种做法

一种是线段树维护一次方程系数,一种是线段树维护矩阵

准备都写一写


维护系数

首先把式子推出来 $$CS=B\times \sum\limits_{i=1}^np_i\times(f_{i-1}+1)$$

\[f_i=(t+p_i-p_i\times t)\times f_{i-1}+p_i
\]

发现 \(f_i\) 是关于 \(f_{i-1}\) 的一次函数 \(y=kx+b\) 形式,可以线段树维护这个 \(k\) 和 \(b\)。

还有 \(p_i\times(f_{i-1}+1)=p_i\times f_{i-1}+p_i\) 也是个一次函数,也能用线段树维护。

具体都维护些什么呢?

在区间 \([L,R]\) 维护 \(k_1,b_1\) 表示 \(f_R=k_1f_{L-1}+b_1\) 和 \(k_2,b_2\) 表示 \(CS[L,R]=k_2f_{L-1}+b_2\)

观察到这样做的好处就是右区间的 \(L-1\) 等于左区间的 \(R\),这样就资瓷快速合并了。

然后如果一个询问是 \([l,r]\),那求出来的 \(CS[L,R]\) 的 \(b_2\) 就是答案了。

#include<set>
#include<map>
#include<cmath>
#include<queue>
#include<cctype>
#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using std::min;
using std::max;
using std::swap;
using std::vector;
typedef double db;
const int N=500005;
typedef long long ll;
const int mod=998244353;
#define pb(A) push_back(A)
#define pii std::pair<int,int>
#define mp(A,B) std::make_pair(A,B)
#define ls cur<<1
#define rs cur<<1|1
#define lss ls,l,mid,ql,qr
#define rss rs,mid+1,r,ql,qr int n,q,t,tb,A,B,b1[N<<2],b2[N<<2];
int p[N],sump[N<<2],k1[N<<2],k2[N<<2]; struct Node{
int sump,k1,b1,k2,b2;
friend Node operator+(Node x,Node y){
Node z;
z.sump=(x.sump+y.sump)%mod;
z.k1=(ll)y.k1*x.k1%mod;
z.b1=((ll)y.k1*x.b1%mod+y.b1)%mod;
z.k2=((ll)y.k2*x.k1%mod+x.k2)%mod;
z.b2=((ll)y.k2*x.b1%mod+y.b2+x.b2)%mod;
return z;
}
}; int getint(){
int X=0,w=0;char ch=0;
while(!isdigit(ch))w|=ch=='-',ch=getchar();
while( isdigit(ch))X=X*10+ch-48,ch=getchar();
if(w) return -X;return X;
} int ksm(int x,int y){
int ans=1;
while(y){
if(y&1) ans=(ll)ans*x%mod;
x=(ll)x*x%mod;y>>=1;
} return ans;
} int inv(int x){
return ksm(x,mod-2);
} void pushup(int cur){
k1[cur]=(ll)k1[rs]*k1[ls]%mod;
b1[cur]=((ll)k1[rs]*b1[ls]%mod+b1[rs])%mod;
k2[cur]=((ll)k2[rs]*k1[ls]%mod+k2[ls])%mod;
b2[cur]=((ll)k2[rs]*b1[ls]%mod+(ll)b2[rs]+b2[ls])%mod;
sump[cur]=((ll)sump[ls]+sump[rs])%mod;
} void build(int cur,int l,int r){
if(l==r){
k1[cur]=((ll)t+p[l]-(ll)p[l]*t%mod+mod)%mod;
b1[cur]=p[l];
k2[cur]=p[l];
b2[cur]=p[l];
sump[cur]=p[l];
return;
} int mid=l+r>>1;
build(ls,l,mid);build(rs,mid+1,r);
pushup(cur);
} void modify(int cur,int l,int r,int ql,int qr,int c){
if(l==r){
k1[cur]=((ll)t+c-(ll)c*t%mod+mod)%mod;
b1[cur]=k2[cur]=b2[cur]=sump[cur]=c;
return;
} int mid=l+r>>1;
ql<=mid?modify(lss,c):modify(rss,c);
pushup(cur);
} Node query(int cur,int l,int r,int ql,int qr){
if(ql<=l and r<=qr) return (Node){sump[cur],k1[cur],b1[cur],k2[cur],b2[cur]};
int mid=l+r>>1;
if(qr<=mid) return query(lss);
if(ql>mid) return query(rss);
return query(lss)+query(rss);
} signed main(){
freopen("omeed.in","r",stdin);freopen("omeed.out","w",stdout);
getint();n=getint(),q=getint(),t=getint(),tb=getint(),A=getint(),B=getint();
t=(ll)t*inv(tb)%mod;
for(int i=1;i<=n;i++){
int x=getint(),y=getint();
p[i]=(ll)x*inv(y)%mod;
} build(1,1,n);
while(q--){
if(getint()==0){
int pos=getint(),a=getint(),b=getint(),c=(ll)a*inv(b)%mod;
p[pos]=c;modify(1,1,n,pos,pos,c);
} else{
int l=getint(),r=getint();
Node ans=query(1,1,n,l,r);
printf("%d\n",((ll)ans.sump*A%mod+(ll)ans.b2*B%mod)%mod);
}
} return 0;
}

维护矩阵

上面的DP式子已经列出来了。

我们现在要解决的这样一类问题:我们知道DP式子,现在要求多次从不同的起点开始做DP

求出DP的转移矩阵然后线段树上动态DP就好了

姿势不对会T一个点懒得卡常了

#include<set>
#include<map>
#include<cmath>
#include<queue>
#include<cctype>
#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using std::min;
using std::max;
using std::swap;
using std::vector;
typedef double db;
const int N=500005;
const int mod=998244353;
typedef long long ll;
#define pb(A) push_back(A)
#define pii std::pair<Mat,int>
#define mp(A,B) std::make_pair(A,B)
#define ls cur<<1
#define rs cur<<1|1
#define lss ls,l,mid,ql,qr
#define rss rs,mid+1,r,ql,qr int n,q,t,ta,A,B;
int p[N],sum[N<<2]; struct Mat{
int a[4][4]; void clear(){
memset(a,0,sizeof a);
} void init(){
clear();
for(int i=1;i<=3;i++) a[i][i]=1;
} friend Mat operator*(Mat x,Mat y){
Mat z;z.clear();
for(int i=1;i<=3;i++)
for(int k=1;k<=3;k++)
for(int j=1;j<=3;j++)
z.a[i][j]=((ll)z.a[i][j]+(ll)x.a[i][k]*y.a[k][j]%mod)%mod;
return z;
}
}cs[N<<2]; void pushup(int cur){
cs[cur]=cs[ls]*cs[rs];
sum[cur]=(sum[ls]+sum[rs])%mod;
} void build(int cur,int l,int r){
if(l==r){
cs[cur].a[1][1]=((ll)p[l]+t-(ll)p[l]*t%mod+mod)%mod;cs[cur].a[1][2]=(ll)B*p[l]%mod;cs[cur].a[1][3]=0;
cs[cur].a[2][1]=0;cs[cur].a[2][2]=1;cs[cur].a[2][3]=0;sum[cur]=p[l];
cs[cur].a[3][1]=p[l];cs[cur].a[3][2]=(ll)p[l]*(B+A)%mod;cs[cur].a[3][3]=1;
return;
} int mid=l+r>>1;
build(ls,l,mid);build(rs,mid+1,r);pushup(cur);
} void modify(int cur,int l,int r,int ql,int qr,int c){
if(l==r){
cs[cur].a[1][1]=((ll)c+t-(ll)c*t%mod+mod)%mod;cs[cur].a[1][2]=(ll)B*c%mod;cs[cur].a[1][3]=0;
cs[cur].a[2][1]=0;cs[cur].a[2][2]=1;cs[cur].a[2][3]=0;sum[cur]=c;
cs[cur].a[3][1]=c;cs[cur].a[3][2]=(ll)c*(B+A)%mod;cs[cur].a[3][3]=1;return;
} int mid=l+r>>1;
ql<=mid?modify(lss,c):modify(rss,c);
pushup(cur);
} Mat query(int cur,int l,int r,int ql,int qr){
if(ql<=l and r<=qr) return cs[cur];
int mid=l+r>>1;
if(qr<=mid) return query(lss);
if(ql>mid) return query(rss);
return query(lss)*query(rss);
} int getint(){
int X=0,w=0;char ch=0;
while(!isdigit(ch))w|=ch=='-',ch=getchar();
while( isdigit(ch))X=X*10+ch-48,ch=getchar();
if(w) return -X;return X;
} int ksm(int a,int y){
int ans=1;
while(y){
if(y&1) ans=(ll)ans*a%mod;
a=(ll)a*a%mod;y>>=1;
} return ans;
} int inv(int x){
return ksm(x,mod-2);
} signed main(){
getint(),n=getint(),q=getint(),t=getint(),ta=getint(),A=getint(),B=getint();
t=(ll)t*inv(ta)%mod;
for(int i=1;i<=n;i++){
int x=getint(),y=getint();
p[i]=(ll)x*inv(y)%mod;
} build(1,1,n);
while(q--){
if(getint()==0){
int pos=getint(),x=getint(),y=getint(),z=(ll)x*inv(y)%mod;
modify(1,1,n,pos,pos,z);
} else{
int l=getint(),r=getint();
printf("%d\n",query(1,1,n,l,r).a[3][2]);
}
} return 0;
}

[JZOJ5837] Omeed的更多相关文章

  1. [JZOJ] 5837.Omeed

    先摆出来这个式子 \[ score=A\sum S_i+B\sum S_i\times f(i) \] 先研究\(f\)函数(也就是Combo函数) 显然的有 \[ f(i)=P_i(f(i-1)+1 ...

  2. 20210824 Prime,Sequence,Omeed

    考场 T1 貌似是 luogu 上原题 T2 计数,想起了这题和这题,但没有 \(n^2\) 一档的分...准备打个表 T3 期望 DP,但暴力是 \(O(qn)\) 的,发现 \(combo\) 的 ...

  3. Omeed 线段树

    目录 题面 题解 代码 题面 2.12 - - - 题解 大概还是挺妙的? 首先基础分和连击分互不干扰,所以可以分开统计. 基础分的统计比较简单,等于: \[A \sum_{i = l}^{r} p_ ...

  4. 题解 Omeed

    传送门 差了一点没想到正解-- 首先单次询问的 \(O(n)\) 写法很好想,考虑如何优化 首先基础分区间求和即可 然后那个连击分的话,是一个关于 \(f_i\) 和 \(f_{i-1}\) 的柿子 ...

  5. [考试总结]noip模拟47

    感觉自己放弃题目还是过于容易. 其实第一题不是很难,但是自己拿了一个暴力就走人了.. 然后其实简单优化一下子就有不少分数. 然后第二题的本质不同的子序列个数的方程没有推出来,如果推出来就会直接有 \( ...

  6. Noip模拟47 2021.8.25

    期望得分:55+24+53 实际得分:0+0+3 乐死 累加变量清零了吗? 打出更高的部分分暴力删了吗? 样例解释换行你看见了吗? T1 Prime 打出55分做法没删原来的暴力,结果就轻松挂55分 ...

  7. 2021.8.24考试总结[NOIP47]

    T1 prime 发现只需筛小于等于$mid(\sqrt r,k)$的质数,之后用这些质数筛掉区间内不合法的数即可. $code:$ 1 #include<bits/stdc++.h> 2 ...

随机推荐

  1. 2018-10-19,下午4点拿到京东offer

    今天12点没有问题就签.岗位,地点都很满足.京东居然对我这么好,保证以后不做损坏京东利益的事写给以后的自己.

  2. Oracle中为表设置自动增长的标识列

    -- Create sequence create sequence INNERID minvalue 1 maxvalue 99999999999999 start with 1 increment ...

  3. [c#.net]未能加载文件或程序集“”或它的某一个依赖项。系统找不到指定的文件

    问题是这样嘀: 项目采用了三层架构和工厂模式,并借鉴了PetShop的架构,因为这个项目也是采用分布式的数据库,目前只有三个数据库,主要出于提高访问性能考虑. 原来是按照网上对PetShop的介绍来给 ...

  4. boost--内存池

    boost的内存池实现了一个快速.紧凑的内存分配和管理器,使用它可以完全不用考虑delete释放问题.常用的boost内存池有pool.object_pool.singleton_pool. 1.po ...

  5. Spring通过在META-INF/spring.handlers中的属性进行配置文件解析

    在Spring的入口函数refresh()之中进行的. AbstractApplicationContext ConfigurableListableBeanFactory beanFactory = ...

  6. 理解React组件的生命周期

    本文作者写作的时间较早,所以里面会出现很多的旧版ES5的时代的方法.不过,虽然如此并不影响读者理解组件的生命周期.反而是作者分为几种不同的触发机制来解释生命周期的各个方法,让读者更加容易理解涉及到的概 ...

  7. TabLoaout简单框架

    import android.os.Bundle; import android.support.annotation.Nullable; import android.support.design. ...

  8. 第十节:详细讲解一下Java多线程,随机文件

    前言 大家好,给大家带来第十节:详细讲解一下Java多线程,随机文件的概述,希望你们喜欢 多线程的概念 线程的生命周期 多线程程序的设计 多线程的概念 多线程的概念:程序是静态的,进程是动态的.多进程 ...

  9. apache环境之困扰,Rewrite导致无法加载多个不同的.html文件

    又是一个项目,为访问多个纯静态html页面h5游戏页,能够做一些简单分享和跳转即可.原本是一个简单得不能的项目,但是却多生了事端. 我按照apache的惯例,将文件上传到服务器的DocumentRoo ...

  10. jdk-tomcat-jenkens 安装

    1--安装JDK 下载JDK放到你需要的目录,解压,然后添加环境变量 2--安装tomcat 从官方网站下载tomcat的安装包,然后解压 启动tomcat , TOMCAT的默认端口是8080,要记 ...