Mas 的童年

题目描述:不知道传送门有没有用?

反正就是对于每个前缀序列求一个断点,使得断点左右两个区间的 分别的异或和 的和最大

分析

jzoj 原题?

但是我 TM 代码没存账号也过期了啊!

然后只能绞尽脑汁思考做法,然后在左边 bzt 大仙各种提示下终于打出了一个 n log n 的做法

考虑暴力肯定非常可做于是我们来个暴力完事

我们想想当前的前缀序列的异或和为 x 的情况下,我们考虑从高位到低位去使贡献最大化:

那么假设当前处理的位置是 i :

如果 x 的第 i 位为 0 ,并且能从前面找到某个断点使得左边右边两个区间的当前位都是 1 ,那这样肯定是最优的,我们给找寻的断点加上这个限制

另外,我们考虑是从高位处理到低位的,我们在考虑 y 位的贡献时前面位的贡献必须已经最大化了,那么我们要在前面高位的限制下寻找新的断点,找到就增加答案并增加限制,找不到就拉倒什么也别干

如果 x 的第 i 位为 1 ,那么怎么分都没关系,贡献都是 1 ,限制就比较小

然后我们江当前的前缀作为一个断点,江它的所有真子集都变为 1 表示存在这样的断点

但是这样的复杂度是 \(n m\) 的, m 就是 \(a_i\) 的数据范围(说白了就是 \(n^2\))

那么怎么优化呢?

我们考虑从断点更新方面无法优化了,那么考虑这个子集能不能预处理

预处理的方法就是对于每种状态我们记录一下它最早可以达到的时间点,然后寻找的时候直接看当前限制状态记录的时间点是否在当前时间的前面,如果是的话就说明存在满足条件的断点

这样做有什么优化到的地方么? 我们发现可以用高维前缀和传递信息,复杂度$ n log m $

于是高维前缀和就弄掉

//by Judge
#include<cstdio>
#include<cstring>
#include<iostream>
#define Rg register
#define fp(i,a,b) for(Rg int i=(a),I=(b)+1;i<I;++i)
#define fd(i,a,b) for(Rg int i=(a),I=(b)-1;i>I;--i)
#define lowbit(x) (x&-x)
#define ll long long
using namespace std;
const int M=2e6+3;
#ifndef Judge
#define getchar() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
#endif
char buf[1<<21],*p1=buf,*p2=buf;
inline bool cmin(int& a,int b){return a>b?a=b,1:0;}
inline bool cmax(int& a,int b){return a<b?a=b,1:0;}
inline int read(){ int x=0,f=1; char c=getchar();
for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;
for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;
} char sr[1<<21],z[20];int CCF=-1,Z;
inline void Ot(){fwrite(sr,1,CCF+1,stdout),CCF=-1;}
inline void print(int x,char chr=' '){
if(CCF>1<<20)Ot();if(x<0)sr[++CCF]=45,x=-x;
while(z[++Z]=x%10+48,x/=10);
while(sr[++CCF]=z[Z],--Z);sr[++CCF]=chr;
} int n,mx,a[M],pre[M],tmp;
int main(){ n=read(),memset(pre,0x3f,sizeof pre);
fp(i,1,n) a[i]=read()^a[i-1]; fd(i,n,1) pre[a[i]]=i,cmax(mx,a[i]);
fd(i,mx,1) fp(j,0,20) if((i>>j)&1) cmin(pre[i^(1<<j)],pre[i]);
fp(i,1,n){ tmp=0;
fd(j,20,0) if(((a[i]>>j)&1)^1)
if(pre[tmp|(1<<j)]<=i) tmp|=1<<j;
print((a[i]^tmp)+tmp);
} return sr[CCF]='\n',Ot(),0;
}

Z的家乡

题目描述:不知道传送门有没有用?

反正就是给你一棵树每次在树上加一个叶子,然后要你在每次加入的节点后 dfs

序的最大字典序...(讲不大清楚?)

然后就是一眼看出这玩意儿超级像 LCT 中的 access ,每次都把当前节点所在的链优先级提到最大了,但是不会做...

没搞懂,抄了篇比较短的代码感觉还不错...就当是锻炼手速了!

//by Judge
#include<bits/stdc++.h>
#define fp(i,a,b) for(register int i=(a),I=(b)+1;i<I;++i)
using namespace std;
const int mod=998244353;
const int M=1e5+3;
typedef int arr[M];
#ifndef Judge
#define getchar() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
#endif
char buf[1<<21],*p1=buf,*p2=buf;
inline int mul(int a,int b){return 1ll*a*b%mod;}
inline int inc(int a,int b){return a+b>=mod?a+b-mod:a+b;}
inline int read(){ int x=0,f=1; char c=getchar();
for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;
for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;
} char sr[1<<21],z[20];int CCF=-1,Z;
inline void Ot(){fwrite(sr,1,CCF+1,stdout),CCF=-1;}
inline void print(int x,char chr='\n'){
if(CCF>1<<20)Ot();if(x<0)sr[++CCF]=45,x=-x;
while(z[++Z]=x%10+48,x/=10);
while(sr[++CCF]=z[Z],--Z);sr[++CCF]=chr;
} int n,pw[M];
struct data{ int siz,ans;
data operator +(data y)const{
return (data){siz+y.siz,inc(ans,mul(y.ans,pw[siz]))};
}
};
struct seg{ int l,r; data x; };
namespace seg_T{ int x,tot,root[M],stk[M*20]; data y; seg t[M*20];
#define mid ((l+r)>>1)
#define lch t[rt].l
#define rch t[rt].r
inline int New(){return *stk?stk[(*stk)--]:++tot;}
inline void R(int& rt){t[rt].l=t[rt].r=0,t[rt].x=(data){0,0},stk[++*stk]=rt,rt=0;}
inline void add(int l,int r,int& rt){ if(!rt) rt=New(); if(l==r) return t[rt].x=y,void();
if(x<=mid) add(l,mid,lch); else add(mid+1,r,rch); t[rt].x=t[rch].x+t[lch].x;
}
inline void del(int l,int r,int& rt){ if(l==r) return R(rt);
if(x<=mid) del(l,mid,lch); else del(mid+1,r,rch);
t[rt].x=t[rch].x+t[lch].x; if(!t[rt].x.siz) R(rt);
}
inline void add(int rt,int _x,data _y){x=_x,y=_y,add(1,n,root[rt]);}
inline void del(int rt,int _x){x=_x,del(1,n,root[rt]);}
inline data get(int x){return t[root[x]].x;}
} using namespace seg_T;
struct LCT{ arr f,l,r,id; data t[M];
inline bool isroot(int x){return l[f[x]]^x&&r[f[x]]^x;}
inline void pushup(int x){ id[x]=r[x]?id[r[x]]:x,t[x]=t[r[x]]+get(x)+(data){1,x}+t[l[x]];}
inline void rotate(int x){ int y=f[x],z=f[y];
if(l[z]==y) l[z]=x; else if(r[z]==y) r[z]=x; f[x]=z,f[y]=x;
if(l[y]==x) f[l[y]=r[x]]=y,r[x]=y; else f[r[y]=l[x]]=y,l[x]=y;
pushup(y),pushup(x);
}
inline void splay(int x){
for(int y,z;!isroot(x);rotate(x))
if(y=f[x],z=f[y],!isroot(y))
rotate((l[z]==y)^(l[y]==x)?x:y);
}
inline void init(){fp(i,1,n) pushup(i);}
inline void link(int u,int v){ int p1,p2,x;
for(f[x=v]=u,p1=0;x;x=f[x]){ splay(x);
if(r[x]) add(x,id[r[x]],t[r[x]]);
if(r[x]=p1,p1&&p1^v) del(x,p2);
p1=x,p2=id[x],pushup(x);
} splay(v),print(t[v].ans);
}
}T;
int main(){ n=read(),pw[0]=1;
fp(i,1,n) pw[i]=mul(pw[i-1],257933);
T.init(); fp(i,2,n) T.link(read(),i);
return Ot(),0;
}

战棋游戏

题目描述:不知道传送门有没有用?

反正就是不可做

分析

NTT,不会,我只会打板子

代码?算了吧...(比赛比完都没几个 A 的)

NOI2019 SX 模拟赛 no.5的更多相关文章

  1. NOIP模拟赛20161022

    NOIP模拟赛2016-10-22 题目名 东风谷早苗 西行寺幽幽子 琪露诺 上白泽慧音 源文件 robot.cpp/c/pas spring.cpp/c/pas iceroad.cpp/c/pas ...

  2. 大家AK杯 灰天飞雁NOIP模拟赛题解/数据/标程

    数据 http://files.cnblogs.com/htfy/data.zip 简要题解 桌球碰撞 纯模拟,注意一开始就在袋口和v=0的情况.v和坐标可以是小数.为保险起见最好用extended/ ...

  3. NOIP模拟赛-2018.11.5

    NOIP模拟赛 好像最近每天都会有模拟赛了.今天从高二逃考试跑到高一机房,然而高一也要考试,这回好像没有拒绝的理由了. 今天的模拟赛好像很有技术含量的感觉. T1:xgy断句. 好诡异的题目,首先给出 ...

  4. NOIP模拟赛-2018.10.22

    模拟赛 今天第一节课是历史,当然是不可能上的,一来到机房发现今天高二考试... 老师说以后可能还要给高一考...那还不如现在跟着做好了,毕竟在学长学姐中垫底显得没那么丢人 这套题风格挺奇怪的...为什 ...

  5. noip模拟赛 蒜头君救人

    分析:之前的一道模拟赛题是dp+dfs,这道题是dp+bfs. 我们设f[stu][i][j]为当前状态为stu,走到(i,j)的答案,考虑怎么设计stu,每个人的状态有3种:要么在原地,要么被背着, ...

  6. NOIP模拟赛 17.10.10

    初次见面(firstmeet)[题目背景]雾之湖边,静得可怕.露米娅出神凝望.黑白连衣裙,像极了绽放的墨黑和洁白的莲.身边的雾之湖,倒映着血色天空.酒红的双眸,映照一切.低声浅笑,双臂伸直,她悄无声息 ...

  7. NOI模拟赛 Day1

    [考完试不想说话系列] 他们都会做呢QAQ 我毛线也不会呢QAQ 悲伤ING 考试问题: 1.感觉不是很清醒,有点困╯﹏╰ 2.为啥总不按照计划来!!! 3.脑洞在哪里 4.把模拟赛当作真正的比赛,紧 ...

  8. NOIP第7场模拟赛题解

    NOIP模拟赛第7场题解: 题解见:http://www.cqoi.net:2012/JudgeOnline/problemset.php?page=13 题号为2221-2224. 1.car 边界 ...

  9. contesthunter暑假NOIP模拟赛第一场题解

    contesthunter暑假NOIP模拟赛#1题解: 第一题:杯具大派送 水题.枚举A,B的公约数即可. #include <algorithm> #include <cmath& ...

随机推荐

  1. springboot中通用mapper结合mybatis generator的使用

    通用mapper就是指的是  tk.mybatis  包下的.这个是通用mapper就是说自动生成的dao层需要继承这个框架提供的mapper类.而我们之前用的org.mybatis这个最开始是普通的 ...

  2. HDU 1030(三角数阵 数学)

    题意是问在给定的三角形数阵中从一个数到另一个数所要跨过的边数. 最初的时候很迷,除了发现每层的数字个数与层数间的关系和每层数最后一个数与层数的关系外什么也没看出来,打算先求出数字所在的层数,然后计算到 ...

  3. 认证加密算法php hash_hmac和java hmacSha1的问题

    public class Test{ public static void main(String[] args) throws Exception { String postString = &qu ...

  4. GCC编译器原理(二)------编译原理一:ELF文件(1)

    二.ELF 文件介绍 2.1 可执行文件格式综述 相对于其它文件类型,可执行文件可能是一个操作系统中最重要的文件类型,因为它们是完成操作的真正执行者.可执行文件的大小.运行速度.资源占用情况以及可扩展 ...

  5. mybatis-configxml样板

    <?xml version="1.0" encoding="UTF-8" ?><!DOCTYPE configuration PUBLIC & ...

  6. [C++]2-6 排列

    /* 排列(Permutation) 用1,2,3,...,9组成3个三位数abc,def和ghi,每个数字恰好使用一次,要求abc:def:ghi = 1:2:3. 按照"abc def ...

  7. [C++]线性链表之顺序表<一>

    顺序表中数据元素的存储地址是其序号的线性函数,只要确定了存储顺序表的起始地址(即 基地址),计算任意一个元素的存储地址的时间是相等的,具有这一特点的存储结构称为[随机存储]. 使用的基本数据结构:数组 ...

  8. jdk下载地址

    1.5   http://www.oracle.com/technetwork/java/javasebusiness/downloads/java-archive-downloads-javase5 ...

  9. 第25月第17天 django rest framwork authentication /tmp/mysql.sock

    1.authentication https://www.django-rest-framework.org/api-guide/authentication/#authentication 2.dj ...

  10. 网络XSS攻击和CSRF攻击原理及防范

    网络XSS攻击和CSRF攻击原理及防范 原文地址:http://www.freebuf.com/articles/web/39234.html 随着Web2.0.社交网络.微博等等一系列新型的互联网产 ...