题解不按顺序给出
目录：
ycz的妹子
题目背景：
(ycz)有很多很多的妹子（(ycz)：瞎说）
题目描述：
机房神犇(ycz)有n个青梅竹马，她们分别住在1~n号城市中。小时候的她们美丽可爱，但是由于女大十八变，有些妹子的颜值发生了变化，但是十分重感情的(ycz)神犇不忍心抛弃她们，于是记录下来了她们颜值变化的值，我们用(C\, x\, y)表示第(x)个城市的妹子的颜值下降了(y)（(y)有可能是负数）。长大之后的(ycz)非常有魅力，有许多妹子被(ycz)迷得神魂颠倒，我们用(I\, x\, y)表示第(x)个城市有一个妹子喜欢上了(ycz)，她的颜值为(y)（(y)有可能是负数，但是(ycz)来者不拒）。但在中途有一些妹子和(ycz)吵架了，于是就分手了，我们用(D\, x)表示第(x)个妹子和(ycz)分手了。
最近神犇(ycz)要去全国各地找他的妹子们，为了方便计算，我们珂以把(ycz)的妹子所在的城市当作是一条直线，并且挨在一起。神犇(ycz)由于忙于和他的妹子们联系此时已经很累了，于是交给你一个这样的任务：他想知道他在某个时间去找他的所有妹子们珂以获得多大的愉悦度，这个愉悦度为他找的妹子的颜值数，你要做的就是求出这个愉悦度之和（注意长大后妹子们的颜值可能为负数/滑稽）。
注意：每个城市只允许有一个妹子，也就是说后来喜欢上(ycz)的妹子会赶走之前这个城市喜欢(ycz)的妹子（一城不容二女）。
输入格式：
第一行两个正整数(n)和(m) ((1<=n<=100000))
第二行为(n)个整数(a_i)，表示小时候(ycz)的青梅竹马的颜值((1<=a_i<=10^9))
接下来(m)行，每行为一条信息，每条信息可能是下面的一种:
(C\, x\, y)表示第(x)个城市的妹子的颜值下降了(y)
(I\, x\, y)表示在第(x)个城市有一个颜值为(y)的妹子迷上了(ycz)
(D\, x)表示第(x)个妹子和(ycz)分手了
(Q)表示(ycz)现在想知道如果现在去找他所有的妹子们珂以获得多大的愉悦度
说明：妹子们居住的城市编号最大为(5 imes 10^5)
输出格式：
对于每一个(Q)输出一个整数
样例输入：
5 10
1 2 3 4 5
Q
C 3 2
Q
I 6 6
Q
D 4
Q
C 4 2
I 7 9
Q
样例输出：
15
13
19
15
22
提示与说明：
对于30%的数据(1<=n,m<=10)
对于70%的数据(1<=n,m<=1000)
对于100%的数据(1<=n,m<=100000,1<=a_i,y<=10^9)
题解：
30分：我也不知道怎么拿啊...(暴力分xiajb给的)
70分：直接暴力就好了...
code:点此跳过此代码
/*Program from Luvwgyx*/
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=100010;
int n,m,a[maxn];char s[10];
int read(){
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
int main(){
n=read();m=read();
for(int i=1;i<=n;i++)a[i]=read();
for(int i=1;i<=m;i++){
scanf("%s",s+1);
if(s[1]=='C'){
int x=read(),y=read();
a[x]-=y;
}
if(s[1]=='I'){
int x=read(),y=read();
a[x]=y;
}
if(s[1]=='D'){
int cnt=0,pos,x=read();
for(int i=1;i<=maxn-10;i++)
if(a[i]){cnt++;if(cnt==x){pos=i;break;}}
a[pos]=0;
}
if(s[1]=='Q'){
int ans=0;
for(int i=1;i<=maxn-10;i++)
ans+=a[i];
printf("%d
",ans);
}
}
return 0;
}
100分：
线段树中存颜值和以及妹子的数量就好了，具体的代码去看。
PS：应Wolfycz墙裂要求，此处贴上他的代码（好丑啊）我代码在下面，泥萌觉得哪个好看就看哪个吧QAQ
code:点此跳过此代码
/*program from Wolfycz*/
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define inf 0x7f7f7f7f
using namespace std;
typedef long long ll;
typedef unsigned int ui;
typedef unsigned long long ull;
inline int read(){
int x=0,f=1;char ch=getchar();
for (;ch<'0'||ch>'9';ch=getchar()) if (ch=='-') f=-1;
for (;ch>='0'&&ch<='9';ch=getchar()) x=(x<<1)+(x<<3)+ch-'0';
return x*f;
}
inline void print(int x){
if (x>=10) print(x/10);
putchar(x%10+'0');
}
const int N=1e5;
int val[(N<<1)+10];
struct Segment{
#define ls (p<<1)
#define rs (p<<1|1)
struct node{
int cnt;ll sum;
void insert(int _sum,int _cnt){sum=_sum,cnt=_cnt;}
node(){sum=cnt=0;}
}tree[(N<<3)+10];
friend node operator +(const node &x,const node &y){
node z;
z.sum=x.sum+y.sum;
z.cnt=x.cnt+y.cnt;
return z;
}
void build(int p,int l,int r){
if (l==r){
tree[p].insert(val[l],(bool)val[l]);
return;
}
int mid=(l+r)>>1;
build(ls,l,mid),build(rs,mid+1,r);
tree[p]=tree[ls]+tree[rs];
}
void change(int p,int l,int r,int x,int v){
if (l==r){
tree[p].insert(v,1);
return;
}
int mid=(l+r)>>1;
if (x<=mid) change(ls,l,mid,x,v);
else change(rs,mid+1,r,x,v);
tree[p]=tree[ls]+tree[rs];
}
void insert(int p,int l,int r,int x,int v){
if (l==r){
tree[p].sum-=v;
return;
}
int mid=(l+r)>>1;
if (x<=mid) insert(ls,l,mid,x,v);
else insert(rs,mid+1,r,x,v);
tree[p]=tree[ls]+tree[rs];
}
void Delete(int p,int l,int r,int x){
if (l==r){
tree[p].insert(0,0);
return;
}
int mid=(l+r)>>1;
if (x<=tree[ls].cnt) Delete(ls,l,mid,x);
else Delete(rs,mid+1,r,x-tree[ls].cnt);
tree[p]=tree[ls]+tree[rs];
}
}Tree;
char s[2];
int main(){
int n=read(),m=read();
for (int i=1;i<=n;i++) val[i]=read();
Tree.build(1,1,N<<1);
for (int i=1;i<=m;i++){
scanf("%s",s);
if (s[0]=='C'){
int x=read(),y=read();
Tree.insert(1,1,N<<1,x,y);
}
if (s[0]=='I'){
int x=read(),y=read();
Tree.change(1,1,N<<1,x,y);
}
if (s[0]=='D'){
int x=read();
Tree.Delete(1,1,N<<1,x);
}
if (s[0]=='Q') printf("%lld
",Tree.tree[1].sum);
}
return 0;
}
/*Program from Luvwgyx*/
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const ll maxn=2e5+10;
ll n,m,a[maxn];char s[10];
struct node{ll cnt,sum;}tree[maxn<<2];
ll read(){
ll x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
void update(ll k){
tree[k].cnt=tree[k<<1].cnt+tree[k<<1|1].cnt;
tree[k].sum=tree[k<<1].sum+tree[k<<1|1].sum;
return;
}
void build(ll k,ll l,ll r){
if(l==r){
tree[k].sum=a[l];
if(a[l]!=0)tree[k].cnt++;
return ;
}
ll mid=(l+r)>>1;
build(k<<1,l,mid);
build(k<<1|1,mid+1,r);
update(k);
}
void change(ll k,ll l,ll r,ll x,ll y){
if(l==r){tree[k].sum-=y;/*prllf("
%d
",l);*/return ;}
ll mid=(l+r)>>1;
if(x<=mid)change(k<<1,l,mid,x,y);
else change(k<<1|1,mid+1,r,x,y);
update(k);
}
void add(ll k,ll l,ll r,ll x,ll y){
if(l==r){tree[k].sum=y;tree[k].cnt=1;/*prllf("
%d
",l);*/return ;}
ll mid=(l+r)>>1;
if(x<=mid)add(k<<1,l,mid,x,y);
else add(k<<1|1,mid+1,r,x,y);
update(k);
}
void del(ll k,ll l,ll r,ll x){
if(l>r)return ;//prllf("%d %d
",l,r);
if(l==r&&x==tree[k].cnt){tree[k].sum=0;tree[k].cnt--;return ;}
ll mid=(l+r)>>1;//prllf("%d
",tree[k<<1].cnt);
if(x<=tree[k<<1].cnt)del(k<<1,l,mid,x);
else del(k<<1|1,mid+1,r,x-tree[k<<1].cnt);
update(k);
}
/*
void dfs(ll k,ll l,ll r){
if(l==r){if(tree[k].sum!=0)prllf("%d:%d ",l,tree[k].sum);return ;}
ll mid=(l+r)>>1;
dfs(k<<1,l,mid);dfs(k<<1|1,mid+1,r);
}
*/
int main(){
//freopen("girl10.in","r",stdin);
//freopen("girl10.out","w",stdout);
n=read();m=read();
for(ll i=1;i<=n;i++)a[i]=read();
build(1,1,maxn-10);
while(m--){
scanf("%s",s);
if(s[0]=='C'){ll x=read(),y=read();change(1,1,maxn-10,x,y);}
if(s[0]=='I'){ll x=read(),y=read();add(1,1,maxn-10,x,y);}
if(s[0]=='D'){ll x=read();del(1,1,maxn-10,x);}
if(s[0]=='Q')printf("%lld
",tree[1].sum);
//dfs(1,1,maxn-10);puts("");
}
return 0;
}
/*
2 11
18 4
C 2 -1
C 1 0
I 19 -17
I 8 -1
D 3
I 19 12
Q
D 1
C 2 -2
D 2
Q
*/
lty loves 96！
题目背景：
众所周知，(lty)非常喜欢(96)这两个数字（想歪的现在马上面壁去），更甚于复读（人本复）！
题目描述：
由于爱屋及乌，因此，(lty)对于那些含有(96)的数也非常喜欢，而这里的含有不是一般的含有，而是具有以下性质的含有（三条都需要满足）：
- 这个数为一个(N)位数，且没有前置零
- 这个数中至少要出现(M)次(9)和(6)（例：(986996)中出现了(5)次，(9)出现了(3)次，(6)出现了(2)次，共计(5)次）
- 这个数存在任意连续的三位(A),(B),(C)，满足下面任意一条
- (A+B+C)为(9)或(6)
- ((A^2+B^2))%(C)为(9)或(6)，如果(C)为(0)，则该条件视为不满足
输入格式：
一行，两个数(N),(M)。
输出格式：
一个数，表示这样的数的个数。
样例输入：
3 1
样例输出：
452
提示与说明：
对于10%的数据，(N<=6)，
对于40%的数据，(N<=18)，
对于100%的数据，(N<=50)，(0<=M<=N)。
题解：
很简单的一道数位dp。
记录(f[i][a][b][j][k](a,bleq 9;i,jleq N;k=0 or 1))，表示一个(i)位数，第(i)位和第(i-1)位分别为(a,b)，这个数中有(j)个6和9，这个数是否出现过题目中描述的两个条件，有则(k=1)，没有则(k=0)。
转移：
当(a,b,c)不满足题目中条件时：
否则
(f)数组的类型？
如果你开(int)类型，你可以获得(10)分的好成绩。
如果你开(long long)类型，你可以获得(40)分的好成绩。
如果你开普通的未压位高精度类型，你可以获得(0)分的好成绩。
如果你使用任意压位高精度类型，你可以获得(100)分的好成绩。
emmmm，好像有(dalao)用(int128)水了(85pts)来着....
code:点此跳过此代码
/*Program from wym*/
#include<bits/stdc++.h>
const int maxn=50;
const int digit=9;
const int base=1000000000;
struct bignum{
int v[maxn/digit+2];
int clear(){
memset(v,0,sizeof v);
return 0;
}
int print() {
printf("%d",v[v[0]]);
for(int i=v[0]-1; i>0; --i){
int k=1;
for(int j=1; j<digit; ++j){
k*=10;
if(v[i]<k)putchar('0');
}printf("%d",v[i]);
}
puts("");
return 0;
}
};
bignum operator +(bignum a,bignum b){
bignum ans; ans.clear();
ans.v[0]=std::max(a.v[0],b.v[0]);
int r=0;
for(int i=1; i<=ans.v[0]; ++i){
ans.v[i]=a.v[i]+b.v[i]+r;
r=ans.v[i]/base;
ans.v[i]%=base;
}
if(r!=0)ans.v[++ans.v[0]]=r;
return ans;
}
bignum operator +=(bignum &a,bignum b){
a=a+b;
return a;
}
bignum f[maxn+1][10][10][maxn+1][2],ans;
int n,m;
int main(){
scanf("%d%d",&n,&m);
for(int i=0; i<=n; ++i)
for(int a=0; a<=9; ++a)
for(int b=0; b<=9; ++b)
for(int j=0; j<=n; ++j){
f[i][a][b][j][0].clear();
f[i][a][b][j][1].clear();
}
f[0][0][0][0][0].v[0]=1;
f[0][0][0][0][0].v[1]=1;
for(int i=1; i<=n; ++i)
for(int a=0; a<=9; ++a)
for(int b=0; b<=9; ++b)
for(int j=0; j<=i; ++j)
for(int c=0; c<=9; ++c){
int last_j=j-(((a==6)||(a==9))?1:0);
if((last_j<0)||(last_j>=i))continue;
if((i>=3)&&(((a+b+c==9)||(a+b+c==6))||((c!=0)&&(((a*a+b*b)%c==9)||((a*a+b*b)%c==6)))))
f[i][a][b][j][1]+=f[i-1][b][c][last_j][0]+f[i-1][b][c][last_j][1];
else {
f[i][a][b][j][0]+=f[i-1][b][c][last_j][0];
f[i][a][b][j][1]+=f[i-1][b][c][last_j][1];
}
}
ans.clear();
for(int i=1; i<=9; ++i)
for(int j=0; j<=9; ++j)
for(int k=m; k<=n; ++k)
ans+=f[n][i][j][k][1];
ans.print();
return 0;
}
mzf的考验
题目背景：
(mzf)立志要成为一个豪杰，当然，他也是一个(OIer)。
他希望自己除了会(OI)之外还会各种东西，比如心理学、吉他、把妹等等。
为了让自己有更大的魅力，他不驼背，不熬夜，整天锻炼，双目炯炯有神，是我们机房最不像(OIer)的人。
然而，在与我们格格不入若干天并且将《易经》研究透彻之后，承受不住我们对他另类的言论，他爆发了。
机房在那一刹那仿佛天塌地陷，世界末日。
题目描述：
八卦有乾、坤、震、巽、坎、离、艮、兑；
两两组合，一上一下，形成了六十四卦，每卦六爻，一共三百八十四爻。
爻分阴阳，阳爻性属阳刚，阴爻性属阴柔。天下之大，无奇不有。千奇百怪，皆出此处。
(mzf)研究透彻了易经之后，画出了(n)个奇怪的图案。他说那是他改进出来的更强大的卜卦体系。
每一个图案有二十行，每一行要么是阴爻((0))，要么是阳爻((1))，作为一个(OIer)，我们可以将卦象看成一个个二进制串；
他将(n)个图案画在了符纸上，然后进行(m)次操作：
操作1：翻转区间([l,r])的图案，比如((3,1,2,5))变成((5,2,1,3))；
操作2：(mzf)画地为卦，将([l,r])之间的卦象都异或上新画的那个卦象；
操作3：(mzf)会询问机房里的其他人([l,r])之间卦象代表的二进制数权值和。
如果不能正确回答每个操作(3)，那么机房风水格局将会改变，我们都将...！
由于(mzf)疯狂之下将我们都捆♂绑♂了起来，所以只能求求你来帮我们解决这个问题。
输入格式：
第一行两个正整数：(n)，(m)（(n)为序列长度，(m)为操作个数）
第二行(n)个正整数：(a[i]) (用(10)进制数表示每个卦象)((1<=i<=n))
接下来(m)行：每行首先一个正整数(opt)表示操作类型
- (opt==1)：两个正整数：(l)，(r)。请翻转区间([l,r])；
- (opt==2)：三个正整数：(l)，(r)，(d)。请将区间([l,r])中的所有卦象都异或卦象(d)。((0<=d<=10^5))
- (opt==3)：两个正整数：(l)，(r)。请查询区间([l,r])的卦象权值和。
输出格式：
对于每个 (opt==3) 的情况，输出一行答案。
样例输入：
8 9
4 6 2 1 7 9 10 2
1 1 4
3 1 6
2 4 5 2
3 1 6
2 1 5 8
3 1 6
2 5 7 10
3 4 7
3 1 10
样例输出：
29
29
69
24
59
提示与说明：
对于20%的数据，(n<=1000)，(m<=1000)
对于另外20%的数据，不存在操作(1)
对于另外20%的数据，保证(n)为(2)的次幂，且在操作(1)中，保证(l=i imes(2^j)+1,r=(i+1) imes(2^j))，其中(i)，(j)为任意值
对于100%的数据，(n<=10^5)，(m<=5 imes 10^4)，(1<=l<=r<=n)，(0<=d<2^{20})
题解：
方法一：
暴力模拟，期望得分：20。
code:点此跳过此代码
/*Program from Luvwgyx*/
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn=1010;
int n,m,a[maxn],b[maxn];
int read(){
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
int main(){
n=read();m=read();
for(int i=1;i<=n;i++)a[i]=read();
for(int i=1;i<=m;i++){
int opt=read(),l=read(),r=read(),d;
if(opt==1){
for(int i=l;i<=r;i++)b[i]=a[i];
for(int i=r;i>=l;i--)a[i]=b[l+r-i];
memset(b,0,sizeof(0));
}else if(opt==2){
d=read();
for(int i=l;i<=r;i++)a[i]^=d;
}else if(opt==3){
ll ret=0;
for(int i=l;i<=r;i++)ret+=a[i];
printf("%lld
",ret);
}
}
return 0;
}
方法二：
不存在操作1，那我们可以考虑用线段树来做。
普通的线段树：期望得分：20。
结合暴力，期望得分：40。
代码就不放了，自己去写一写吧
方法三：
有翻转操作了，但是这个反转的区间有限制，我们会发现这个区间其实就是线段树上的区间的左右端点，我们只需要找到这个区间，然后翻转就好了。
线段树+拆位异或...期望得分：40分
结合暴力，期望得分：60分
线段树内维护权值和以及二进制下数位的值，然后拆位异或，注意一下标记下传之类的就好了。
都是板子，直接上代码看吧...
code:点此跳过此代码
/*Program from lichenxi*/
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define mid ((s[x].l+s[x].r)>>1)
#define size (s[x].r-s[x].l+1)
int n,m,opt,x,y,c[21],ans[21],f[21];
long long d[300001],z,sum,w[21],a[300001];
struct oo{int l,r,rev,t[21],ls,rs;long long _XOR;}s[1200001];
int read(){
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
void gettwo(int *a,long long v)
{
int now=0;
while(v)
{
if(v&1)a[now]=1;
v>>=1;now++;
}
}
void update(int x){for(int i=0;i<20;i++)s[x].t[i]=s[s[x].ls].t[i]+s[s[x].rs].t[i];}
void build(int x,int l,int r)
{
s[x].l=l,s[x].r=r;
if(l==r){gettwo(s[x].t,d[l]);return ;}
s[x].ls=x<<1,s[x].rs=x<<1|1;
build(x<<1,l,mid);build(x<<1|1,mid+1,r);
update(x);
}
void pushrev(int x)
{
int ca=s[s[x].ls].r-s[s[x].ls].l,ka=s[s[x].rs].r-s[s[x].rs].l;
s[s[x].ls].l=s[x].l;s[s[x].ls].r=s[s[x].ls].l+ca;
s[s[x].rs].l=s[s[x].ls].r+1;s[s[x].rs].r=s[s[x].rs].l+ka;
}
void pushdown_rev(int x)
{
s[s[x].ls].rev^=1,s[s[x].rs].rev^=1,s[x].rev^=1;
swap(s[x].ls,s[x].rs);
pushrev(x);
}
void push_XOR(int *c,int x){for(int i=0;i<20;i++)if(c[i])s[x].t[i]=size-s[x].t[i];}
void pushdown_XOR(int x)
{
s[s[x].ls]._XOR^=s[x]._XOR;s[s[x].rs]._XOR^=s[x]._XOR;
for(int i=0;i<20;i++)f[i]=0;gettwo(f,s[x]._XOR);
push_XOR(f,s[x].ls);push_XOR(f,s[x].rs);
s[x]._XOR=0;
}
void putrev(int x,int l,int r)
{
pushrev(x);
if(s[x].rev)pushdown_rev(x);
if(s[x]._XOR)pushdown_XOR(x);
if(s[x].l==l&&s[x].r==r)
{
swap(s[x].ls,s[x].rs);
s[s[x].ls].rev^=1;s[s[x].rs].rev^=1;
int ca=s[s[x].ls].r-s[s[x].ls].l,ka=s[s[x].rs].r-s[s[x].rs].l;
s[s[x].ls].l=s[x].l;s[s[x].ls].r=s[s[x].ls].l+ca;
s[s[x].rs].l=s[s[x].ls].r+1;s[s[x].rs].r=s[s[x].rs].l+ka;
return ;
}
if(l<=mid)putrev(s[x].ls,l,r);
if(r>mid)putrev(s[x].rs,l,r);
update(x);
}
void change(int x,int l,int r,long long d)
{
pushrev(x);
if(s[x].rev)pushdown_rev(x);
if(s[x]._XOR)pushdown_XOR(x);
if(l<=s[x].l&&r>=s[x].r){push_XOR(c,x),s[x]._XOR^=d;return ;}
if(l<=mid)change(s[x].ls,l,r,d);
if(r>mid)change(s[x].rs,l,r,d);
update(x);
}
void get(int x,int l,int r)
{
pushrev(x);
if(s[x].rev)pushdown_rev(x);
if(s[x]._XOR)pushdown_XOR(x);
if(l<=s[x].l&&r>=s[x].r)
{
for(int i=0;i<20;i++)ans[i]+=s[x].t[i];
return ;
}
if(l<=mid)get(s[x].ls,l,r);
if(r>mid)get(s[x].rs,l,r);
update(x);
}
int main()
{
// freopen("data.in","r",stdin);
// freopen("data.out","w",stdout);
w[0]=1;
for(int i=1;i<22;i++)w[i]=w[i-1]<<1;
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)scanf("%lld",&d[i]);
if(n<=1000&&m<=1000){
while(m--){
int opt=read(),l=read(),r=read(),x;
if(opt==1){
for(int i=l;i<=r;i++)a[i]=d[i];
for(int i=r;i>=l;i--)d[i]=a[l+r-i];
memset(a,0,sizeof(a));
}else if(opt==2){
x=read();
for(int i=l;i<=r;i++)d[i]^=x;
}else if(opt==3){
long long ret=0;
for(int i=l;i<=r;i++)ret+=d[i];
printf("%lld
",ret);
}
}
return 0;
}
build(1,1,n);
while(m--)
{
scanf("%d%d%d",&opt,&x,&y);
if(opt==1)putrev(1,x,y);
if(opt==2)
{
scanf("%lld",&z);
for(int i=0;i<20;i++)c[i]=0;
gettwo(c,z);change(1,x,y,z);
}
if(opt==3)
{
for(int i=0;i<20;i++)ans[i]=0;get(1,x,y);
sum=0;for(int i=0;i<20;i++)sum+=ans[i]*w[i];
printf("%lld
",sum);
}
}
}
方法四：
当操作(1)没有限制后，我们可以想到用平衡树，开(20)棵平衡树就可以(TLE)了（常数小的说不定卡的过）。每次翻转的时候我们需要翻转(20)棵平衡树，这样非常不优。那么我们就把每个数的二进制存到平衡树的节点上，每次异或的时候再拿出来用，翻转操作也只要做一次就好。理论时间复杂度(O(20 m logn))，常数小的基本上可以过
至于有(dalao)写的是(splay)却(55)或(60)啥的我也不知道啥情况啊....
我真的没有想卡常啊...可能是我的(splay)常数比较小吧...说不定用(fhq treap)常数小一些，但是我不会写
期望得分：100分。
code:点此跳过此代码
/*program from Wolfycz*/
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define inf 0x7f7f7f7f
using namespace std;
typedef long long ll;
typedef unsigned int ui;
typedef unsigned long long ull;
inline int read(){
int x=0,f=1;char ch=getchar();
for (;ch<'0'||ch>'9';ch=getchar()) if (ch=='-') f=-1;
for (;ch>='0'&&ch<='9';ch=getchar()) x=(x<<1)+(x<<3)+ch-'0';
return x*f;
}
inline void print(int x){
if (x>=10) print(x/10);
putchar(x%10+'0');
}
const int N=5e5;
int v[N+10],g[20];
void Extract(int x,int *cnt){for (int i=0;i<20;i++) cnt[i]=x&1,x>>=1;}
struct Splay{
#define ls(x) tree[x][0]
#define rs(x) tree[x][1]
#define T(x) (rs(f[x])==x)
ll sum[N+10];
int tree[N+10][2],f[N+10],size[N+10],cnt[N+10][20],Xor[N+10],val[N+10],A[20];
bool rev[N+10];
int root,len;
void updata(int x){
size[x]=size[ls(x)]+size[rs(x)]+1;
sum[x]=sum[ls(x)]+sum[rs(x)]+val[x];
for (int i=0;i<20;i++) cnt[x][i]=cnt[ls(x)][i]+cnt[rs(x)][i]+((val[x]>>i)&1);
}
void move(int x){
int fa=f[x],son=tree[x][T(x)^1];
tree[x][T(x)^1]=fa;
tree[fa][T(x)]=son;
if (son) f[son]=fa;
f[x]=f[fa];
if (f[x]) tree[f[x]][T(fa)]=x;
f[fa]=x;
updata(fa),updata(x);
}
void splay(int x){
while (f[x]){
if (f[f[x]]) T(f[x])==T(x)?move(f[x]):move(x);
move(x);
}
root=x;
}
void Add_rev(int x){
if (!x) return;
swap(ls(x),rs(x));
rev[x]^=1;
}
void Add_xor(int x,int v){
if (!x) return;
Extract(v,A);
ll res=0;
for (int i=0;i<20;i++){
if (A[i]) cnt[x][i]=size[x]-cnt[x][i];
res+=1ll*cnt[x][i]*g[i];
}
sum[x]=res;
val[x]^=v,Xor[x]^=v;
}
void pushdown(int x){
if (rev[x]){
Add_rev(ls(x));
Add_rev(rs(x));
rev[x]=0;
}
if (Xor[x]){
Add_xor(ls(x),Xor[x]);
Add_xor(rs(x),Xor[x]);
Xor[x]=0;
}
}
int find(int x,int i){
pushdown(i);
if (size[ls(i)]+1==x) return i;
if (x<=size[ls(i)]) return find(x,ls(i));
return find(x-size[ls(i)]-1,rs(i));
}
void reverse(int l,int r){
int x=find(l,root),y=find(r+2,root);
splay(x),splay(y);
if (f[x]!=root) move(x);
Add_rev(rs(x));
// updata(x),updata(y);
}
void _xor(int l,int r,int v){
int x=find(l,root),y=find(r+2,root);
splay(x),splay(y);
if (f[x]!=root) move(x);
Add_xor(rs(x),v);
updata(x),updata(y);
}
ll Query(int l,int r){
int x=find(l,root),y=find(r+2,root);
splay(x),splay(y);
if (f[x]!=root) move(x);
return sum[rs(x)];
}
void build(int fa,int l,int r,int &x){
if (l>r) return;
int mid=(l+r)>>1;
x=++len,f[x]=fa,val[x]=v[mid];
if (l==r){
size[x]=1,sum[x]=val[x];
Extract(val[x],cnt[x]);
return;
}
build(x,l,mid-1,ls(x));
build(x,mid+1,r,rs(x));
updata(x);
}
void init(int n){
len=2,root=1;
rs(1)=size[1]=2,f[2]=size[2]=1;
build(2,1,n,ls(2));
updata(2),updata(1);
}
}Tree;
int main(){
g[0]=1;
for (int i=1;i<20;i++) g[i]=g[i-1]<<1;
int n=read(),m=read();
for (int i=1;i<=n;i++) v[i]=read();
Tree.init(n);
for (int i=1;i<=m;i++){
int type=read(),l=read(),r=read();
if (type==1) Tree.reverse(l,r);
if (type==2) Tree._xor(l,r,read());
if (type==3) printf("%lld
",Tree.Query(l,r));
}
return 0;
}
hby与tkw的基情
题目背景：
基情恒久远，一对永流传。
(hby)和(tkw)是一对好基友，他们经常在一起做♂游♂戏 (/滑稽)。
题目描述：
他们喜欢玩字符串游戏，尤其喜欢玩回文串。每次(hby)会给出一个数(n)，那么(tkw)就需要给出(Ans=sumlimits_{i=1}^ni imes s[i] imes[i\%2])
其中(s[i])代表长度为(i)的数字回文串的个数，最后面是(bool)表达式
不过由于(tkw)最近学(ycz)找妹纸去了，于是他就将这个问题交给了你，如果你不能在(1s)内答出来，那么(hby)和(tkw)的基情将会破裂！（不过那样tkw就可以安心地找妹纸了）
由于答案会非常大，所以你只需要输出答案(\%10^9+7)的值即可
输入格式:
第一行一个整数(T)。
接下来(T)行，每行一个数(n)。
输出格式：
共(T)行，每行代表一个答案
样例输入：
2
1
3
样例输出：
26
2054
提示与说明：
对于10%的数据：(n<=5)
对于另外20%的数据：(sum n<=10^7)
对于另外20%的数据：(T=1)
对于100%的数据：(T<=5 imes 10^5,n<=10^9)
题解：
10%：随便玩玩就好
30%：发现:(s[i]=s[i-2] * 26) ，(s[1]=26)，递推即可
50%：由于只有一组数据，可以用矩阵快速幂优化以上的递推。
100%：开了(O2)说不定矩阵快速幂就过了。。。
然后我们要求的式子其实是
其中(m=frac{n+1}{2})
这东西显然一个差比数列，裂项相减一下
直接快速幂即可
/*program from Wolfycz*/
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define inf 0x7f7f7f7f
using namespace std;
typedef long long ll;
typedef unsigned int ui;
typedef unsigned long long ull;
inline int read(){
int x=0,f=1;char ch=getchar();
for (;ch<'0'||ch>'9';ch=getchar()) if (ch=='-') f=-1;
for (;ch>='0'&&ch<='9';ch=getchar()) x=(x<<1)+(x<<3)+ch-'0';
return x*f;
}
inline void print(int x){
if (x>=10) print(x/10);
putchar(x%10+'0');
}
const int p=1e9+7,inv=2.8e8+2,inv2=5.6e8+5;
int mlt(int a,int b){
int res=1;
for (;b;b>>=1,a=1ll*a*a%p) if (b&1) res=1ll*res*a%p;
return res;
}
int main(){
for (int T=read();T;T--){
int n=(read()+1)>>1;
int Ans=(1ll*(26+1ll*(2*n-1)*mlt(26,n+1)%p-2ll*(mlt(26,n+1)-26)%p*inv%p)*inv%p+p)%p;
printf("%d
",Ans);
}
return 0;
}
抓住czx
题目背景：
蒟蒻(lty)出了一道题，但是由于太弱了，所以希望喜欢鸽子的(czx)来帮他写一个(std)。由于(czx)又放鸽子去了，所以没有写(std)。蒟蒻(lty)觉得受到了学长的鄙视，所以决定去(czx)放鸽子的地方找他。
题意简述：
(czx)放鸽子的地方是一个公园，公园珂以看作是由(n)个点(m)条边组成的无向图(保证无自环），(lty)将从公园的入口（(b)号节点）进去寻找(czx)，而(czx)会在(a_i)个单位时间时变化位置到第(x)个节点去，在此之前(lty)已经知道了(czx)的具体位置和接下来他位置的变化方案，蒟蒻(lty)现在想知道他至少需要花多少时间找到(czx)。
输入格式：
第一行四个整数(n),(m),(b),(e),(b)和(e)的意义如题面所示。
接下来(m)行，每行三个整数(x,y,z)，表示(x)到(y)之间有一条长为(z)的边。
第(m+1)行一个整数(T)，表示(czx)位置变化的次数。
接下来(T)行，每行两个整数(a_i)和(x)，表示(czx)将在第(a_i)个单位时间时移动到第(x)个点上去。
输出格式：
一个整数表示最短所需时间。
样例输入：
6 9 1 6
1 2 1
1 3 3
1 4 4
2 3 2
3 6 6
4 5 6
2 5 9
3 5 7
5 6 2
3
10 3
8 5
9 2
样例输出：
9
样例解释：
在开始的时候就直接走到(2)号节点，然后等到(czx)过来。总花费时间(9)个单位时间。
数据范围：
对于30%的数据，(n<=100),(m<=1000),(T<=100)
对于另外30%的数据，(T=0)
对于100%的数据，(n<=10^5),(m<=5 imes10^5),(T<=10^5)
数据结果保证在(int)范围内
题解：
关于(SPFA)：它还活着（本来想卡的。。。）
30%：瞎(jb)弄吧，其实我也不知道这个部分分有啥意义。。。
另外30%：单源最短路板子。（其实是数据太大放不上。。。强行减小数据大小）
100%：首先求一遍单源最短路，如果我们出发去找(czx)，但是半路中(czx)就(TP)（(MC)玩家）走了，那么我们不如最开始就找之后的某个点。这样就可以找到一个点，在(czx)下次(TP)前抓住他即可。
/*program from Wolfycz*/
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define inf 0x7f7f7f7f
using namespace std;
typedef long long ll;
typedef unsigned int ui;
typedef unsigned long long ull;
inline int read(){
int x=0,f=1;char ch=getchar();
for (;ch<'0'||ch>'9';ch=getchar()) if (ch=='-') f=-1;
for (;ch>='0'&&ch<='9';ch=getchar()) x=(x<<1)+(x<<3)+ch-'0';
return x*f;
}
inline void print(int x){
if (x>=10) print(x/10);
putchar(x%10+'0');
}
const int N=1e5,M=1e6;
int pre[(M<<1)+10],now[N+10],child[(M<<1)+10],val[(M<<1)+10];
int h[N+10],deep[N+10];
bool vis[N+10];
int n,m,S,tot;
struct S1{
int x,T;
void insert(int _x,int _T){x=_x,T=_T;}
bool operator <(const S1 &a)const{return T<a.T;}
}A[N+10];
void join(int x,int y,int z){pre[++tot]=now[x],now[x]=tot,child[tot]=y,val[tot]=z;}
void SPFA(int x){
int head=0,tail=1;
memset(deep,63,sizeof(deep));
h[1]=x,vis[x]=1,deep[x]=0;
while (head!=tail){
if (++head>N) head=1;
int Now=h[head];
for (int p=now[Now],son=child[p];p;p=pre[p],son=child[p]){
if (deep[son]>deep[Now]+val[p]){
deep[son]=deep[Now]+val[p];
if (!vis[son]){
if (++tail>N) tail=1;
vis[h[tail]=son]=1;
}
}
}
vis[Now]=0;
}
}
int main(){
n=read(),m=read(),S=read(),A[0].insert(read(),0);
for (int i=1;i<=m;i++){
int x=read(),y=read(),z=read();
join(x,y,z),join(y,x,z);
}
SPFA(S);
int T=read();
for (int i=1;i<=T;i++){
int x=read(),y=read();
A[i].insert(y,x);
}
A[++T].insert(0,inf);
sort(A,A+1+T);
for (int i=0;i<=T;i++){
if (deep[A[i].x]<A[i+1].T){
printf("%d
",max(deep[A[i].x],A[i].T));
break;
}
}
return 0;
}