Luvwgyx的娱乐场-题解

题解不按顺序给出

目录:

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\times 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\n",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\n",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("\n%d\n",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("\n%d\n",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\n",l,r);
    if(l==r&&x==tree[k].cnt){tree[k].sum=0;tree[k].cnt--;return ;}
    ll mid=(l+r)>>1;//prllf("%d\n",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\n",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,b\leq 9;i,j\leq N;k=0\ or\ 1)\),表示一个\(i\)位数,第\(i\)位和第\(i-1\)位分别为\(a,b\),这个数中有\(j\)个6和9,这个数是否出现过题目中描述的两个条件,有则\(k=1\),没有则\(k=0\)

转移:

\(a,b,c\)不满足题目中条件时:

\[ f[i][a][b][j][0]+=f[i-1][b][c][j_{last}][0],f[i][a][b][j][1]+=f[i-1][b][c][j_{last}][1] \]

否则
\[ f[i][a][b][j][1]+=f[i-1][b][c][j_{last}][0]+f[i-1][b][c][j_{last}][1] \]
\(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\)表示操作类型

  1. \(opt==1\):两个正整数:\(l\)\(r\)。请翻转区间\([l,r]\)
  2. \(opt==2\):三个正整数:\(l\)\(r\)\(d\)。请将区间\([l,r]\)中的所有卦象都异或卦象\(d\)\((0<=d<=10^5)\)
  3. \(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\times(2^j)+1,r=(i+1)\times(2^j)\),其中\(i\)\(j\)为任意值

对于100%的数据,\(n<=10^5\)\(m<=5\times 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\n",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\n",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\n",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\n",Tree.Query(l,r));
    }
    return 0;
}

hby与tkw的基情

题目背景:

基情恒久远,一对永流传。

\(hby\)\(tkw\)是一对好基友,他们经常在一起做♂游♂戏 (/滑稽)。

题目描述:
他们喜欢玩字符串游戏,尤其喜欢玩回文串。每次\(hby\)会给出一个数\(n\),那么\(tkw\)就需要给出\(Ans=\sum\limits_{i=1}^ni\times s[i]\times[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\times 10^5,n<=10^9\)

题解:

10%:随便玩玩就好

30%:发现:\(s[i]=s[i-2] * 26\)\(s[1]=26\),递推即可

50%:由于只有一组数据,可以用矩阵快速幂优化以上的递推。

100%:开了\(O2\)说不定矩阵快速幂就过了。。。
然后我们要求的式子其实是
\[\sum\limits_{i=1}^m (2\times i-1)\times 26^i\]
其中\(m=\frac{n+1}{2}\)
这东西显然一个差比数列,裂项相减一下
\[Ans=\sum\limits_{i=1}^m(2\times i-1)\times 26^i\]
\[26Ans=\sum\limits_{i=1}^m(2\times i-1)\times 26^{i+1}\]
\[25Ans=26+(2n-1)\times 26^{n+1}-2\times\sum\limits_{i=1}^m26^i\]
\[25Ans=26+(2n-1)\times 26^{n+1}-2\times\dfrac{26^{n+1}-26}{25}\]
直接快速幂即可

/*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\n",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\times10^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\n",max(deep[A[i].x],A[i].T));
            break;
        }
    }
    return 0;
}

转载于:https://www.cnblogs.com/Luvwgyx/p/9538102.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值