bzoj 4499: 线性函数

0没有逆元！

0没有逆元！

0没有逆元！

重要的事情说三遍！

把query的式子展开就是

首先是一个看上去很对的做法：

另d[i]=PI[j=i to n] k[i]

那么答案就是

用线段树维护d，维护b[i-1]*d[i]的区间和就可以了。最后答案*d[r+1]^(-1)

但是，k可以为0，0没有逆元，所以不能这么做。

再仔细一看，第一个式子，之间生用线段树维护不就行了？水题我却做了这么久。

   
   

    
    

     
     #include<cstdio>
    
    
    
    

     
     #include<cstdlib>
    
    
    
    

     
     #include<cstring>
    
    
    
    

     
     #include<algorithm>
    
    
    
    

     
      
    
    
    
    

     
     #define ll long long
    
    
    
    

     
     #define inf 1e9
    
    
    
    

     
     #define eps 1e-8
    
    
    
    

     
     #define md 1000000007
    
    
    
    

     
     #define N 200010
    
    
    
    

     
     using namespace std;
    
    
    
    

     
     struct Tr { ll sum,k;} tr[4*N];
    
    
    
    

     
     ll b[N],k[N];
    
    
    
    

     
      
    
    
    
    

     
     Tr update(Tr l,Tr r)
    
    
    
    

     
     {
    
    
    
    

     
     	return (Tr){(l.sum*r.k+r.sum)%md,l.k*r.k%md};
    
    
    
    

     
     }
    
    
    
    

     
      
    
    
    
    

     
     void build(int i,int l,int r)
    
    
    
    

     
     {
    
    
    
    

     
     	if (l==r) { tr[i]=(Tr){b[l],k[l]}; return;}
    
    
    
    

     
     	int mid=(l+r)>>1;
    
    
    
    

     
     	build(i<<1,l,mid); build(i<<1|1,mid+1,r);
    
    
    
    

     
     	tr[i]=update(tr[i<<1],tr[i<<1|1]);
    
    
    
    

     
     }
    
    
    
    

     
      
    
    
    
    

     
     void modify(int i,int l,int r,int x,ll k,ll b)
    
    
    
    

     
     {
    
    
    
    

     
     	if (l==r) { tr[i]=(Tr){b,k}; return;}
    
    
    
    

     
     	int mid=(l+r)>>1;
    
    
    
    

     
     	if (x<=mid) modify(i<<1,l,mid,x,k,b);
    
    
    
    

     
     	  else modify(i<<1|1,mid+1,r,x,k,b);
    
    
    
    

     
     	tr[i]=update(tr[i<<1],tr[i<<1|1]);
    
    
    
    

     
     }
    
    
    
    

     
      
    
    
    
    

     
     Tr query(int i,int l,int r,int ql,int qr)
    
    
    
    

     
     {
    
    
    
    

     
     	if (ql<=l&&r<=qr) return tr[i];
    
    
    
    

     
     	int mid=(l+r)>>1;
    
    
    
    

     
     	if (qr<=mid) return query(i<<1,l,mid,ql,qr);
    
    
    
    

     
     	if (mid+1<=ql) return query(i<<1|1,mid+1,r,ql,qr);
    
    
    
    

     
     	return update(query(i<<1,l,mid,ql,qr),query(i<<1|1,mid+1,r,ql,qr));
    
    
    
    

     
     }
    
    
    
    

     
      
    
    
    
    

     
     int main()
    
    
    
    

     
     {
    
    
    
    

     
     	int n,m; char st[5];
    
    
    
    

     
     	scanf("%d%d",&n,&m);
    
    
    
    

     
     	for (int i=1;i<=n;i++) scanf("%lld%lld",&k[i],&b[i]);
    
    
    
    

     
     	build(1,1,n);
    
    
    
    

     
     	for (int i=1;i<=m;i++)
    
    
    
    

     
     	{
    
    
    
    

     
     		scanf("%s",st);
    
    
    
    

     
     		if (st[0]=='M')
    
    
    
    

     
     		{
    
    
    
    

     
     			int x; ll K,B;
    
    
    
    

     
     			scanf("%d%lld%lld",&x,&K,&B);
    
    
    
    

     
     			modify(1,1,n,x,K,B);
    
    
    
    

     
     		}
    
    
    
    

     
     		else
    
    
    
    

     
     		{
    
    
    
    

     
     			int l,r; ll x;
    
    
    
    

     
     			scanf("%d%d%lld",&l,&r,&x);
    
    
    
    

     
     			Tr t=query(1,1,n,l,r);
    
    
    
    

     
     			printf("%lld\n",(t.k*x+t.sum)%md);
    
    
    
    

     
     		}
    
    
    
    

     
     	}
    
    
    
    

     
     	return 0;
    
    
    
    

     
     }