【FCS NOI2018】福建省冬摸鱼笔记 day5
第五天,也是讲课的最后一天。
数据结构专题,讲师:杨志灿
他的blog我似乎找不到了……以前肯定是在百度博客里面。但是现在百度博客消失了。
PPT做的很有感觉,说了很多实用的技巧。
我觉得其实是收获最大的一天,因为听懂了XD
中午划水
下午的题也非常良心,然而@ghostfly233和@Melacau说他们做过原题???
问题就是非常卡常!n=10^5的nlogn题目,400组询问,4秒?6亿复杂度啊!
然而就是跑过了,我也没办法,树状数组复杂度不满……
矩阵快速幂也卡……非常难受。
【T1】
题意:有\(1\leq a\leq b\leq c\leq n\),问多少个\((a,b,c)\)三元组满足\(a+b^2\equiv c^3(mod\;k)\)。
题解:枚举b,发现a+b^2的范围是b^2+1到b^2+b是连续的一段,考虑计算这段中的c大于等于b的c^3的个数。
显然用树状数组维护,复杂度有点问题,可是就是AC了XD。代码非常的优美。
#include<cstdio>
#include<cstring>
int n,k,bit[100001];
long long ans;
inline void I(int i){for(++i;i<=k;i+=i&-i)++bit[i];}
inline int Q(int i){int s=0;for(++i;i;i-=i&-i)s+=bit[i];return s;}
int main(){
freopen("exclaim.in","r",stdin);
freopen("exclaim.out","w",stdout);
int T; scanf("%d",&T);
for(int t=1;t<=T;++t){
scanf("%d%d",&n,&k); ans=0;
memset(bit,0,sizeof bit);
int sum=0;
for(int i=n;i;--i){
I(1ll*i*i*i%k); ++sum;
long long l=1ll*i*i+1, r=1ll*i*i+i;
ans+=1ll*(r/k-l/k)*sum+Q(r%k)-Q(l%k-1);
}
printf("Case %d: %lld\n",t,ans);
}
return 0;
}
【T2】
题意:你玩一个多结局RPGgalgame(???),游戏流程类似树形结构,你总共通关了k次,每个场景你都会获得该场景的价值(可能还有好感度和浪费手纸?),但是多次经历该场景就只有一次价值了。
价值都是正整数。
问他怎样攻略能获得最大价值?
题解:有一个贪心思路:每次都走最大价值的路线,这样一定最终价值也最大。
如何证明?
我们先证明最大的一定要选吧,这样后面的归纳一下就完了。
先分类讨论:
假设我们没有选择最大的,那么我们从最大的这条路线往回回溯,遇到第一个有选取的场景,那么走这个场景的这条(任意一条)路线,都可以换成走最大的路线。
为什么?假设走到这个场景的路线只有一条,那么显然更优,因为前面的价值都一样,换成这条会更优。
假设有多条,那么其中任意一条换掉之后,多了最优的这一条的价值收益,并且自己这条的收益还不用完全减掉,因为还有其他的路线走过,所以比上一种情况还来得优。
那么我们就证明了,最大的路线一定要选,接下来把最大的路线的收益删去,变成0,那么这时再选取最大的还是一样的。于是证明完毕。
那么如何处理呢?包括查询最大收益和修改。
其实这就是一个线段树能够维护的:
把节点按照DFS序排序,一个点的所有孩子就是一段区间。
把所有的叶子节点(结局)按照DFS序扔进线段树,初始值就是路线收益,查询时直接查最大值,修改时,我们知道一个点的孩子对应一段区间,所以区间减法即可。
查询k次,每次O(1),修改最多n次,每次O(log 叶子节点数)。
那么这题就算做完了。
#include<cstdio>
#include<algorithm>
using namespace std;
int n,k,a[200001],fa[200001],t[200001],cnt;
int h[200001],to[200001],nxt[200001],tot;
inline void ins(int x,int y){nxt[++tot]=h[x];to[tot]=y;h[x]=tot;}
int L[200001],R[200001],vis[200001];
long long dat[524289],lazy[524289],ans;
int pos[524289];
void dfs(int u){
if(!h[u]) {L[u]=R[u]=++cnt, t[cnt]=u; return;}
L[u]=n; for(int i=h[u];i;i=nxt[i]) dfs(to[i]), L[u]=min(L[u],L[to[i]]), R[u]=max(R[u],R[to[i]]);
}
inline void pushdown(int i){
lazy[i<<1]+=lazy[i];
lazy[i<<1|1]+=lazy[i];
dat[i<<1]+=lazy[i];
dat[i<<1|1]+=lazy[i];
lazy[i]=0;
}
void build(int i,int l,int r){
pos[i]=l;
if(l==r) return;
int mid=l+r>>1;
build(i<<1,l,mid), build(i<<1|1,mid+1,r);
}
void Ins(int i,int l,int r,int a,int b,int v){
if(a<=l&&r<=b) {dat[i]+=v; lazy[i]+=v; return;}
if(r<a||b<l) return;
pushdown(i); int mid=l+r>>1;
Ins(i<<1,l,mid,a,b,v); Ins(i<<1|1,mid+1,r,a,b,v);
dat[i]=max(dat[i<<1],dat[i<<1|1]);
pos[i]=dat[i<<1]>=dat[i<<1|1]?pos[i<<1]:pos[i<<1|1];
}
int main(){
freopen("game.in","r",stdin);
freopen("game.out","w",stdout);
scanf("%d%d",&n,&k); k=min(k,n);
for(int i=1;i<=n;++i) scanf("%d",a+i);
for(int i=1,x,y;i<n;++i) scanf("%d%d",&x,&y), ins(x,y), fa[y]=x;
dfs(1);
build(1,1,cnt);
for(int i=1;i<=n;++i) Ins(1,1,cnt,L[i],R[i],a[i]);
for(int i=1;i<=k;++i){
int p=t[pos[1]];
ans+=dat[1];
int sum=0;
while(p&&!vis[p]) Ins(1,1,cnt,L[p],R[p],-a[p]), vis[p]=1, sum+=a[p], p=fa[p];
}
printf("%lld",ans);
return 0;
}
【T3】
题意:求有向图中长度小于k的环的个数,没有自环。
题解:简单矩阵,有著名定理支持:邻接矩阵的k次方就是从一个节点出发,经过k步到达另一个节点的路径数。
即求\(A^1+A^2+\cdots+A^{k-1}\)。
简单矩阵快速幂即可,加一个稀疏矩阵优化,跑的飞快。根本不需要分治。
#include<cstdio>
#include<cstring>
int n,n2,k,p,ans;
struct Mat{int fk[200][200];}M1,M2;
inline Mat operator*(Mat p1,Mat p2){
Mat p3;memset(p3.fk,0,sizeof p3.fk);
for(register int i=0;i<n2;++i)
for(register int k=0;k<n2;++k) if(p1.fk[i][k])
for(register int j=0;j<n2;++j)
p3.fk[i][j]=(p3.fk[i][j]+1ll*p1.fk[i][k]*p2.fk[k][j])%p;
return p3;
}
int main(){
freopen("tour.in","r",stdin);
freopen("tour.out","w",stdout);
scanf("%d",&n); n2=n<<1;
char ch;
for(register int i=0;i<n;++i) M1.fk[i][i]=M2.fk[i][i]=M2.fk[i+n][i+n]=1;
for(register int i=n;i<n2;++i) for(register int j=0;j<n;++j) (ch=getchar())=='Y'?M1.fk[i][j]=M1.fk[i][j+n]=1:(ch!='N'?--j:0);
scanf("%d%d",&k,&p); --k;
if(p==1){puts("0");return 0;}
while(k){
if(k&1) M2=M2*M1;
k>>=1; M1=M1*M1;
}
for(register int i=0;i<n;++i) ans=(ans+M2.fk[i+n][i])%p;
printf("%d",ans);
return 0;
}
【FCS NOI2018】福建省冬摸鱼笔记 day5的更多相关文章
- 【FCS NOI2018】福建省冬摸鱼笔记 day1
省冬的第一天. 带了本子,笔,一本<算法导论>就去了.惊讶于为什么同学不带本子记笔记. 他们说:“都学过了.”,果然这才是巨神吧. 第一天:数论,讲师:zzx 前几页的课件挺水,瞎记了点笔 ...
- 【FCS NOI2018】福建省冬摸鱼笔记 day6【FJOI 2018】福建省选混分滚蛋记 day1
记录一下day6发生的事情吧. 7:30 到达附中求索碑,被人膜,掉RP. 7:50 进考场,6楼的最后一排的最左边的位置,世界上最角落的地方,没有任何想法. 发现电脑时间和别人不一样,赶快调了一下. ...
- 【FCS NOI2018】福建省冬摸鱼笔记 day3
第三天. 计算几何,讲师:叶芃(péng). dalao们日常不记笔记.@ghostfly233说他都知道了,就盼着自适应辛普森积分. 我计算几何基础不好……然而还是没怎么讲实现,感觉没听什么东西进去 ...
- 【FCS NOI2018】福建省冬摸鱼笔记 day2
第二天. 同学还是不带本子记笔记.dalao. 第二天:图论,讲师:@ExfJoe 全程划水,前面都讲水算法[虽然我可能已经忘记了]什么最短路,Tarjan,最小生成树,2SAT,差分约束啥的,我现在 ...
- 【FCS NOI2018】福建省冬摸鱼笔记 day4
第四天. 动态规划专题,讲师:闫神 讲了一些DP优化技巧,然而思想难度好大啊……根本没想到能优化那地步,连DP方程都没有呢. 不过有几题我还是想明白了. 讲了单调队列,决策单调性,四边形不等式,斜率优 ...
- Hash 日常摸鱼笔记
本篇文章是Hash在信息学竞赛中的应用的学习笔记,分多次更新(已经有很多坑了) 一维递推 首先是Rabin-Karp,对于一个长度为\(m\)的串\(S\) \(f(S)=\sum_{i=1}^{m} ...
- 可持久化Treap 赛前摸鱼笔记
1.基本结构 随机化工具 unsigned int SEED = 19260817; //+1s inline int Rand(){ SEED=SEED*1103515245+12345; retu ...
- [摸鱼]cdq分治 && 学习笔记
待我玩会游戏整理下思绪(分明是想摸鱼 cdq分治是一种用于降维和处理对不同子区间有贡献的离线分治算法 对于常见的操作查询题目而言,时间总是有序的,而cdq分治则是耗费\(O(logq)\)的代价使动态 ...
- 【转】让Chrome化身成为摸鱼神器,利用Chorme运行布卡漫画以及其他安卓APK应用教程
下周就是十一了,无论是学生党还是工作党,大家的大概都会有点心不在焉,为了让大家更好的心不在焉,更好的在十一前最后一周愉快的摸鱼,今天就写一个如何让Chrome(google浏览器)运行安卓APK应用的 ...
随机推荐
- BZOJ5323 JXOI2018游戏(线性筛+组合数学)
可以发现这个过程非常类似埃氏筛,将在该区间内没有约数的数定义为质数,那么也就是求每种方案中选完所有质数的最早时间之和. 于是先求出上述定义中的质数个数,线性筛即可.然后对每个最短时间求方案数,非常显然 ...
- python selenium2 有关cookie操作实例及如何绕开验证码
1.先看一下cookie是啥 cookie是访问web时服务器记录在用户本地的一系列用户信息(比如用户登录信息),以便对用户进行识别 from selenium import webdriver im ...
- 【刷题】BZOJ 3732 Network
Description 给你N个点的无向图 (1 <= N <= 15,000),记为:1-N. 图中有M条边 (1 <= M <= 30,000) ,第j条边的长度为: d_ ...
- struts2练习时犯的错误(2016年11月4日)
1.Tomcat启动时报错 严重: 文档无效: 找不到语法. at (null:3:8) org.xml.sax.SAXParseException; systemId: file:/F:/Progr ...
- jsp中的js中获取项目路径的方法
在jsp中加上 <% String path = request.getContextPath(); String basePath = request.getScheme()+":/ ...
- Windows下CVSNT安装配置
首先要说明:X64下安装此软件会报 “cvs [login aborted]: WIN-4H9CRJO1TRA\Administrator: Switch to user failed due to ...
- python学习笔记__反射
反射 # 通过字符串的形式对对象中的成员进行操作(获取/查找/添加/删除). 操作的内置函数: 1.获取 getattr(object, name) # 去对象object中获取name的内容 c ...
- 如何在servlet的监听器中使用spring容器的bean
另外补充下:在web Server容器中,无论是Servlet,Filter,还是Listener都不是Spring容器管理的,因此我们都无法在这些类中直接使用Spring注解的方式来注入我们需要的对 ...
- 【题解】Tree-String Problem Codeforces 291E AC自动机
Prelude 传送到Codeforces:(/ω\)--- (/ω•\) Solution 很水的一道题. 对查询的串建出来AC自动机,然后树上随便跑跑就行了. 为什么要写这篇题解呢? 我第一眼看到 ...
- Gogs安装配置(快速搭建版)转载
gogs官网 oschina gogs介绍 一句话描述: 一款极易搭建的自助 Git 服务. 环境 centos7:golang+mysqldb+git 安装配置环境 yum install mysq ...