Codeforces Round #740 (Div. 2) F. Top-Notch Insertions 线段树 / 平衡树 + 组合数学

传送门

文章目录

  • 题意:
  • 思路:

题意:

思路:

考虑最终的序列是什么鸭子的,首先序列肯定单调不降,也就是a1≤a2≤a3≤...≤ana_1\le a_2\le a_3\le ...\le a_na1a2a3...an,显然不可能都是≤\le号,因为如果插入的话是有可能产生<<<号的。假设我们现在有xxx<<< 号,那么应该有(n+n−1−xn)\binom{n+n-1-x}{n}(nn+n1x)个序列,因为我们将ai≤ai+1a_i\le a_{i+1}aiai+1转换成ai<ai+1+1a_i< a_i+1+1ai<ai+1+1,也就是扩展了一下值域,最多扩展了n−1−xn-1-xn1x个值域,也就是目前值域是[1,n∗2−1−x][1,n*2-1-x][1,n21x],从中选nnn个数构成一个严格递增的序列方案就是(n∗2−1−xn)\binom{n*2-1-x}{n}(nn21x),那么现在问题就是如何知道xxx

首先插入一组最多能产生一个<<<号,因为有可能若干数插到了同一个数后面,比如这个例子n=3,m=2,(2,1),(3,2)n=3,m=2,(2,1),(3,2)n=3,m=2,(2,1),(3,2),这样在a1a_1a1前面插入了两个数a2,a3a_2,a_3a2,a3,能确定a1>a2,a1>a3a_1>a_2,a_1>a_3a1>a2,a1>a3,但是不能确定a2a_2a2a3a_3a3的关系,所以这样的话我们就认为a2,a3a_2,a_3a2,a3之间是一个≤\le 号。我们考虑倒着来处理这mmm组,首先最终序列长度是nnn,假设现在有(xi,yi)(x_i,y_i)(xi,yi),那么先找到第yiy_iyi个位置ppp和第yi+1y_i+1yi+1的位置qqq,此时将ppp删除,将qqq标记一下代表有个数在他前面有<<<即可。至于为什么是yi,yi+1y_i,y_{i+1}yi,yi+1而不是yi−1,yiy_{i-1},y_{i}yi1,yi,是因为如果yi=1y_i=1yi=1的时候,yi−1=0y_i-1=0yi1=0,这个时候在线段树上二分就会出错。

在线段树上二分即可。

如果正着来,可以建一个基于权值的平衡树,让后模拟插入即可。

由于题目的问题,线段树不能每次建树。。。

// Problem: F. Top-Notch Insertions
// Contest: Codeforces - Codeforces Round #740 (Div. 2, based on VK Cup 2021 - Final (Engine))
// URL: https://codeforces.com/contest/1561/problem/F
// Memory Limit: 512 MB
// Time Limit: 3000 ms
// 
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#include<random>
#include<cassert>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid ((tr[u].l+tr[u].r)>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=1000010,mod=998244353,INF=0x3f3f3f3f;
const double eps=1e-6;int n,m;
PII p[N];
struct Node {int l,r;int cnt;
}tr[N<<2];void pushup(int u) {tr[u].cnt=tr[L].cnt+tr[R].cnt;
}void build(int u,int l,int r) {tr[u]={l,r};if(l==r) {tr[u].cnt=1;return;}build(L,l,Mid); build(R,Mid+1,r);pushup(u);
}void change(int u,int pos,int x) {if(tr[u].l==tr[u].r) {tr[u].cnt=x;return;}if(pos<=Mid) change(L,pos,x);else change(R,pos,x);pushup(u);
}int query(int u,int cnt) {if(tr[u].l==tr[u].r) return tr[u].l;if(cnt<=tr[L].cnt) return query(L,cnt);else return query(R,cnt-tr[L].cnt); 
}LL fun[N],inv[N];LL qmi(LL a,LL b) {LL ans=1;while(b) {if(b&1) ans=ans*a%mod;a=a*a%mod;b>>=1;}return ans%mod;
}LL C(int a,int b) {if(a<0||b<0||a<b) return 0;return fun[a]*inv[b]%mod*inv[a-b]%mod;
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);fun[0]=1;for(int i=1;i<N;i++) fun[i]=fun[i-1]*i%mod;inv[N-1]=qmi(fun[N-1],mod-2);for(int i=N-2;i>=0;i--) inv[i]=inv[i+1]*(i+1)%mod;build(1,1,N-1);int _; scanf("%d",&_);while(_--) {scanf("%d%d",&n,&m);for(int i=1;i<=m;i++) scanf("%d%d",&p[i].X,&p[i].Y);set<int>s; vector<int>v;for(int i=m;i>=1;i--) {int a=query(1,p[i].Y),b=query(1,p[i].Y+1);change(1,a,0); s.insert(b); v.pb(a);}for(auto x:v) change(1,x,1);int cnt=s.size();printf("%lld\n",C(n*2-1-cnt,n));}return 0;
}
/**/
// Problem: F. Top-Notch Insertions
// Contest: Codeforces - Codeforces Round #740 (Div. 2, based on VK Cup 2021 - Final (Engine))
// URL: https://codeforces.com/contest/1561/problem/F
// Memory Limit: 512 MB
// Time Limit: 3000 ms
// 
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#include<random>
#include<cassert>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid ((tr[u].l+tr[u].r)>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=1000010,mod=998244353,INF=0x3f3f3f3f;
const double eps=1e-6;int n,m;
int root,tot;
int x,y;
struct Node {int l,r;int size,fa,val,rank,tag;
}tr[N];int newnode(int v) {tr[++tot].val=v,tr[tot].size=1,tr[tot].rank=rand();tr[tot].l=tr[tot].r=tr[tot].tag=0;return tot;
}void pushup(int u) {if(!u) return;tr[u].size=tr[tr[u].l].size+tr[tr[u].r].size+1;// tr[u].val=max(tr[tr[u].l].val,tr[tr[u].r].val);
}void pushdown(int u) {if(!u) return;int tag=tr[u].tag; tr[u].tag=0;if(tr[u].l) tr[tr[u].l].tag+=tag,tr[tr[u].l].val+=tag;if(tr[u].r) tr[tr[u].r].tag+=tag,tr[tr[u].r].val+=tag;
}void split(int u,int k,int &x,int &y) {if(!u) { x=y=0; return; }pushdown(u);if(tr[u].val<=k) x=u,split(tr[u].r,k,tr[u].r,y);else y=u,split(tr[u].l,k,x,tr[u].l);pushup(u);
}int merge(int u,int v) {if(!u||!v) return u+v;if(tr[u].rank<tr[v].rank) {pushdown(u);tr[u].r=merge(tr[u].r,v);pushup(u);return u;}else {pushdown(v);tr[v].l=merge(u,tr[v].l);pushup(v);return v;}
}int find(int u,int x) {if(!u) return 0;if(tr[u].val==x) return 1;pushdown(u);if(x<tr[u].val) return find(tr[u].l,x);else return find(tr[u].r,x);
}LL fun[N],inv[N];LL qmi(LL a,LL b) {LL ans=1;while(b) {if(b&1) ans=ans*a%mod;a=a*a%mod;b>>=1;}return ans%mod;
}LL C(int a,int b) {if(a<0||b<0||a<b) return 0;return fun[a]*inv[b]%mod*inv[a-b]%mod;
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);fun[0]=1;for(int i=1;i<N;i++) fun[i]=fun[i-1]*i%mod;inv[N-1]=qmi(fun[N-1],mod-2);for(int i=N-2;i>=0;i--) inv[i]=inv[i+1]*(i+1)%mod;int _; scanf("%d",&_);while(_--) {root=tot=0;scanf("%d%d",&n,&m);for(int i=1;i<=m;i++) {int a,b; scanf("%d%d",&a,&b);int flag=find(root,b);split(root,b-1,x,y);if(y) tr[y].val++,tr[y].tag++;root=merge(x,flag? y:merge(newnode(b+1),y));}printf("%lld\n",C(n*2-1-tot,n));}return 0;
}
/**/

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314546.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

沃尔玛招聘.NET软件工程师

做为码农, 你可能会知道BAT, 微软, google, oracle, facebook等公司, 然而你知道沃尔玛吗? 在最新公布的2019年世界500强里(http://www.fortunechina.com/fortune500/c/2019-07/22/content_339535.htm?沃尔玛领先阿里巴巴181名.百度不是世界500强沃尔玛领先T公司236名沃尔玛领…

E:Three Blocks Palindrome(hard and easy)(树状数组 ? 前缀和?)

Three Blocks Palindrome (hard version) 思路 考虑到每个数字的范围是12001 ~ 2001 200&#xff0c;于是我们可以通过枚举两侧的元素来寻找最优答案。 我们有一个贪心策略&#xff0c;两侧都以我们枚举的元素作为结尾点&#xff0c;假如我们当前枚举的数字是1&#xff0c;于…

Codeforces Round #731 (Div. 3) G. How Many Paths? dfs + 拓扑 + 思维

传送门 题意&#xff1a; 给你一张nnn个点mmm条边的图&#xff0c;让你对每个点确定一个编号&#xff0c;规则如下&#xff1a; (1)(1)(1) 对于不能到的点编号为000。 (2)(2)(2) 对于只有一条路径能到这个点的点编号为111。 (3)(3)(3) 对于有不止一条路径能到这个点的点编号为…

.NET Core 3.0深入源码理解HttpClientFactory之实战

写在前面前面两篇文章透过源码角度&#xff0c;理解了HttpClientFactory的内部实现&#xff0c;当我们在项目中使用时&#xff0c;总会涉及以下几个问题&#xff1a;HttpClient超时处理以及重试机制HttpClient熔断器模式的实现HttpClient日志记录与追踪链接下来我们将从使用角度…

F:Maximum White Subtree(树形dp)

Maximum White Subtree 思路 如果考虑其覆盖范围只会到其子树上&#xff0c;不会到其父节点上的话(假设的情况)&#xff0c;这道题就非常好写了&#xff0c;就是一个简单的自底向上传递的树形dpdpdp。所以我们还要考虑的就是连接其父节点&#xff0c;因此我们只需要再进行一个…

树莓派 + Windows IoT Core 搭建环境监控系统

前言&#xff1a;Windows IoT 是微软为嵌入式开发板设计的一种物联网操作系统&#xff0c;运行Windows UWP&#xff08;C# 开发&#xff09;&#xff0c;可以设计出丰富的交互界面&#xff0c;驱动GPIO&#xff0c;连接一些传感器做有意思的事&#xff0c;本文详细介绍如何使用…

CF1547F Array Stabilization (GCD version) st表 + 尺取/二分

传送门 题意&#xff1a; 思路&#xff1a; 容易发现&#xff0c;我们将所有aaa都除上所有aaa的gcdgcdgcd&#xff0c;实际上就是让你求一个最小的lenlenlen&#xff0c;对于所有iii&#xff0c;gcd(ai,ai1,...,ailen−1)1gcd(a_i,a_{i1},...,a_{ilen-1})1gcd(ai​,ai1​,...…

Gauss高斯消元——模板

就是线性代数的初等行变化&#xff1a; 倍加。倍乘。交换行。 #include <bits/stdc.h> #define mp make_pair #define pb push_backusing namespace std;typedef long long ll; typedef pair<int, int> pii; typedef unsigned long long ull;const double pi ac…

.net持续集成测试篇之Nunit文件断言、字符串断言及集合断言

使用前面讲过的方法基本上能够完成工作中的大部分任务了,然而有些功能实现起来还是比较麻烦的,比如说字符串相等性比较不区分大小写,字符串是否匹配某一正则规则,集合中的每一个(某一个)元素是否符合特定规则等,Nunit提供了一些特定的方法用来实现一些用普通断言比较难以实现的…

Educational Codeforces Round 90 (Rated for Div. 2)(A, B, C, D, E)

Educational Codeforces Round 90 (Rated for Div. 2) Donut Shops 思路 分三种情况&#xff1a; ac/ba c / bac/b这个时候两个的单价是相同的&#xff0c;如果b1b 1b1&#xff0c;也就是aca cac&#xff0c;无论买多少数量的东西&#xff0c;这两个的价格都是一样的&…

Technocup 2020 - Elimination Round 2 E. Rock Is Push dp

传送门 文章目录题意&#xff1a;思路题意&#xff1a; 给你一个nmnmnm的方格&#xff0c;你初始在(1,1)(1,1)(1,1)点&#xff0c;有些位置有箱子&#xff0c;你可以走到某个位置向你的方向推动这个箱子&#xff0c;箱子不能出界&#xff0c;问你走到(n,m)(n,m)(n,m)有多少种方…

HttpClientFactory 使用说明 及 对 HttpClient 的回顾和对比

在 C# 中&#xff0c;平时我们在使用 HttpClient 的时候&#xff0c;会将 HttpClient 包裹在 using 内部进行声明和初始化&#xff0c;如&#xff1a;using(var httpClient new HttpClient()){ }至于为什么&#xff1f;无外乎是&#xff1a;项目代码中就是这样写的&#xf…

Codeforces Round #653 (Div. 3)(A, B, C, D, E1详解)

Codeforces Round #653 (Div. 3) Required Remainder Thinking(binary search) 既然是找最大值问题&#xff0c;我又懒得去推式子&#xff0c;于是我直接就上了一个二分&#xff0c;二分写法比结论稍微繁琐了一点吧&#xff0c;但是还是挺好想的。 根据题意&#xff0c;我们…

记录用友T+接口对接的心酸历程

前言&#xff1a;公司的业务主要是对接财务系统做单据传输或者凭证处理的&#xff0c;难免少不了和各大财务软件做数据对接&#xff0c;其中当然是必须通过接口来传递数据了。于是乎&#xff0c;用友T的版本来了&#xff0c;对接的工作自然是我来做&#xff0c;可没想到就是这样…

P6327 区间加区间sin和 线段树 + 数学

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个长度为nnn的序列aaa&#xff0c;有mmm次操作&#xff0c;每次操作分两种类型&#xff1a; (1)l,r,v(1)l,r,v(1)l,r,v&#xff0c;将al,al1,...,ara_l,a_{l1},...,a_ral​,al1​,...,ar​分别加上vvv…

集群故障处理之处理思路以及听诊三板斧(三十四)

前言 本篇主要分享一些处理故障和问题绝招&#xff0c;比如听诊三板斧&#xff1a;1&#xff09;查看日志2&#xff09;查看资源详情和事件3&#xff09;查看资源配置&#xff08;YAML&#xff09;如果还是不太好分析&#xff0c;那就祭出神器——kubectl-debug。最后&…

AC Automaton

简单复习一下字符串&#xff0c;顺便存个板子。 // Problem: P3808 【模板】AC自动机&#xff08;简单版&#xff09; // Contest: Luogu // URL: https://www.luogu.com.cn/problem/P3808 // Memory Limit: 512 MB // Time Limit: 1000 ms // // Powered by CP Editor (http…

selenium搜狗搜图简单操作(爬取任意关键字的图片)

Picture preview Steps 这里以搜狗搜图作为一个简单的例子&#xff1a; 第一步就是获取我们的需求了&#xff0c;输入需要下载的图片名字以及需要下载的图片数量。 接下来我们需要在搜狗搜图中输入用户的需求&#xff0c;透过xpathxpathxpath锁定搜索框&#xff0c;然后输入…

Microsoft REST API指南

经过3个月的碎片时间的翻译和校验&#xff0c;由长沙.NET技术社区翻译的英文原文文档《Microsoft REST API指南 》已经翻译完成&#xff0c;现刊载前十一章如下&#xff0c;欢迎大家点击“查看原文”按钮&#xff0c;查看指南的完整内容。PS&#xff1a;内容很长&#xff0c;全…

D. Salary Changing(二分,前缀和,贪心,瞎搞)

Salary Changing Thinking 这道题第一思路就是二分&#xff0c;模拟了一下样例&#xff0c;感觉好像行于是就开始写。 对于二分&#xff0c;我们一定是二分中位数是什么&#xff0c;二分的边界对我们来说是非常重要的&#xff0c;所以我们在二分前有必要确认我们的二分边界&…