首页
学习
活动
专区
圈层
工具
发布
首页
学习
活动
专区
圈层
工具
MCP广场
社区首页 >专栏 >AT2005 [AGC003E] Sequential operations on Sequence

AT2005 [AGC003E] Sequential operations on Sequence

作者头像
yzxoi
发布2022-09-28 20:41:54
发布2022-09-28 20:41:54
58400
代码可运行
举报
文章被收录于专栏:OIOI
运行总次数:0
代码可运行

AT2005 [AGC003E] Sequential operations on Sequence

yzxoi

2022-08-21 (Updated: 2022-09-03)

oi

思维

题目链接:AGC003E

给定一个初始为 1\sim n 的序列,给定 q 次操作,每次操作把数组长度变为 q_i,新增/减少的数为上一个操作后的数组的重复。问 q 次操作后,1\sim n 每个数出现次数。

n,q \leq 10^5,q_i\leq 10^{18}

Tutorial

容易观察得到一个小性质:如果 q_i > q_{i+1}q_i 操作一定是无效的。

于是现在操作序列变为递增的。

考虑一次操作一定是把所有 q_i 对应的贡献先乘上 \lfloor \frac{q_i}{q_{i-1}} \rfloor,再处理 q_i \bmod q_{i-1} 部分的边角。

f_i 表示该次序列对最终序列的贡献,于是最后倒推一遍即可。

注意到边角部分一定会出现在 q_j < q_i \bmod q_{i-1},q_i \bmod q_{i-1} > q_{j+1}q_j 数量加入贡献,剩余部分继续拆分。

由于每次拆一定会 /2,于是差分一下,复杂度是 O(N\log ^2 N)

Solution

代码语言:javascript
代码运行次数:0
运行
复制
#include<bits/stdc++.h>
#define Tp template<typename Ty>
#define Ts template<typename Ty,typename... Ar>
#define W while
#define I inline
#define RI register int
#define LL long long
#define Cn const
#define CI Cn int&
#define gc getchar
#define D isdigit(c=gc())
#define pc(c) putchar((c))
using namespace std;
namespace Debug{
	Tp I void _debug(Cn char* f,Ty t){cerr<<f<<'='<<t<<endl;}
	Ts I void _debug(Cn char* f,Ty x,Ar... y){W(*f!=',') cerr<<*f++;cerr<<'='<<x<<",";_debug(f+1,y...);}
	Tp ostream& operator<<(ostream& os,Cn vector<Ty>& V){os<<"[";for(Cn auto& vv:V) os<<vv<<",";os<<"]";return os;}
	#define gdb(...) _debug(#__VA_ARGS__,__VA_ARGS__)
}using namespace Debug;
namespace FastIO{
	Tp I void read(Ty& x){char c;int f=1;x=0;W(!D) f=c^'-'?1:-1;W(x=(x<<3)+(x<<1)+(c&15),D);x*=f;}
	Ts I void read(Ty& x,Ar&... y){read(x),read(y...);}
	Tp I void write(Ty x){x<0&&(pc('-'),x=-x,0),x<10?(pc(x+'0'),0):(write(x/10),pc(x%10+'0'),0);}
	Tp I void writeln(Cn Ty& x){write(x),pc('\n');}
}using namespace FastIO;
Cn int N=1e5+10;
int n,q;LL a[N],f[N],g[N];
I void Sol(LL d,LL w){
	RI j=upper_bound(a+1,a+n+1,d)-a-1;
	if(!j) g[1]+=w,g[d+1]-=w;else f[j]+=d/a[j]*w,Sol(d%a[j],w);
}
int main(){
	RI i,p;LL x;for(read(p,q),a[n=1]=p,i=1;i<=q;i++){read(x);W(n&&a[n]>=x) n--;a[++n]=x;}
	for(f[n]=1,i=n-1;i>=1;i--) f[i]+=a[i+1]/a[i]*f[i+1],Sol(a[i+1]%a[i],f[i+1]);
	for(g[1]+=f[1],g[a[1]+1]-=f[1],i=1;i<=p;i++) writeln(g[i]+=g[i-1]);return 0;
}
本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
原始发表:2022-08-21 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • AT2005 [AGC003E] Sequential operations on Sequence
    • Tutorial
    • Solution
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档