牛客网CSP-S提高组赛前集训营Round4
牛客网CSP-S提高组赛前集训营
标签(空格分隔): 题解 算法 模拟赛
| 题目 | 描述 | 做法 |
|---|---|---|
| \(BSOJ6377\) | 求由\(n\)长度的数组复制\(k\)次的数组里每个连续子序列出现数字种类的和 | 对每一种颜色计算贡献区间 |
| \(BSOJ6378\) | 求树上两条长度分别为\(p\),\(q\)不相交路径对数 | 换根背包\(dp~or~\)\(LCA\)+差分 |
| \(BSOJ6379\) |
\(T1\)
有一个长为\(n×k\)的数组,它是由长为\(n\)的数组\(A_1\),\(A_2\),...,\(A_n\)重复k次得到的。
定义这个数组的一个区间的权值为它里面不同的数的个数,现在,你需要求出对于这个数组的每个非空区间的权值之和。
答案对\(10^9+7\)取模。
注意到计算每一个区间的影响是很难的,因为我们能表示一个区间颜色种类数的方法是最快\(O(lenlog_2len)\)的
而这道题\(len\le 10^{14}\)让我们放弃
这又是一个套路了,既然我们不好直接求,那么拆分问题
给一个相似的例子
求\(a_1,a_2\cdots,a_n\)两两异或的和
肯定正着做不好做,但考虑到,这个问题在\(a_i\in\{0,2^k\}\)时就很好做
因此想到拆位
对于这道题也是一样,考虑计算每种颜色贡献的区间
仿佛也不好做,但正难则反,计算每种颜色不会贡献的区间很简单
那就是相邻的两个颜色间的所有区间
对\(n\)的序列我们可以轻易用一次扫描\(O(n)\)通过记录每个颜色上次出现的位置\(last_i\)利用组合数求出
但对\(n*k\)的呢

考虑有些部分重复求了
红色的被夹在区间中间的部分有\(K\)个
绿色蓝色各有\(1\)个
绿蓝拼在一起的有\(K-1\)个
因此我们只需要求\(n/2n\)的区间把色块拼起来就可以了
最后说明一下
长度\(\le len\)区间个数\(\frac{len(len+1)}{2}\)
#include<bits/stdc++.h>
#define re register
#define N 100001
#define INF 0x3f3f3f3f
#define mod 1000000007ll
using namespace std;
inline char nc(){
static char buf[1048576],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,1048576,stdin),p1==p2)?EOF:*p1++;
}
#define getchar nc
template<typename _int>
inline void read(re _int& x){
re char opt;re _int flag=1,res=0;
while((opt=getchar())<'0'||opt>'9')if(opt=='-')flag=-1;
while(opt>='0'&&opt<='9'){res=(res<<3)+(res<<1)+opt-'0';opt=getchar();}
x=res*flag;
}
typedef long long ll;
int a[N<<1],pos[N],t[N],tot,l[N<<1];
ll ans,n,k;
inline ll f(re ll x){return x*(x+1)%mod*500000004ll%mod;}
inline void Read(void){
re int i;read(n);read(k);
for(i=1;i<=n;++i){read(a[i]);t[i]=a[i];}
sort(t+1,t+n+1);tot=unique(t+1,t+n+1)-(t+1);
for(i=1;i<=n;++i){a[i]=lower_bound(t+1,t+tot+1,a[i])-t;a[i+n]=a[i];}
}
inline void Solve(void){
re int i;re ll self,sum,len,left,right;
ans=f(1ll*n*k%mod)*tot%mod;
self=0;
for(i=1;i<=n;++i){
l[i]=pos[a[i]];pos[a[i]]=i;
if(l[i]+1==i||!l[i])continue;len=1ll*i-l[i]-1ll;
self=(self+f(len))%mod;
}
for(i=1;i<=tot;++i)pos[i]=0;
left=0;
for(i=1;i<=n;++i){
l[i]=pos[a[i]];pos[a[i]]=i;
if(l[i]||l[i]+1==i)continue;len=1ll*i-l[i]-1ll;
left=(left+(1ll*len*(len+1ll)>>1ll))%mod;
}
for(i=1;i<=tot;++i)pos[i]=n+1;
right=0;
for(i=n;i>=1;--i){
l[i]=pos[a[i]];pos[a[i]]=i;
if(l[i]!=n+1||l[i]-1==i)continue;len=1ll*l[i]-i-1ll;
right=(right+(1ll*len*(len+1ll)>>1ll))%mod;
}
for(i=1;i<=tot;++i)pos[i]=0;
sum=0;
for(i=1;i<=(n<<1);++i){
l[i]=pos[a[i]];pos[a[i]]=i;
if(l[i]+1==i)continue;len=1ll*i-l[i]-1ll;
sum=(sum+(1ll*len*(len+1ll)>>1ll))%mod;
}
sum=((sum-self*2ll%mod-left)%mod+mod)%mod;
ans=(((ans-1ll*k*self%mod+mod)%mod-sum*(k-1ll)%mod-left-right)%mod+mod)%mod;
printf("%lld\n",(ans%mod+mod)%mod);
}
int main(void){
Read();Solve();
return 0;
}
\(T2\)
\(有一棵n个点的树和两个整数p, q,求满足以下条件的四元组(a, b, c, d)的个数:\)
\(1.1≤a, b, c, d≤n。\)
\(2.点a到点b的经过的边数为p。\)
\(3.点c到点d的经过的边数为q。\)
\(4.不存在一个点,它既在点a到点b的路径上,又在点c到点d的路径上。\)
注意到不相交的情况太多(我考场上写换根写疯了)
首先考虑反求问题,不相交的路径数=所有路径数-相交路径数
考虑相交有哪些情况

我们把第三种颜色互换一下就变成了第二种
我们发现绿色的链总是过黑链的\(lca\)的
我们考虑枚举这个\(lca\)
他可以是两条只在子树内的链或者是一条子树内的链和一条从子树内(可以不进)到子树外的链
因此我们只需要算过一个点\(x\)子树内的链长为\(p/q\)方案数\(f_{x,p/q}\)以及从子树内(可以不进)到子树外的链长为\(p/q\)方案数\(g_{x,p/q}\)
考虑一条链拆成两条
设从\(x\)出发,子数内单链长为\(k\)的方案\(sf_{x,k}\)
则\(\displaystyle{sf_{x,k}=\sum_{y\in son_x}}sf_{y,k-1}\)
设从\(x\)出发,子数外单链长为\(k\)的方案\(gf_{x,k}\)
则\(\displaystyle{sg_{x,k}=sg_{fa_x,k-1}+sf_{fa_x,k-1}-sf_{x,k-2}}\)
就是来自父亲上面的,父亲下面不过自己的
合并的时候是类似背包的,我们把前面的儿子(\(sf'\))与\(y(sf)\)合并
因此\(\displaystyle{f_{x,p/q}=\sum_{k\in[0,p/q)}sf'_{x,k}\times sf_{x,p/q-k-1}}\)
\(g\)更简单
\(\displaystyle{g_{x,p/q}=\sum_{k\in[0,p/q]}sf_{x,k}\times sg_{x,p/q-k}}\)
\(T3\)
留坑
牛客网CSP-S提高组赛前集训营Round4的更多相关文章
- 牛客CSP-S提高组赛前集训营1
牛客CSP-S提高组赛前集训营1 比赛链接 官方题解 before:T1观察+结论题,T2树形Dp,可以换根或up&down,T3正解妙,转化为图上问题.题目质量不错,但数据太水了~. A-仓 ...
- 牛客CSP-S提高组赛前集训营3
A 货物收集 显然是一个二分答案的题. #include<iostream> #include<cstdio> #include<cstring> #include ...
- 牛客CSP-S提高组赛前集训营3 赛后总结
货物收集 二分答案.复杂度\(O(n\log n)\). 货物分组 用费用提前计算的思想,考虑用一个新的箱子来装货物会发生什么. 显然费用会加上后面的所有货物的总重. \(60\)分的\(O(n^2) ...
- 牛客CSP-S提高组赛前集训营2 ———— 2019.10.31
比赛链接 期望得分:100+20+20 实际得分:40+20+30 awa cccc T1 :基于贪心的思路,然后开始爆搜(雾 那必然是会死的,好吧他就是死了 #include<iostrea ...
- 牛客CSP-S提高组赛前集训营1———2019.10.29 18:30 至 22:00
期望得分:100+0+10 实际得分:40+0+0 考炸了... T1:题目链接 究竟为什么会这样,,, 仔细研读我的丑代码 发现... 枯辽.... #include<cstdio> # ...
- 牛客CSP-S提高组赛前集训营2 T2沙漠点列
原题链接 算法不难,比赛的时候就和cyc大佬一起yy了正解,不过因为交的时候比较急(要回寝室惹),我有两数组开错大小直接爆到50,cyc大佬则只把文件输入关了一半,直接爆零(╯ ̄Д ̄)╯┻━┻ 要尽量 ...
- 20191029 牛客CSP-S提高组赛前集训营1
前一个小时看这几道题感觉要爆零 A. 仓鼠的石子游戏 分析一下发现a[i]>1a[i]>1a[i]>1时后先手必输,a[i]=1a[i]=1a[i]=1时先手必赢 然后直接看1的个数 ...
- 牛客CSP-S提高组赛前集训营4 赛后总结
复读数组 分成 3 种区间算答案: 一个块内的区间 两个块交界处,长度小于块长的区间 长度不小于块长的区间 对于第三种区间,容易发现每个区间的权值一样,只需要算出个数即可. 对于前两种空间,我的思路是 ...
- 牛客CSP-S提高组赛前集训营5 赛后总结
A.无形的博弈 心理题. 答案为\(2^n\),可感性理解结论的正确性. #include<bits/stdc++.h> #define LL long long const LL Mod ...
随机推荐
- 简易版php文件上传_超详细详解
1.php简化版的图片上传(没有各种验证) 1234 <form action="" enctype="multipart/form-data" meth ...
- memcpy() 实现循环缓冲区数据的读写
使用memcpy()函数做循环缓冲区的读写 首先对mencpy函数做个简单的介绍 下面是 memcpy() 函数的声明 void *memcpy(void *str1, const void *str ...
- IAR_STM32_CCM内存使用
在IAR中,硬件环境为STM32F464单片机,额外的CCM不能用DMA访问,但可以作为内部额外扩充RAM使用 修改文件中的内容为: define symbol __ICFEDIT_region_RO ...
- 【题解】【网络流24题】汽车加油行驶问题 [P4009] [Loj6223]
[题解][网络流24题]汽车加油行驶问题 [P4009] [Loj6223] 传送门:汽车加油行驶问题 \([P4009]\) \([Loj6223]\) [题目描述] 给出一个 \(N \times ...
- 类、闭包、monand
类.闭包.monand 三者的相同点都是封装了数据与操作. 类:定义数据类型,倾向于数据的抽象:集合数据为操作提供上下文. 闭包:对环境变量和操作的封装抽象.操作基于环境变量,核心问题是环境变量的管理 ...
- N(C)O(S)I(P)P 2019 退役记
N(C)O(S)I(P)P 2019 退役记 day-4 今天下午老师突然咕了,于是一下午欢乐时光 今天上午考试T3线段树维护个区间加,区间乘 一遍过编译,一遍过样例(第一次,俺比较弱(虽然也发现和暴 ...
- 一.B/S架构和C/S架构
1.B/S架构 Browser-Server, 浏览器和服务器架构.包含客户端浏览器.web应用服务器.数据库服务器的软件系统.用户只需要一个浏览器就可以访问服务.系统更新的时候,只需要更新服务端, ...
- PyTorch 安装 报错,原因是pip 不是64位的。
原因: import pip._internal print(pip._internal.pep425tags.get_supported()) 换位64位的python版本. import pip. ...
- Java 之 Set 接口
一.Set 概述 java.util.Set 接口继承 collection 接口,它与 Collection 接口中的方法基本一致,并没有对 Collection 接口进行功能上的扩充,只是比 Co ...
- Redux 进阶之 react-redux 和 redux-thunk 的应用
1. react-redux React-Redux 是 Redux 的官方 React 绑定库. React-Redux 能够使你的React组件从Redux store中读取数据,并且向 stor ...