題目地址

/*
* @Author: hannibal
* @Date:   2018-08-07 10:42:26
* @Last Modified by:   hannibal
* @Last Modified time: 2018-08-07 17:08:44
*/
#pragma GCC optimize(2)
#include<bits/stdc++.h>
using namespace std;
typedef unsigned long long  ll;
const ll mod = 10007;
const ll maxn = 1e5+10;
ll Case = 1, n, m;
struct
 
 node{
	ll l, r;
	ll add, mul, set;
	ll p1, p2, p3;
	ll mid(){return (l+r)/2;}
}tr[maxn<<2];
void pushup(ll rt) {
	tr[rt].p1 = (tr[rt<<1].p1 + tr[rt<<1|1].p1)%mod;
	tr[rt].p2 = (tr[rt<<1].p2 + tr[rt<<1|1].p2)%mod;
	tr[rt].p3 = (tr[rt<<1].p3 + tr[rt<<1|1].p3)%
 
mod;
}
void caladd(ll rt) {
	ll len1 = (tr[rt<<1].r-tr[rt<<1].l+1)%mod;
	ll len2 =( tr[rt<<1|1].r-tr[rt<<1|1].l+1)%mod;
	tr[rt<<1].add = (tr[rt<<1].add+tr[rt].add%mod);tr[rt<<1|1].add = (tr[rt<<1|1].add+tr[rt].add)%mod;
	ll temp = (tr[rt].add*tr[rt].add%
 
mod)*tr[rt].add%mod;
	tr[rt<<1].p3 =(tr[rt<<1].p3+(temp*len1)%mod+3*tr[rt].add*((tr[rt<<1].p2+tr[rt<<1].p1*tr[rt].add)%mod)%mod)%mod;
	tr[rt<<1|1].p3 = (tr[rt<<1|1].p3+(temp*len2)%mod+3*tr[rt].add*((tr[rt<<1|1].p2+tr[rt<<1|1].p1*tr[rt].add)%mod)%mod)%mod;
	tr[rt<<1].p2 = (tr[rt<<1].p2+(tr[rt].add*tr[rt].add%mod*len1%mod)+2*tr[rt<<1].p1*tr[rt].add%mod)%mod;
	tr[rt<<1|1].p2 = (tr[rt<<1|1].p2+(tr[rt].add*tr[rt].add%mod*len2%mod)+2*tr[rt<<1|1].p1*tr[rt].add%mod)%mod;
	tr[rt<<1].p1 = (tr[rt<<1].p1+len1*tr[rt].add%mod)%mod;
	tr[rt<<1|1].p1 = (tr[rt<<1|1].p1+len2*tr[rt].add%mod)%mod;
	tr[rt].add = 0;
}
void calmul(ll rt) {
	ll x = tr[rt].mul;
	tr[rt<<1].mul = (x*tr[rt<<1].mul)%mod;
	tr[rt<<1|1].mul = (x*tr[rt<<1|1].mul)%mod;
	if(tr[rt<<1].add) tr[rt<<1].add = (tr[rt<<1].add*x)%mod;
	if(tr[rt<<1|1].add) tr[rt<<1|1].add = (tr[rt<<1|1].add*x)%mod;
	ll temp = (x*x)%mod*x%mod;
	tr[rt<<1].p1 = (tr[rt<<1].p1*x)%mod;
	tr[rt<<1|1].p1 = (tr[rt<<1|1].p1*x)%mod;
	tr[rt<<1].p2 = (tr[rt<<1].p2*x%mod*x)%mod;
	tr[rt<<1|1].p2 = (tr[rt<<1|1].p2*x%mod*x)%mod;
	tr[rt<<1].p3 = (tr[rt<<1].p3*temp%mod)%mod;
	tr[rt<<1|1].p3 = (tr[rt<<1|1].p3*temp%mod)%mod;
	tr[rt].mul = 1;
}
void calset(ll rt) {
	ll len1 = (tr[rt<<1].r-tr[rt<<1].l+1)%mod;
	ll len2 = (tr[rt<<1|1].r-tr[rt<<1|1].l+1)%mod;
	tr[rt<<1].set = tr[rt<<1|1].set = tr[rt].set;
	tr[rt<<1].add = tr[rt<<1|1].add = 0;
	tr[rt<<1].mul = tr[rt<<1|1].mul = 1;
	ll temp = (tr[rt].set*tr[rt].set)%mod*tr[rt].set%mod;
	tr[rt<<1].p1 = len1*(tr[rt].set%mod)%mod;
	tr[rt<<1|1].p1 = len2*(tr[rt].set%mod)%mod;
	tr[rt<<1].p2 = len1*(tr[rt].set*tr[rt].set%mod)%mod;
	tr[rt<<1|1].p2 = len2*(tr[rt].set*tr[rt].set%mod)%mod;
	tr[rt<<1].p3 = len1*temp%mod;
	tr[rt<<1|1].p3 = len2*temp%mod;
	tr[rt].set = 0;
}
void pushdown(ll rt) {
	if(tr[rt].set) calset(rt);
	if(tr[rt].mul != 1) calmul(rt);
	if(tr[rt].add) caladd(rt);
}
void build(ll rt, ll l, ll r){
	tr[rt].l = l; tr[rt].r = r;tr[rt].mul = 1;
	tr[rt].add =  tr[rt].set = 0;
	if(l == r) {tr[rt].p1 = tr[rt].p2 = tr[rt].p3 = 0;return;}
	ll mid = tr[rt].mid();
	build(rt<<1, l, mid);build(rt<<1|1, mid+1, r);
	pushup(rt);
}
void update(ll rt, ll l, ll r, ll flag, ll c) {
	if(tr[rt].l == l && tr[rt].r == r) {
		if(flag == 1){
			tr[rt].add += c;
			ll temp = c*c%mod*c%mod*(r-l+1)%mod;
			tr[rt].p3 = (tr[rt].p3+temp+3*c*((tr[rt].p2+tr[rt].p1*c)%mod))%mod;
			tr[rt].p2 = (tr[rt].p2+(c*c%mod*(r-l+1)%mod)+2*tr[rt].p1*c%mod)%mod;
			tr[rt].p1 = (tr[rt].p1+(r-l+1)*c%mod)%mod;
		}
		else if(flag == 2) {
			tr[rt].mul = (tr[rt].mul*c)%mod;
			if(tr[rt].add) tr[rt].add = (tr[rt].add*c)%mod;
			tr[rt].p1 = (tr[rt].p1*c)%mod;
			tr[rt].p2 = (tr[rt].p2*c%mod*c)%mod;
			tr[rt].p3 = (tr[rt].p3*c%mod*c%mod*c)%mod;
		}
		else if(flag == 3) {
			tr[rt].set = c; tr[rt].add = 0; tr[rt].mul = 1;
			tr[rt].p1 = (r-l+1)%mod*c%mod;
			tr[rt].p2 = (r-l+1)%mod*c%mod*c%mod;
			tr[rt].p3 = (r-l+1)%mod*c%mod*c%mod*c%mod;
		}
		return;
	}
	pushdown(rt);
	ll mid = tr[rt].mid();
	if(r <= mid) update(rt<<1, l, r, flag, c);
	else if(l > mid) update(rt<<1|1, l, r, flag, c);
	else update(rt<<1, l, mid, flag, c), update(rt<<1|1, mid+1, r, flag, c);
	pushup(rt);
}
ll query(ll rt, ll l, ll r, ll p) {
	if(tr[rt].l == l && tr[rt].r == r) {
		if(p == 1) return tr[rt].p1%mod;
		else if(p == 2) return tr[rt].p2%mod;
		else return tr[rt].p3%mod;
	}
	pushdown(rt);
	ll mid = tr[rt].mid();
	if(r <= mid) return query(rt<<1, l, r, p);
	else if(l > mid) return query(rt<<1|1, l, r, p);
	else return (query(rt<<1, l, mid, p)+query(rt<<1|1,mid+1, r, p)+mod)%mod;
}
void solve() {
	build(1, 1, n);
	for(ll i = 1; i <= m; i++){
		ll q, x, y, p;
		scanf("%llu%llu%llu%llu", &q, &x, &y, &p);
		if(q == 4) printf("%llu\n", query(1, x, y, p)%mod);
		else update(1, x, y, q, p);
	}
}
int main() {
    ios::sync_with_stdio(false);
#ifndef ONLINE_JUDGE
    freopen("in.txt", "r",  
 

            
  
           

              
              
            
            
相關推薦
			   
            
            
            
 

    

    
    
線段樹維護區間(平方和，立方和)修改區間(加，賦值，乘)
     
 
							
							
							題目地址
/*
* @Author: hannibal
* @Date:   2018-08-07 10:42:26
* @Last Modified by:   hannibal
* @Last Modified time: 2018-08-07 17:08: 

  
 

    

    
    
hdu4578線段樹維護平方和,立方和(加,乘,賦值)或者珂朵莉樹
     
 quic   date   ++   long   pda   lower   bool   代碼   fin   一開始我用線段樹去做,結果debug了一個下午和晚上QAQ,後來學了珂朵莉樹,發現這題原來可以這麽簡單的寫,發出兩個代碼對比一下

#include<bits/stdc++.h> 

  
 

    

    
    
hdu4578線段樹維護平方和,立方和(加,乘,賦值)
     
 
                #include<bits/stdc++.h>
#define ll  unsigned long long
using namespace std;
const ll INFINITE = INT_MAX;
const ll MAXNUM = 100005*4; 

  
 

    

    
    
線段樹維護區間最大連續和
     
 
							
							
							在我的理解中，線段樹就是線段組成的樹，越靠近根部的層，線段的長度越長，一個結構體（樹的每一個點的左邊界，右邊界以及要維護的值）和三個函式（build——拆的過程），update（更新的過程，更新的時候要先二分找到最底層的區間，然後需要向上將所有的區間都更新了）， 

  
 

    

    
    
[Noi2016]區間[離散化+線段樹維護+決策單調性]
     
 fin   include   efi   cmp   http   說明   int   min   unique   
4653: [Noi2016]區間
Time Limit: 60 Sec  Memory Limit: 256 MBSubmit: 621  Solved: 329[Submit][ 

  
 

    

    
    
FJUT3568 中二病也要敲程式碼（線段樹維護區間連續最值）題解
     
 題意：有一個環，有1~N編號，m次操作，將a位置的值改為b，問你這個環當前最小連續和多少（不能全取也不能不取） 
思路：用線段樹維護一個區間最值連續和。我們設出兩個變數Lmin，Rmin，Mmin表示區間左邊最小連續和，右邊最小連續和，區間最小連續和，顯然這可以通過這個方式更新維護。 
 
現在我們已經可以維 

  
 

    

    
    
jsk Star War （線段樹維護區間最小最大值 + 二分）
     
  
 
 
 Description 
 公元20XX年，人類與外星人之間的大戰終於爆發。 
 現有一個人類軍團，由n名士兵組成，第i個士兵的戰鬥力值對應一個非負整數ai (1 \leq i \leq n1≤i≤n)。 
 有一天，某個戰力爆表的外星人NaN單獨向地球人宣戰，已知它的戰力值為k (1 \leq 

  
 

    

    
    
ACM-ICPC 2018 徐州賽區網路預賽 H. Ryuji doesn't want to study （線段樹維護字首和的字首和）
     
  
  
  
 題意  很簡答就是每次查詢 (L,R)的時候 ，查詢 
 
  
   
    a[L]∗(R−L+1)+a[L+1]∗(R−L−1)+a[L+2]∗(R−L−2)....a[R]∗1
   
    
     
      a
     
     
      [
     
    

  
 

    

    
    
最敏捷的機器人（線段樹維護區間最值）
     
 題面： 
Wind設計了很多機器人。但是它們都認為自己是最強的，於是，一場比賽開始了……機器人們都想知道誰是最敏捷的，於是它們進行了如下一個比賽。首先，他們面前會有一排共n個數，它們比賽看誰能最先把每連續k個數中最大和最小值寫下來，當然，這些機器人運算速度都很，它們比賽的是誰寫得快。但是Wind也想知道答案， 

  
 

    

    
    
ACM-ICPC 2018 徐州賽區網路預賽 H. Ryuji doesn't want to study （線段樹維護字首和的字首和）
     
 
							
							
							題意 
很簡答就是每次查詢 (L,R)的時候 ，查詢 a[L]∗(R−L+1)+a[L+1]∗(R−L−1)+a[L+2]∗(R−L−2)....a[R]∗1a[L]∗(R−L+1)+a[L+1]∗(R−L−1)+a[L+2]∗(R−L−2)....a[R]∗1 

  
 

    

    
    
HDU 3974 和 POJ 3321（線段樹維護DFS序）
     
  
 
 先介紹下DFS序 
 樹的dfs序就是用來維護一系列樹上的問題的，這類問題主要是解決一棵樹上的所有後代結點資訊的更改和祖先結點有關，主要先通過dfs來記錄一個樹的每一個頂點的出入時間戳，來控制它子樹上的所有結點的狀態的更新。 
 用in陣列記錄每個結點入棧的時間，out陣列記錄出棧時間（如果一個結點 

  
 

    

    
    
Fast Arrangement （線段樹 維護區間最大值 lazy標記）
     
 
								
								            
							
							
							Fast Arrangement
題意：
有一列火車同一時刻只能承載K個人，然後有Q個人要買票（給出Q組區間表示要買票的區間），（注意： 火車行駛區間為a-b的範圍），問那個人可以成功買到票（先到先得） 

  
 

    

    
    
（線段樹維護字首和）Ryuji doesn't want to study
     
 
                Ryuji is not a good student, and he doesn't want to study. But there are n books he should learn, each book has its knowledge a[i]a[i].

U 

  
 

    

    
    
HDU 5726 GCD （線段樹維護區間gcd）
     
 
								
								            
						
                
GCD
Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)
Total Submi 

  
 

    

    
    
hdu1540（線段樹維護連續區間模型）
     
 
                
/*
translation:
	打地道戰，n個村莊用地道連成一條直線，鬼子有時破壞掉一個村莊，這時八路軍要修好這個據點。現在要求詢問任意一個村莊，
	得出這個村莊現在與幾個村莊相連，包括它本身。
solution:
	線段樹維護連續子區間
	將村莊抽象成一個點，正常下值為 

  
 

    

    
    
區間MEX 線段樹維護mex陣列
     
 
							
							
							NKOJ 4254 區間MEX

問題描述


  給你一個長度為n的數列，元素編號1到n，第i個元素值為Ai。現在有m個形如(L,R)的提問，你需要回答出區間[L,R]的mex值。即求出區間[L,R]中沒有出現過的最小的非負整數。


輸入格式


  第一行 

  
 

    

    
    
P4198 樓房重建 （線段樹維護區間上升子序列）
     
 
                題目描述



小A的樓房外有一大片施工工地，工地上有N棟待建的樓房。每天，這片工地上的房子拆了又建、建了又拆。他經常無聊地看著窗外發呆，數自己能夠看到多少棟房子。

為了簡化問題，我們考慮這些事件發生在一個二維平面上。小A在平面上(0,0)點的位置，第i棟樓房可以用一條連線 

  
 

    

    
    
[永不止步-2017]_區間第K大-線段樹維護
     
 
                
	把線段樹的結點的資料域設定為vector型別即可別的操作為區間更新模板
	思路就是這樣runtime error暫時沒改對
	
	#include<bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f 

  
 

    

    
    
Codeforces Round #426 (Div. 2) D. The Bakery(線段樹維護dp)
     
 src   lap   codeforce   blank   com   date   close   scanf   logs   題目鏈接： Codeforces Round #426 (Div. 2) D. The Bakery
題意：
給你n個數，劃分為k段，每段的價值為這一段不同的數的個數，問如何 

  
 

    

    
    
HDU 6155 Subsequence Count 線段樹維護矩陣
     
 input   ans   ons   tom   thml   other   family   卡常   子序列   Subsequence Count
Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 256000/256000 K (J