「Wdsr-1」笨蛋结构
题目背景
众所周知,琪露诺是笨蛋。
题目描述
琪露诺希望维护一个长度为 $n$ 的整数序列 $a$,初始值都为 $0$。
现在琪露诺想要进行 $q$ 次操作,每次选择序列中的一段区间 $[s,s+l-1]$ 并给出两个数字 $w,k$,使对所有的 $i \in [1,l]$,$a_{s+i-1}$ 加上 $w\times i^k$ 。
琪露诺不希望 $k$ 很大,因此她给出了一个整数 $m$,满足 $0\le k\le m$。
为了不让头脑简单的琪露诺感到困惑,你只需要输出 依次进行完所有操作后,序列中的每个数字对 $2^{64}$ 取模(即 $\text{unsigned long long}$ 自然溢出)后的结果即可。
为了帮助你更好的理解题意,这里给出一段伪代码:
$$\def{\b}#1{\textbf{ #1 }}\def{\t}#1{\text{ #1 }}\def{\s}{\quad}
\def{\l}{\underline{\kern{300pt}}\cr[-10pt]}
\def{\r}{\overline{\underline{\kern{300pt}}}}
\begin{aligned}
&\r\cr&\b{Algorithm:}\t{An easy structure}\cr[-13pt]&\l\cr
&\begin{aligned}
\t{1.}&\b{input}n,m,q \cr
\t{2.}&\b{for}i=1\b{to} q \b{do} \cr
\t{3.}&\s\b{input} s,l,w,k \cr
\t{4.}&\s\b{for} j=1 \b{to} l \b{do}\cr
\t{5.}&\s\s a[s+j-1] \gets a[s+j-1]+w\times \t{pow}(j,k) \cr
\t{6.}&\s\b{end}\cr
\t{7.}&\b{end}\cr
\t{8.}&\b{for} i=1 \b{to} n \b{do}\cr
\t{9.}&\s\b{output} a[i]\cr
\t{10.}&\b{end}\cr
\end{aligned}\cr[-12pt]
&\r\end{aligned}
%Made by @离散小波变换° .
%You can find his contributions by searching "JoesSR".
$$
其中 $\rm pow(a,b)$ 的含义为 $a^b$。
输入输出格式
输入格式
请调用下方代码的 `input(n,m,q,S,L,W,K)` 来输入 $n,m,q,s_i,l_i,w_i,k_i$。下标**从一开始**。
其中 $s,l,w,k$ 的含义与题目描述保持一致。
输出格式
请调用下方代码的 `output(n,R)` 进行输出。其中 $R_i$ 为所有操作结束后的数列,下标**从一开始**。
输入输出样例
输入样例 #1
10 0 5 233
输出样例 #1
6942214367
输入样例 #2
1000 9 500 6666
输出样例 #2
7636746723064426256
说明
#### 样例一说明
生成的数据为:
```plain
10 0 5
7 1 1558211206 0
1 3 401324017 0
4 5 235225636 0
6 4 2137131141 0
1 2 3791175968 0
```
它的结果是:
```plain
4192499985 4192499985 401324017 235225636 235225636 2372356777 3930567983 2372356777 2137131141 0
```
---
#### 数据生成&数据输出
```cpp
typedef unsigned long long u64;
typedef unsigned int u32;
u32 MT[624],idx;
void _init(u32 seed){
MT[0]=seed; idx=0; for(int i=1;i<624;++i)
MT[i]=(0x6c078965*(MT[i-1]^((MT[i-1])>>30)+i));
}
void _gene(){
for(int i=0;i<624;++i){
int x=MT[i]&0x80000000+(MT[(i+1)%624]&0x7fffffff);
MT[i]=MT[(i+397)%624]^(x>>1);
if(x&2)MT[i]^=0x9908b0df;
}
}
u32 _calc(){
if(!idx) _gene(); int x=MT[idx];
x^=x>>11,x^=(x<<7)&(0x9d2c5680);
x^=(x<<15)&0xefc60000,x^=x>>18;
idx=(idx+1)%624; return x;
}
u64 _get(){u64 ret=_calc()*_calc(); return ret;}
u64 _get(u64 _l,u64 _r){return _get()%(_r-_l+1ull)+_l;}
void input(int &_n,int &_m,int &_q,int *_S,int *_L,u64 *_W,int *_K){
u32 seed; scanf("%d%d%d%u",&_n,&_m,&_q,&seed); _init(seed); int i=1;
if(_n>100) for(;i<=_q/4;++i){
int _a=_get(1,_n-100),_b=_get(_a+_m,_a+_m+1),_l=_b-_a+1,_k=_get(0,_m);
u64 _w=_get(); _S[i]=_a,_L[i]=_l,_W[i]=_w,_K[i]=_k;
}
if(_n>100) for(;i<=_q/2;++i){
int _a=_get(1,100),_b=_get(_n-100,_n),_l=_b-_a+1,_k=_get(0,_m);
u64 _w=_get(); _S[i]=_a,_L[i]=_l,_W[i]=_w,_K[i]=_k;
}
for(;i<=_q;++i){
int _a=_get(1,_n),_b=_get(1,_n); if(_a>_b) swap(_a,_b);
int _l=_b-_a+1,_k=_get(0,_m); u64 _w=_get();
_S[i]=_a,_L[i]=_l,_W[i]=_w,_K[i]=_k;
}
}
void output(int n,u64 *R){
u64 ret=n^_get(); for(int i=1;i<=n;i++) ret^=_get()+R[i];
printf("%llu\n",ret);
}
```
其中,调用 `input()` 读入数据;调用 `output()` 输出数据。
**请勿在任何时候调用除了`input`和`output`外的函数,且这两个函数只能调用一次**。
---
#### 数据范围
共 $20$ 个测试点,满足如下条件:
$$\def\arraystretch{1.5}\begin{array}{|c|c|c|c|}\hline
\textbf{编号} & n & m & q \\ \hline
[1,3] & \le 3\times 10^3 & =9 & \le 3\times 10^3 \\\hline
[4,5] & \le 3\times 10^5 & =0 & \le 3\times 10^5 \\\hline
[6,9] & \le 3\times 10^5 & =1 & \le 3\times 10^5 \\\hline
[10,13] & \le 3\times 10^5 & =2 & \le 3\times 10^5 \\\hline
[14,16] & \le 3\times 10^5 & =9 & \le 3\times 10^5 \\\hline
[17,20] & \le 5\times 10^5 & =9 & \le 1\times 10^6 \\\hline
\end{array}$$
其中,$[l,r]$ 表示编号为 $l,l+1,\cdots,r-1,r$ 的测试点。
对于 $100\%$ 的数据,满足 $1\le l_i \le l_i+s_i-1 \le n,0\le k_i\le m,0 \le w\le 2^{64}-1$。