洛谷noip 模拟赛 day1 T3
T7983 大芳的逆行板载
题目背景
大芳有一个不太好的习惯:在车里养青蛙。青蛙在一个n厘米(11n毫米s)的Van♂杆子上跳来跳去。她时常盯着青蛙看,以至于突然逆行不得不开始躲交叉弹。有一天他突发奇想,在杆子上每1厘米为一个单位,瞎涂上了墨水,并且使用mOgic,使青蛙跳过之处墨水浓度增加x。当然,他还会闲着无聊滴几滴墨水再涂♂抹均匀。
他现在无时无刻都想知道,第l厘米到第r厘米墨水的浓度是多少?
哦不!等等,他现在找到了一个计算器,可以输入几个数字与x,计算他们的x次幂和,所以。。。他想知道的是第l厘米到第r厘米墨水的浓度的x次幂和是多少?
题目描述
大芳有3种舰长技能骚操作
续:把青蛙放到第l厘米处,戳青蛙使其跳至r。效果:第l厘米至第r厘米墨水浓度增加x
- 抚♂摸:擦干杆子某一部分,重新滴加墨水并抹匀。效果:使第l厘米至第r厘米墨水浓度都变成x
最后一种是:
- 压线逆行,将车流看做⑨弹幕找安定点,掏出计算器,大喊板载后计算:
第l厘米至第r厘米墨水浓度的x次幂和是几何?记得答案要
模100000000710000000071000000007
输入输出格式
输入格式:
第一行nnn和mmm,表示杆子长n厘米,大芳要进行m次骚操作。
第二行nnn个数字,表示初始墨水浓度。第i个数字为第i厘米墨水浓度
接下来每行4个数字,依次为:操作编号(1、2或3),lll,rrr,xxx
输出格式:
每次进行3操作,输出一行表示答案
记得膜模1000000007
输入输出样例
5 5 19844 14611 26475 4488 6967 2 1 3 15627 2 1 2 30113 2 3 5 14686 2 5 5 32623 3 1 2 8
466266421
说明
kkk表示询问的幂的大小,也就是操作3对应的xxx。
对于20%的数据,满足n,m≤1000n,m\leq 1000n,m≤1000
对于另外20%的数据,满足k≤1k\leq 1k≤1
对于另外20%的数据,满足k≤2k\leq 2k≤2
对于另外20%的数据,满足n,m≤50000n,m\leq 50000n,m≤50000
对于100%的数据,满足n,m≤100000,0≤k≤10n,m\leq 100000,0\leq k \leq 10n,m≤100000,0≤k≤10
操作1,2对应的x≤109+7x\le 10^9+7x≤109+7
比赛挂了.......闲着没事放下T2 来艹这个 T3
因为想到了 二项式定理 (x+c)^n=sigma i=0..n (c^i)*C(n,i)*(x^(n-i))
然后就觉得可能可以写了...
果然代码能力不足写挂了 赛后调了一个多小时发现时标记之间的关系弄错了............惨啊
——————————————————————————————————————————
这道题呢 涉及到区间加 区间覆盖 区间幂求和
区间覆盖其实很好弄
主要在于区间加 但是他求的幂只有(1——10) 所以我们可以暴力维护一波(1——10)的幂 -用二项式定理
这个不懂的百度学学咯 复杂度也就只有 100logn 能接受
当然推的时候注意先推幂的次数比较大的 其他都是正常的线段树操作了
当然别忘了覆盖标记遇到加法标记要把加法标记清零啊
#include<cstdio> #include<cstring> #include<algorithm> #include<cmath> #define LL long long using namespace std; const int M=1<<18,mod=1000000007; LL read(){LL ans=0,f=1,c=getchar();while(c<'0'||c>'9'){if(c=='-') f=-1; c=getchar();}while(c>='0'&&c<='9'){ans=ans*10+(c-'0'); c=getchar();}return ans*f; } LL C[15][15],c[15],q[15],num[M],w; int h,n,m,L,R; struct node{LL f,tag;LL m[11]; }tr[M]; void calc_f(int x,int l,int r,LL v){LL sz=r-l+1;tr[x].m[1]=sz*v%mod;for(int i=2;i<=10;i++) tr[x].m[i]=(tr[x].m[i-1]*v)%mod; } void push(int x,int k,LL v){LL sum=1;for(int i=1;i<=k;i++) sum=sum*v%mod,tr[x].m[k]=(tr[x].m[k]+sum*C[k][i]%mod*tr[x].m[k-i]%mod)%mod; } void calc_tag(int x,LL v){for(int i=10;i>=1;i--) push(x,i,v); } void up(int x){int l=x<<1,r=x<<1^1;for(int i=0;i<=10;i++) tr[x].m[i]=(tr[l].m[i]+tr[r].m[i])%mod; } void down(int x,int l,int r){int ll=x<<1,rr=x<<1^1;LL mid=(l+r)>>1;if(tr[x].f!=-1){tr[ll].f=tr[x].f; tr[rr].f=tr[x].f;calc_f(ll,l,mid,tr[x].f);calc_f(rr,mid+1,r,tr[x].f);tr[x].f=-1; tr[ll].tag=tr[rr].tag=0;}if(tr[x].tag){if(l==r) return ;calc_tag(ll,tr[x].tag);calc_tag(rr,tr[x].tag);tr[x<<1].tag+=tr[x].tag; tr[x<<1^1].tag+=tr[x].tag; tr[x].tag=0;} } void modify_tag(int x,int l,int r){if(L<=l&&r<=R){tr[x].tag+=w;if(l==r){tr[x].m[1]+=w;for(int i=2;i<=10;i++) tr[x].m[i]=tr[x].m[i-1]*tr[x].m[1]%mod;}else calc_tag(x,w);return ;}down(x,l,r);int mid=(l+r)>>1;if(L<=mid) modify_tag(x<<1,l,mid);if(R>mid) modify_tag(x<<1^1,mid+1,r);up(x); } void modify_f(int x,int l,int r){if(L<=l&&r<=R){tr[x].f=w;tr[x].tag=0;calc_f(x,l,r,w);return ;}down(x,l,r);int mid=(l+r)>>1;if(L<=mid) modify_f(x<<1,l,mid);if(R>mid) modify_f(x<<1^1,mid+1,r);up(x); } LL query(int x,int l,int r,int k){if(L<=l&&r<=R) return tr[x].m[k];down(x,l,r);LL mid=(l+r)>>1;LL sum=0;if(L<=mid) sum=sum+query(x<<1,l,mid,k);if(R>mid) sum=sum+query(x<<1^1,mid+1,r,k);return sum%mod; } void build(int x,int l,int r){tr[x].f=-1;if(l==r){tr[x].m[0]=1;tr[x].m[1]=num[l]; for(int i=2;i<=10;i++) tr[x].m[i]=(tr[x].m[i-1]*tr[x].m[1])%mod; return ;}int mid=(l+r)>>1;build(x<<1,l,mid);build(x<<1^1,mid+1,r);up(x); } void prepare(){q[0]=1; q[1]=1;for(int i=2;i<=10;i++) q[i]=q[i-1]*i;for(int i=1;i<=10;i++)for(int j=1;j<=i;j++)C[i][j]=q[i]/(q[j]*q[i-j]); } int main() {prepare();n=read(); m=read();for(int i=1;i<=n;i++) num[i]=read();build(1,1,n); for(int i=1;i<=m;i++){h=read();if(h==1){L=read(); R=read(); w=read();modify_tag(1,1,n);}else if(h==2){L=read(); R=read(); w=read();modify_f(1,1,n);}else{L=read(); R=read(); int k=read();printf("%lld\n",query(1,1,n,k));}//printf("[%lld]\n",tr[1].m[1]); }return 0; }
View Code
转载于:https://www.cnblogs.com/lyzuikeai/p/7263132.html
洛谷noip 模拟赛 day1 T3相关推荐
- 洛谷mNOIP模拟赛Day2-入阵曲
题目背景 pdf题面和大样例链接:http://pan.baidu.com/s/1cawM7c 密码:xgxv 丹青千秋酿,一醉解愁肠. 无悔少年枉,只愿壮志狂. 题目描述 小 F 很喜欢数学,但是到 ...
- 【HHHOJ】NOIP模拟赛 捌 解题报告
点此进入比赛 得分: \(30+30+70=130\)(弱爆了) 排名: \(Rank\ 22\) \(Rating\):\(-31\) \(T1\):[HHHOJ260]「NOIP模拟赛 捌」Dig ...
- NOIP模拟赛csy2021/10/30
NOIP模拟赛csy2021/10/30 比赛时间规划 赛后反思与总结 这..总的来说感觉打的很不好,根本没有状态,有一部分原因是今天来晚了,太慌,更多的还是这次题感觉很难o(╥﹏╥)o 比赛时间规划 ...
- NOI.AC NOIP模拟赛 第六场 游记
NOI.AC NOIP模拟赛 第六场 游记 queen 题目大意: 在一个\(n\times n(n\le10^5)\)的棋盘上,放有\(m(m\le10^5)\)个皇后,其中每一个皇后都可以向上.下 ...
- 【noip模拟赛4】Matrix67的派对 暴力dfs
[noip模拟赛4]Matrix67的派对 描述 Matrix67发现身高接近的人似乎更合得来.Matrix67举办的派对共有N(1<=N<=10)个人参加,Matrix67需要把他们安排 ...
- jyzy noip模拟赛5.22-2
不知道哪来的题 jyzy noip模拟赛5.22-2 样例输入 1 2 3 4 样例输出 0.200000000000000 数据 |a|,|b|,|c|,|d|<=1e9 很多大佬迅速想到二分 ...
- PKUSC 模拟赛 day1 下午总结
下午到了机房之后又困又饿,还要被强行摁着看英文题,简直差评 第一题是NOIP模拟赛的原题,随便模拟就好啦 本人模拟功力太渣不小心打错了个变量,居然调了40多分钟QAQ #include<cstd ...
- 【WZOI第二次NOIP模拟赛Day1T2】世界末日 解题报告
[WZOI第二次NOIP模拟赛Day1T2]世界末日 Problem 2 世界末日 (doomsday.pas/c/cpp) 背景 话说CWQ大牛终于打开了那扇神秘大门,但迎接他的不是什么神秘的东西, ...
- NOIP 模拟赛 长寿花 题解
NOIP 模拟赛 长寿花 题解 要放 \(n\) 层物品,第 \(i\) 层有 \(a_i\) 个位置放物品,物品有 \(m\) 中颜色,有约束条件: 同一层两个相邻物品颜色不能相同. 相邻两层颜色集 ...
最新文章
- Word2016如何取消宏已被禁用的提示
- Python基础教程:字典和Json
- 1、MySQL 8.0.20最新版本在Linux上安装
- u-boot.lds文件简介
- 紧急整理了 20 道 Spring Boot 面试题,我经常拿来面试别人
- 银行营业网点管理系统——dao包(BranchesDao)
- 对怀孕的人有害的食物。。。朋友们记住咯!(欢迎转载)
- 站在BERT肩膀上的NLP新秀们(PART II)
- python3爬取网易云歌单数据清洗_如何利用Python网络爬虫爬取网易云音乐歌词
- eclipse中提示“java application configureation name i”
- Linux 设备驱动的固件加载【转】
- 数字电路:数据选择器与译码器
- 前端标准模式与混杂模式的区别
- git 删除历史commit
- 创意几何点线面三角形孟菲斯海报banneer广告设计ai矢量背景素材
- 全概率公式与贝叶斯公式-机器学习
- python123.io同一个号可以同时在不同设备登陆吗-Python 高级 3
- 自我保健很重要:先付钱
- SAP中通过定制开发在寄售结算MRKO中实现发票校验功能
- vue移动端滑动切换图片的一个简单思路