题解 Dove 打扑克
考场上觉得复杂度是假的就没怎么优化,然后考完题解帮我证明了它是真的……
首先合并可以用并查集维护,可以顺便维护出集合的大小
对于操作2,发现如果 \(size_i\) 是确定的,可以用权值线段树很方便的维护出合法的 \(size_j\)的个数
每次只需枚举出现过的 \(size_i\) 即可,所以我觉得复杂度是假的
- 如果存在 \(\sum a_i = n\) ,那么有 \(diff \{a_i\} \leqslant \sqrt n\)
我们有 \(\sum size_i = n\) ,所以不同的size个数只有不超过根号个
所以用个unordered_set维护当前存在的size,线段树查询即可
yysy,我把线段树换成树状数组从1700ms变成了400ms
Code:
#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define N 100010
#define ll long long
#define reg register int
//#define int long long
char buf[1<<21], *p1=buf, *p2=buf;
#define getchar() (p1==p2&&(p2=(p1=buf)+fread(buf, 1, 1<<21, stdin)), p1==p2?EOF:*p1++)
inline int read() {
int ans=0, f=1; char c=getchar();
while (!isdigit(c)) {if (c=='-') f=-f; c=getchar();}
while (isdigit(c)) {ans=(ans<<3)+(ans<<1)+(c^48); c=getchar();}
return ans*f;
}
int n, m;
int fa[N], cnt[N];
bool vis[N];
inline int find(int p) {return fa[p]==p?p:fa[p]=find(fa[p]);}
namespace force{
void solve() {
for (int i=1; i<=n; ++i) fa[i]=i, cnt[i]=1, vis[i]=1;
ll ans;
for (int i=1,x,y,c,f1,f2; i<=m; ++i) {
if (read()&1) {
x=read(); y=read();
f1=find(x), f2=find(y);
if (f1!=f2) {
cnt[f1]+=cnt[f2];
vis[f2]=0;
fa[f2]=f1;
}
}
else {
ans=0;
c=read();
for (int i=1; i<=n; ++i) if (vis[i])
for (int j=i+1; j<=n; ++j) if (vis[j])
if (abs(cnt[i]-cnt[j])>=c) ++ans;
printf("%lld\n", ans);
}
}
exit(0);
}
}
namespace task1{
int tl[N<<2], tr[N<<2], sum[N<<2];
#define tl(p) tl[p]
#define tr(p) tr[p]
#define sum(p) sum[p]
#define pushup(p) sum(p)=sum(p<<1)+sum(p<<1|1)
void build(int p, int l, int r) {
tl(p)=l; tr(p)=r;
if (l==r) return ;
int mid=(l+r)>>1;
build(p<<1, l, mid);
build(p<<1|1, mid+1, r);
}
void upd(int p, int pos, int dat) {
if (tl(p)==tr(p)) {sum(p)+=dat; return ;}
int mid=(tl(p)+tr(p))>>1;
if (pos<=mid) upd(p<<1, pos, dat);
else upd(p<<1|1, pos, dat);
pushup(p);
}
int query(int p, int l, int r) {
if (l<=tl(p) && r>=tr(p)) return sum(p);
int mid=(tl(p)+tr(p))>>1, ans=0;
if (l<=mid) ans+=query(p<<1, l, r);
if (r>mid) ans+=query(p<<1|1, l, r);
return ans;
}
void solve() {
build(1, 1, n); upd(1, 1, n);
for (int i=1; i<=n; ++i) fa[i]=i, cnt[i]=1, vis[i]=1;
ll ans;
for (int i=1,x,y,c,f1,f2; i<=m; ++i) {
if (read()&1) {
x=read(); y=read();
f1=find(x), f2=find(y);
if (f1!=f2) {
upd(1, cnt[f1], -1); upd(1, cnt[f2], -1);
cnt[f1]+=cnt[f2];
upd(1, cnt[f1], 1);
vis[f2]=0;
fa[f2]=f1;
}
}
else {
ans=0;
c=read();
for (reg i=1,l,r; i<=n; ++i) if (vis[i]) {
if ((l=cnt[i]-c)>=1) ans+=query(1, 1, l)-(c==0);
if ((r=cnt[i]+c+(c==0))<=n) ans+=query(1, r, n);
}
printf("%lld\n", ans/2);
}
}
exit(0);
}
}
namespace task2{
int tl[N<<2], tr[N<<2], sum[N<<2], tot[N], maxn=1;
#define tl(p) tl[p]
#define tr(p) tr[p]
#define sum(p) sum[p]
#define pushup(p) sum(p)=sum(p<<1)+sum(p<<1|1)
void build(int p, int l, int r) {
tl(p)=l; tr(p)=r;
if (l==r) return ;
int mid=(l+r)>>1;
build(p<<1, l, mid);
build(p<<1|1, mid+1, r);
}
void upd(int p, int pos, int dat) {
if (tl(p)==tr(p)) {sum(p)+=dat; return ;}
int mid=(tl(p)+tr(p))>>1;
if (pos<=mid) upd(p<<1, pos, dat);
else upd(p<<1|1, pos, dat);
pushup(p);
}
int query(int p, int l, int r) {
if (l<=tl(p) && r>=tr(p)) return sum(p);
int mid=(tl(p)+tr(p))>>1, ans=0;
if (l<=mid) ans+=query(p<<1, l, r);
if (r>mid) ans+=query(p<<1|1, l, r);
return ans;
}
void solve() {
build(1, 1, n); upd(1, 1, n); tot[1]=n;
for (int i=1; i<=n; ++i) fa[i]=i, cnt[i]=1, vis[i]=1;
ll ans;
for (int i=1,x,y,c,f1,f2; i<=m; ++i) {
if (read()&1) {
x=read(); y=read();
f1=find(x), f2=find(y);
if (f1!=f2) {
upd(1, cnt[f1], -1); upd(1, cnt[f2], -1);
--tot[cnt[f1]]; --tot[cnt[f2]];
cnt[f1]+=cnt[f2];
++tot[cnt[f1]]; maxn=max(maxn, cnt[f1]);
upd(1, cnt[f1], 1);
vis[f2]=0;
fa[f2]=f1;
}
}
else {
ans=0;
c=read();
for (reg i=1,l,r; i<=maxn; ++i) if (tot[i]) {
if ((l=i-c)>=1) ans+=1ll*tot[i]*(query(1, 1, l)-(c==0));
//cout<<"l: "<<l<<' '<<tot[i]*(query(1, 1, l)-(c==0))<<endl;
if ((r=i+c+(c==0))<=n) ans+=1ll*tot[i]*query(1, r, n);
//cout<<"r: "<<r<<' '<<query(1, r, n)<<endl;
}
printf("%lld\n", ans/2);
}
}
exit(0);
}
}
namespace task{
unordered_set<int> s;
unordered_set<int>::iterator sta[N];
int sum[N], tot[N], top;
inline void upd(int i, int dat) {for (; i<=n; i+=i&-i) sum[i]+=dat;}
inline int query(int i) {int ans=0; for (; i; i-=i&-i) ans+=sum[i]; return ans;}
void solve() {
upd(1, n); tot[1]=n; s.insert(1);
for (int i=1; i<=n; ++i) fa[i]=i, cnt[i]=1, vis[i]=1;
ll ans;
for (int i=1,x,y,c,f1,f2; i<=m; ++i) {
if (read()&1) {
x=read(); y=read();
f1=find(x), f2=find(y);
if (cnt[f1]<cnt[f2]) swap(f1, f2);
if (f1!=f2) {
upd(cnt[f1], -1); upd(cnt[f2], -1);
--tot[cnt[f1]]; --tot[cnt[f2]];
cnt[f1]+=cnt[f2];
if (++tot[cnt[f1]]==1) s.insert(cnt[f1]);
upd(cnt[f1], 1);
vis[f2]=0;
fa[f2]=f1;
}
}
else {
ans=0;
c=read();
int l, r;
for (unordered_set<int>::iterator it=s.begin(); it!=s.end(); ++it) {
//cout<<"*it: "<<*it<<endl;
if (!tot[*it]) {sta[++top]=it; continue;}
if ((l=*it-c)>=1) ans+=1ll*tot[*it]*(query(l)-(c==0));
//cout<<"l: "<<l<<' '<<tot[i]*(query(1, 1, l)-(c==0))<<endl;
if ((r=*it+c+(c==0))<=n) ans+=1ll*tot[*it]*(query(n)-query(r-1));
//cout<<"r: "<<r<<' '<<query(1, r, n)<<endl;
}
while (top) s.erase(sta[top--]);
printf("%lld\n", ans/2);
}
}
exit(0);
}
}
signed main()
{
n=read(); m=read();
//if (n<=100) force::solve();
//else if (n<=1000) task1::solve();
//else task2::solve();
task::solve();
return 0;
}
题解 Dove 打扑克的更多相关文章
- 「10.28」Dove 打扑克(链表)·Cicada 与排序(概率)·Cicada 拿衣服(各种数据结构)
A. Dove 打扑克 考场思考半天线段树树状数组,没有什么想法 打完暴力后突然想到此题用链表实现会很快. 因为只有$n$堆,所以设最多有$x$个不同的堆数,那么$x\times (x-1)/2==n ...
- [UPC10525]:Dove打扑克(暴力+模拟)
题目描述 $Dove$和$Cicada$是好朋友,他们经常在一起打扑克来消遣时光,但是他们打的扑克有不同的玩法. 最开始时,牌桌上会有$n$个牌堆,每个牌堆有且仅有一张牌,第$i$个牌堆里里里那个扑克 ...
- NOIP 模拟 $36\; \rm Dove 打扑克$
题解 \(by\;zj\varphi\) 引理 对于一个和为 \(n\) 的数列,不同的数的个数最多为 \(\sqrt n\) 证明: 一个有 \(n\) 个不同的数的数列,和最小就是 \(n\) 的 ...
- 20210811 Dove 打扑克,Cicada 与排序,Cicada 拿衣服
考场 开考感觉 T3 比较可做.T1 看上去不难但毫无思路. 先想了 25min T3,想到一个确定左端点,二分最长的右端点,甚至想到了用猫树维护区间 or and...上厕所回来发现假了,没有单调性 ...
- 晚间测试13 A. Dove 打扑克 vector +模拟
题目描述 分析 这道题比较关键的一点就是要看出最终牌数的种类数不会超过 \(\sqrt{n}\) 种 知道了这个性质我们就可以用 \(vector\) 维护一个有序的序列 \(vector\) 中存放 ...
- Noip模拟36 2021.8.11
刚题的习惯还是改不了,怎么办??? T1 Dove打扑克 考场上打的动态开点线段树+并查集,考后发现自己像一个傻子,并查集就行.. 这几天恶补数据结构疯了 用树状数组维护后缀和,$siz_i$表示编号 ...
- csp-s模拟测试91
csp-s模拟测试91 倒悬吃屎的一套题. $T1$认真(?)分析题意发现复杂度不能带$n$(?),计划直接维护答案,考虑操作对答案的影响,未果.突然发现可以动态开点权值线段树打部分分,后来$Tm$一 ...
- noip模拟36
\(\color{white}{\mathbb{荷花映日,莲叶遮天,名之以:残荷}}\) 今天再次翻车掉出前十 开题看错 \(t1\) 以为操作2的值固定发现是个简单题,然后 \(t2\) 开始大力 ...
- 2021.8.11考试总结[NOIP模拟36]
T1 Dove玩扑克 考场并查集加树状数组加桶期望$65pts$实际$80pts$,考后多开个数组记哪些数出现过,只扫出现过的数就切了.用$set$维护可以把被删没的数去掉,更快. $code:$ 1 ...
随机推荐
- 由ctf来看java报错的危害
很多java报错在我们渗透的时候经常会被发现,但由于没什么用,危害比较低被忽略,开发也很不愿意修改. 但从纵深防御的角度来说,多个小问题的结合就会产生严重的问题.此次遇到的一个ctf题就是一个例子. ...
- UI自动化测试框架Gauge 碰到无法识别Undefined Steps 红色波纹标记
如果碰到无法识别的情况,例如下面的红色波纹,可以试一下: 第一步: 第二步: 不勾选'offline work' 第三部:刷新之后可以重新编译.
- PYTHON3.4.4 升级pip
python -m pip install --upgrade pipupgrade 前面是两个 "-"
- C语言:虚拟地址 和编译模式
所谓虚拟地址空间,就是程序可以使用的虚拟地址的有效范围.虚拟地址和物理地址的映射关系由操作系统决定,相应地,虚拟地址空间的大小也由操作系统决定,但还会受到编译模式的影响.这节我们先讲解CPU,再讲解编 ...
- 家庭账本开发day01
--好长时间没有进行web应用的开发,手有些生疏了,知识点也有些遗忘了,不过抹油关系,边开发边复习边学习新的东西. 今天主要完成了,家庭记账本的项目梳理,如下: 利用layUI模板 修改layUImi ...
- SLAM的数学基础(4):先验概率、后验概率、贝叶斯准则
假设有事件A和事件B,可以同时发生但不是完全同时发生,如以下韦恩图所示: 其中,A∩B表示A和B的并集,即A和B同时发生的概率. 如此,我们很容易得出,在事件B发生的情况下,事件A发生的概率为: 这个 ...
- Day1 Markdown学习!
Markdown学习 标题 一级标题:# (空格)+内容 二级标题:##(空格)+内容 同理可支持到六级标题 字体 Hello,World! 两边两个** 加粗 Hello,World! 两边一个* ...
- war项目依赖war项目
还没有看,立个flag:https://my.oschina.net/u/588379/blog/1817077
- PO封装设计模式 -- Web页面端测试
一.已登录页面 -->新建PO封装的包 -- 以下源码适用于python3以上的版本 代码优化新增 Image -->对操作步骤进行截图 二.basepage 包基础类的封装如下: fro ...
- Winform框架中窗体基类的用户身份信息的缓存和提取
在Winform开发中,有时候为了方便,需要把窗体的一些常规性的数据和操作函数进行封装,通过自定义基类窗体的方式,可以实现这些封装管理,让我们的框架统一化.简单化的处理一些常规性的操作,如这里介绍的用 ...