[NOIP11.1模拟赛]补番报告
Preface
昨天开始补某科学的超电磁炮S 感觉今天就好了点,炮姐赛高
T1 一开始一直想欧拉定理&ex欧拉定理,结果估计70分,数组开小了GG,看了正解发现是我学傻了
T2 一看就是数据结构,之前某次模拟还做过区间位运算线段树但是不敢打,只敲了个前缀和预计50结果有个地方没膜GG
T3 想了蛮久,很有意思的一道题目,考场上画了图发现我们可以将放坐标的那个点看成树根,这样它覆盖了一条近似链的玩意,想搞波DP发现不会,于是打了个骗分的玩意,就是我猜和节点的儿子个数有关,我直接用儿子个数除2向上取整结果...OJ上15,source里30.发现ZZ地没判点的情况,然后发现更ZZ的是我特意判了链的情况,但是我输出的是0
10分就这么走了GG
T1 pow
学傻了的我一直想欧拉定理...结果solution告诉你只要会快速幂就可以A这题了
a既然是定值我们就分块打表啊
由于b最大是1e12,按1e6拆分.
预处理出\({a^1}... {a^{1e6}}\),丢进一个数组,在预处理出\(a^{i \times 1e6}\),丢进另外一个数组
这样所有的b都能拼出来了
还是比较妙的,这告诉我们有的时候去搞旧算法不如想新trick
代码
const int maxn=1000005;
const int inf = 0x7fffffff;
ll a,q,k;
ll b,l,m,c,p;
ll pre[maxn][2],sum[10000005];
inline ll ksm(ll aa,ll cc){
ll ans=1;
while(cc){
if(cc&1)ans=ans*a%p;
aa=aa*aa%p;
cc=cc>>1;
}
return ans;
}
int main(){
FI(pow)FO(pow)
read(a),read(p),read(q),read(k);
read(b),read(l),read(m),read(c);
pre[0][0]=1;
for(ri i=1;i<=1000000;i++){
pre[i][0]=pre[i-1][0]*a%p;
}
pre[1][1]=pre[1000000][0];
for(ri i=2;i<=1000000;i++){
pre[i][1]=pre[i-1][1]*pre[1][1]%p;
}
for(ri i=1;i<=q;i++){
b=(b*m+c)%l;
if(b<=1000000)sum[i]=sum[i-1]^pre[b][0];
else {
ll id=b/1000000;
sum[i]=sum[i-1]^((pre[b%1000000][0]*pre[id][1])%p);
}
}
int kk=k;
while(kk<=q){
printf("%lld\n",sum[kk]);
kk+=k;
}
return 0;
}
T2 seg
英文是seg 中文是tree 就是告诉你用线段树(segment tree)啦
对于区间按位与操作,类比区间取膜的时间复杂度分析(给学弟讲过结果自己还忘了)
每个数最多操作31次,n只有1e5,因此是资瓷的
我们只要维护一个区间按位或判断需不需要修改就好了,其余情况暴力递归修改
2操作就trival了,会线段树的都会
3操作也比较正常,你可以先思考部分分怎么拿,你把这个期望平方和(实际上这个式子就是\(\sum_{i=l}^r (a_i+a_j)^2 (j=l,l+1...r)\))大力展开
发现维护一个区间平方和和区间和就好了
然后有个神坑的地方
就是在判断是否要按位与时的位运算的优先级
if((or_sum[now]&(~dta))==0)return ;
if(or_sum[now]&(~dta)==0)return ;
请您判断上面哪个是对的
答案是第一个,但是你写第二个不会对答案正确性产生影响,但是每次它都会递归下去修改使得时间复杂度大大提高
可见位运算优先级之低
代码
const int maxn=100005;
const int inf=0x7fffffff;
const int P=998244353;
int a[maxn],n,q;
ll sum[maxn<<2];
int eq_sum[maxn<<2];
int or_sum[maxn<<2];
inline void up(int now){
sum[now]=sum[now<<1]+sum[now<<1|1];
eq_sum[now]=(eq_sum[now<<1]+eq_sum[now<<1|1]);
if(eq_sum[now]>P)eq_sum[now]-=P;
or_sum[now]=or_sum[now<<1]|or_sum[now<<1|1];
return ;
}
void build(int now,int l,int r){
if(l==r){
or_sum[now]=sum[now]=a[l];
eq_sum[now]=1ll*a[l]*a[l]%P;
return ;
}
int mid=(l+r)>>1;
build(now<<1,l,mid);
build(now<<1|1,mid+1,r);
up(now);return ;
}
int dta,t,L,R;
void update(int now,int l,int r){
if((or_sum[now]&(~dta))==0)return ;//运算符顺序!!!
if(l==r){
a[l]=a[l]&dta;
or_sum[now]=sum[now]=a[l];
eq_sum[now]=1ll*a[l]*a[l]%P;
return ;
}
int mid=(l+r)>>1;
if(L<=mid)update(now<<1,l,mid);
if(mid<R)update(now<<1|1,mid+1,r);
up(now);return ;
}
ll ans1=0,ans2=0;
void query(int now,int l,int r){
if(L<=l&&r<=R){
ans1+=sum[now];
ans2=(ans2+eq_sum[now]);
if(ans2>P)ans2-=P;
return ;
}
int mid=(l+r)>>1;
if(L<=mid)query(now<<1,l,mid);
if(mid<R)query(now<<1|1,mid+1,r);
return ;
}
int main(){
int x,y,opt;
FO(seg)
//freopen("seg6.in","r",stdin);
//freopen("wtf.out","w",stdout);
read(n);
for(ri i=1;i<=n;i++)read(a[i]);
build(1,1,n);
read(q);
while(q--){
read(opt),read(L),read(R);
ans1=ans2=0;
if(opt==1){
read(dta);
update(1,1,n);
}
if(opt==2){
query(1,1,n);
printf("%lld\n",ans1);
}
if(opt==3){
query(1,1,n);
ans1=ans1%P;
printf("%lld\n",(((ans2<<1)%P*(R-L+1)%P+(ans1<<1)*ans1%P))%P);
}
}
return 0;
}
T3 beacon
一道Topcoder上有趣的题目
我是这么想的,除了一点之外的情况,答案肯定是大与等于1的,我们不妨先钦定一个点,在上面放一个信标,然后以它为根遍历整棵树,显然此时深度相同的点都是非法的
考虑这种情况:有两个兄弟叶节点,那么它们此时是非法的,显然在除这两个叶节点之外的任何一点放置信标这两点还是非法的(它们的深度还是相同),所以我们必须在这两个叶子节点之一放一个信标
稍微拓展一下:假若有三个兄弟叶节点,那么类似的发现你必须在三个叶节点中放两个信标才可以;于是归纳假设发现对于n个兄弟叶节点你必须在之中放n-1个信标
那么对于不是叶节点的点呢?你会发现这时候它们似乎都已经是合法的了
除了链的情况,链的情况下由于最底端只有一个叶子节点不会统计答案,但是实际上你会发现链实际上整体就可以看做一个叶节点处理.
于是按照上面的步骤\(O(N^2)\)就好了
满分做法看不懂,这里给出题解,不知哪位大佬可以帮忙解释一下
如何做到 O(n)? 我们先特判链的情况答案为 1, 然后找到任意一个度数大于 2 的节点, 可以证
明这个点一定不需要放置信标. 于是以这个点作根 O(n) 的贪心即可. 证明如下:
深度相同的点对证明同上, 只考虑深度不同的点对. 如果它们在一颗子树中, 由于度数大于 2 所
以一定有另一颗子树的一个信标把他们区分开. 如果在不同的子树中, 有两种情况:
一个在没放信标的子树中, 一个在放了的子树中. 显然还存在另一个子树放了信标, 由于深度不
同他们会被这个信标区分开.
两个都在放了信标的子树中. 如果根的度数大于 3 则同上. 度数等于 3 时, 如果他们没有被区分
开, 一定是他们先汇集到了一个节点上, 然后走到同一个信标上. 这个点一定是一条奇链的中点, 且
不是根 (由于深度不同), 是在两个子树之一中唯一的. 那么他们走到另一个信标就一定有一个点走
了冤枉路, 既另一个信标可以区分出他们
70分代码
const int maxn=1000005;
const int inf=0x7fffffff;
int n;
struct Edge{
int ne,to;
}edge[maxn<<1];
int h[maxn],num_edge=1;
inline void add_edge(int f,int to){
edge[++num_edge].ne=h[f];
edge[num_edge].to=to;
h[f]=num_edge;
}
int ans=inf,sum=0;
bool is_lef[maxn],on_chain[maxn];
void dfs(int now,int fa){
int v,cnt=0,tot=0;
is_lef[now]=0,on_chain[now]=0;
for(ri i=h[now];i;i=edge[i].ne){
v=edge[i].to;
if(v==fa)continue;
is_lef[now]=1,tot++;
dfs(v,now);
if(!is_lef[v]||on_chain[v])cnt++;
}
if(tot==1&&cnt==1)on_chain[now]=1;
if(cnt>1)sum+=cnt-1;
return ;
}
int deg[maxn];
int main(){
int x,y;
//FO(beacon)
read(n);
if(n==1){puts("0");return 0;}
for(ri i=1;i<n;i++){
read(x),read(y);
add_edge(x,y),add_edge(y,x);
}
for(ri rt=1;rt<=n;rt++){
sum=0;
dfs(rt,0);
ans=min(ans,sum+1);
}
printf("%d\n",ans);
return 0;
}
[NOIP11.1模拟赛]补番报告的更多相关文章
- 【HHHOJ】NOIP模拟赛 捌 解题报告
点此进入比赛 得分: \(30+30+70=130\)(弱爆了) 排名: \(Rank\ 22\) \(Rating\):\(-31\) \(T1\):[HHHOJ260]「NOIP模拟赛 捌」Dig ...
- 【HHHOJ】NOIP模拟赛 玖 解题报告
点此进入比赛 得分: \(100+20+100=220\)(还不错) 排名: \(Rank\ 16\) \(Rating\):\(+20\) \(T1\):[HHHOJ263]「NOIP模拟赛 玖」三 ...
- 冲刺Noip2017模拟赛3 解题报告——五十岚芒果酱
题1 素数 [问题描述] 给定一个正整数N,询问1到N中有多少个素数. [输入格式]primenum.in 一个正整数N. [输出格式]primenum.out 一个数Ans,表示1到N中有多少个素 ...
- 冲刺Noip2017模拟赛2 解题报告——五十岚芒果酱
题1 牛跑步(running) [题目描述] 新牛到部队,CG 要求它们每天早上搞晨跑,从 A 农场跑到 B 农场.从 A 农场到 B 农场中有 n- 个路口,分别标上号,A 农场为 号,B 农场为 ...
- 冲刺Noip2017模拟赛1 解题报告——五十岚芒果酱
题1 国际象棋(chess) [问题描述] 有N个人要参加国际象棋比赛,该比赛要进行K场对弈.每个人最多参加2场对弈,最少参加0场对弈.每个人都有一个与其他人都不相同的等级(用一个正整数来表示).在对 ...
- CH Round #55 - Streaming #6 (NOIP模拟赛day2)解题报告
T1九九归一 描述 萌蛋在练习模n意义下的乘法时发现,总有一些数,在自乘若干次以后,会变成1.例如n=7,那么5×5 mod 7=4,4×5 mod 7=6,6×5 mod 7=2,2×5 mod 7 ...
- 20161005 NOIP 模拟赛 T2 解题报告
beautiful 2.1 题目描述 一个长度为 n 的序列,对于每个位置 i 的数 ai 都有一个优美值,其定义是:找到序列中最 长的一段 [l, r],满足 l ≤ i ≤ r,且 [l, r] ...
- 20161003 NOIP 模拟赛 T2 解题报告
Weed duyege的电脑上面已经长草了,经过辨认上面有金坷垃的痕迹. 为了查出真相,duyege 准备修好电脑之后再进行一次金坷垃的模拟实验. 电脑上面有若干层金坷垃,每次只能在上面撒上一层高度为 ...
- 20161005 NOIP 模拟赛 T3 解题报告
subset 3.1 题目描述 一开始你有一个空集,集合可以出现重复元素,然后有 Q 个操作 1. add s 在集合中加入数字 s. 2. del s 在集合中删除数字 s.保证 s 存在 3. c ...
随机推荐
- Maven Multi-Module Example
Maven Multi-Module - 国内版 Binghttps://cn.bing.com/search?q=Maven+Multi-Module&qs=n&form=QBRE& ...
- linux内核在挂载ramdisk的过程中报错"RAMDISK: incomplete write (10739 != 32768)"如何处理?
1. 原因 ramdisk大小不够 2. 解决方法 在启动变量bootargs中添加参数"ramdisk_size=10000000"即可
- 解决Wireshark安装Npcap组件失败
解决Wireshark安装Npcap组件失败 从Wireshark 3.0开始,Npcap取代Winpcap组件,成为Wireshark默认的网卡核心驱动.由于该组件属于驱动程序,所以安装时候容易 ...
- Python3基础 交换两个变量的值
Python : 3.7.3 OS : Ubuntu 18.04.2 LTS IDE : pycharm-community-2019.1.3 ...
- 算法习题---5.8Unix ls命令(Uva400)
一:题目 输入正整数n以及n个文件名,将这n个文件名按照ASCII优先升序排列,按列优先方式左对齐输出. 注意:文件名最长的为M,那么最右列字符串长度为M,其他列的长度为M+2 注意:一列最多允许出现 ...
- 查看所使用的Linux系统是32位还是64 位的方法
方法一:getconf LONG_BIT # getconf LONG_BIT 1 1 我的Linux是32位!!! 方法二:arch # arch 1 1 显示 i686 就是32位,显示 x86_ ...
- 青葱的岁月 Mybatis JdbcType与Oracle、MySql数据类型对应列表
Mybatis JdbcType Oracle MySql JdbcType ARRAY JdbcType BIGINT BIGINT JdbcType BINARY JdbcTy ...
- iOS popToViewController具体用法
[self.navigationController popToViewController:[self.navigationController.viewControllers objectAtIn ...
- C#操作Access时Parameters集合的使用方法(转)
按照C#操作Sql Server数据库的方式是不一样的,使用时发现占位符数据会混乱. 查阅相关资料及测试后,总结出这样一套使用方法. 方法: public bool TsqlExecute(strin ...
- Docker入门之docker-compose
参考:https://blog.51cto.com/9291927/2310444 一,Docker-compose简介 1,Docker-compose简介 Docker-Compose项目是Doc ...