QBXT模拟赛T3

最新推荐文章于 2021-01-13 15:35:00 发布
最新推荐文章于 2021-01-13 15:35:00 发布
阅读量72 收藏
点赞数
原文链接：http://www.cnblogs.com/bullshit/p/9884569.html
版权

set+树上操作

这道题主要是用set维护每个点的01序列，然后启发式合并

但是我tm（拒绝口吐莲花）调了好几个小时！！！！！！！！！！！！！！！！！！！！！！！！！

主要是set的使用问题：
1、set一定要插入一个极小值，不然可能会越界RE
2、set使用的迭代器还是很好用的，并且迭代器只能++或者--，最好用一下中间变量存一下
3、set的迭代器注意不要用混了！！！
4、set的insert插入一个结构体时注意元素位置不要打错了

代码：

#include<iostream>
#include<cstdio>
#include<queue>
#include<algorithm>
#include<cmath>
#include<set>
#include<cstring>
#define max(a,b) a>b?a:b
#define min(a,b) a<b?a:b
#define maxn 700006
#define mod 1000000007
#define pb push_back
#define all (zmd)
#define pot set<zmd>::iterator
#define rep(i,a,b) for (int i=a;i<=b;++i)
#define erep(i,a) for (int i=head[a];i!=-1;i=e[i].next)
using namespace std;
#define int long long 
struct zmd
{
    int id,l,r;
    bool operator < (const zmd &a) const 
    {
        return r<a.r;
    }
};
set<zmd>st[maxn];
struct hzw
{
    int to,next;    
}e[maxn];
inline int dc(int k)
{
    return (1+k)*k/2;
}
int head[maxn],cur,n,real[maxn],ans,size[maxn],cost[maxn];
inline void add(int a,int b)
{
    e[cur].to=b;
    e[cur].next=head[a];
    head[a]=cur++;
}
inline void FFT(int a,int b)
{   
    for (pot it = st[a].begin();it != st[a].end();++it)
    {
        if (it == st[a].begin())continue;
        if (!it->id) continue;
        int ll = it->l,rr = it->r;//dang qian yao cha ru de qu jian
        
        pot it2 = st[b].lower_bound(all{1,ll,rr}); //fa xian yi ge di yi ge da yu de weizhi 
        int nowl = it2->l,nowr = it2->r; //now
        pot up = ++it2;it2--;
        pot down = --it2;it2++;
        cost[b]-=dc(nowr-nowl+1); 
        if (rr != nowr && ll != nowl) 
        {
            size[b] += 2;
            st[b].erase(it2);
            st[b].insert(all{1,ll,rr});
            cost[b]+=dc(rr-ll+1);
            st[b].insert(all{0,nowl,ll-1});
            cost[b]+=dc(ll-nowl);
            st[b].insert(all{0,rr+1,nowr});
            cost[b]+=dc(nowr-rr);
        }
        else if (ll == nowl && rr == nowr) 
        {
            
            if ( down != st[b].begin()) 
            {
                ll = down->l;
                cost[b]-=dc(down->r-down->l+1);
                st[b].erase(down);
            }
            if (up != st[b].end())
            {
                rr = up->r;
                cost[b]-=dc(up->r-up->l+1);
                st[b].erase(up);
            }
            st[b].erase(it2);
            st[b].insert(all{1,ll,rr});
            cost[b]+=dc(rr-ll+1);
        }
        else 
        {
            if (down != st[b].begin() && nowl == ll) 
            {
                ll = (down)->l;
                cost[b]-=dc(down->r-down->l+1);
                st[b].erase(down);
                st[b].erase(it2);
                cost[b]+=dc(rr-ll+1);
                st[b].insert(all{1,ll,rr});
                cost[b]+=dc(nowr-rr);
                st[b].insert(all{0,rr+1,nowr});
            }
            else if (up != st[b].end() && nowr == rr)
            {
                rr = (up)->r;
                cost[b]-=dc(up->r-up->l+1);
                st[b].erase(up);
                st[b].erase(it2);
                cost[b]+=dc(rr-ll+1);
                st[b].insert(all{1,ll,rr});
                cost[b]+=dc(ll-nowl);
                st[b].insert(all{0,nowl,ll-1});
            }
            else 
            {                   
                st[b].erase(it2);
                if (nowl == ll) 
                {
                    st[b].insert(all{1,ll,rr});
                    st[b].insert(all{0,rr+1,nowr});
                    cost[b]+=dc(rr-ll+1)+dc(nowr-rr);
                }
                else
                {
                    st[b].insert(all{1,ll,rr});
                    st[b].insert(all{0,nowl,ll-1});
                    cost[b]+=dc(rr-ll+1)+dc(ll-nowl);
                }
                
            }
        }
    }
}
inline void dfs(int s,int fa)
{
    int mx = -1,id = s;
    erep(i,s)
    {
        if (e[i].to == fa) continue;
        dfs(e[i].to,s);
        if (size[real[e[i].to]]>mx)
        {
            mx = size[real[e[i].to]]; 
            id = real[e[i].to];
        }
    }
    real[s]=id;
    erep(i,s)
    {
        if (e[i].to==fa) continue;
        if (real[e[i].to] == id) continue;
        FFT(real[e[i].to],id);
        st[real[e[i].to]].clear();
    }
    st[n+1].clear(),st[n+2].clear();
    st[n+1].insert(all{0,1,n}),st[n+1].insert(all{-1,-1,-1});
    st[n+2].insert(all{1,s,s}),st[n+2].insert(all{-1,-1,-1});   
    FFT(n+2,n+1);
    FFT(n+1,id);
    if (s!=1) ans += cost[id];      

}
#undef int
int main()
{
    #define int long long 
    freopen("bye.in","r",stdin);
    freopen("bye.out","w",stdout);
    memset(head,-1,sizeof(head));
    cin>>n;
    rep(i,1,n+2) st[i].insert(all{-1,-1,-1}),st[i].insert(all{0,1,n}),cost[i] = dc(n),size[i] = 1;
    rep(i,1,n-1)
    {
        int a,b;
        scanf("%lld%lld",&a,&b);
        add(a,b);
        add(b,a);
    }
    dfs(1,1);
    cout<<ans;
}

代码2：

#include<iostream>
#include<cstdio>
#include<queue>
#include<algorithm>
#include<cmath>
#include<set>
#include<cstring>
#define max(a,b) a>b?a:b
#define min(a,b) a<b?a:b
#define maxn 700006
#define mod 1000000007
#define pb push_back
#define all (zmd)
#define pot set<zmd>::iterator
#define rep(i,a,b) for (int i=a;i<=b;++i)
#define erep(i,a) for (int i=head[a];i!=-1;i=e[i].next)
using namespace std;
#define int long long
struct zmd {int id,l,r;bool operator < (const zmd &a) const {return r<a.r;}};
set<zmd>st[maxn];
struct hzw {int to,next;} e[maxn];
inline int dc(int k) {return (1+k)*k/2;}
int head[maxn],cur,n,real[maxn],ans,size[maxn],cost[maxn];
inline void add(int a,int b) {e[cur].to=b;e[cur].next=head[a];head[a]=cur++;}
inline void FFT(int a,int b) {
    for (pot it = st[a].begin(); it != st[a].end(); ++it) {
        if (it == st[a].begin())continue;
        if (!it->id) continue;
        int ll = it->l,rr = it->r;//dang qian yao cha ru de qu jian
        pot it2 = st[b].lower_bound(all {1,ll,rr}); //fa xian yi ge di yi ge da yu de weizhi
        int nowl = it2->l,nowr = it2->r; //now
        pot up = ++it2;it2--;pot down = --it2;it2++;cost[b]-=dc(nowr-nowl+1);
        if (rr != nowr && ll != nowl) {
            size[b] += 2;st[b].erase(it2);st[b].insert(all {1,ll,rr});cost[b]+=dc(rr-ll+1);
            st[b].insert(all {0,nowl,ll-1});cost[b]+=dc(ll-nowl);st[b].insert(all {0,rr+1,nowr});cost[b]+=dc(nowr-rr);
        } else if (ll == nowl && rr == nowr) {
            if ( down != st[b].begin()) {ll = down->l;cost[b]-=dc(down->r-down->l+1);st[b].erase(down);}
            if (up != st[b].end()) {rr = up->r;cost[b]-=dc(up->r-up->l+1);st[b].erase(up);}
            st[b].erase(it2);st[b].insert(all {1,ll,rr});cost[b]+=dc(rr-ll+1);
        } else {
            if (down != st[b].begin() && nowl == ll) {
                ll = (down)->l;cost[b]-=dc(down->r-down->l+1);st[b].erase(down);st[b].erase(it2);
                cost[b]+=dc(rr-ll+1);st[b].insert(all {1,ll,rr});cost[b]+=dc(nowr-rr);st[b].insert(all {0,rr+1,nowr});
            } else if (up != st[b].end() && nowr == rr) {
                rr = (up)->r;cost[b]-=dc(up->r-up->l+1);st[b].erase(up);st[b].erase(it2);cost[b]+=dc(rr-ll+1);
                st[b].insert(all {1,ll,rr});cost[b]+=dc(ll-nowl);st[b].insert(all {0,nowl,ll-1});
            } else {
                st[b].erase(it2);
                if (nowl == ll) {st[b].insert(all {1,ll,rr});st[b].insert(all {0,rr+1,nowr});cost[b]+=dc(rr-ll+1)+dc(nowr-rr);} 
                else {st[b].insert(all {1,ll,rr});st[b].insert(all {0,nowl,ll-1});cost[b]+=dc(rr-ll+1)+dc(ll-nowl);}
            }
        }
    }
}
inline void dfs(int s,int fa) {
    int mx = -1,id = s;
    erep(i,s) {if (e[i].to == fa) continue;dfs(e[i].to,s);if (size[real[e[i].to]]>mx) {mx = size[real[e[i].to]];id = real[e[i].to];}}
    real[s]=id;
    erep(i,s) {if (e[i].to==fa) continue;if (real[e[i].to] == id) continue;FFT(real[e[i].to],id);st[real[e[i].to]].clear();}
    st[n+1].clear(),st[n+2].clear();st[n+1].insert(all {0,1,n}),st[n+1].insert(all {-1,-1,-1});st[n+2].insert(all {1,s,s}),st[n+2].insert(all {-1,-1,-1});
    FFT(n+2,n+1);FFT(n+1,id);if (s!=1) ans += cost[id];
}
#undef int
int main() {
#define int long long
    freopen("bye.in","r",stdin);freopen("bye.out","w",stdout);
    memset(head,-1,sizeof(head));cin>>n;
    rep(i,1,n+2) st[i].insert(all {-1,-1,-1}),st[i].insert(all {0,1,n}),cost[i] = dc(n),size[i] = 1;
    rep(i,1,n-1) {int a,b;scanf("%lld%lld",&a,&b);add(a,b);add(b,a);}
    dfs(1,1);cout<<ans;
}

转载于:https://www.cnblogs.com/bullshit/p/9884569.html

weixin_30788239
关注
【模拟赛】qbxt金秋冲刺训练营杯NOIP模拟赛第一场总结
ars4me
09-20 377
前言 虽然是赛后自己测 但也是完完整整的打了一遍模拟赛 断断续续牵扯了两天半的精力
2017.10.1 QBXT 模拟赛
weixin_30838921的博客
10-10 74
题目链接 T1 枚举右端点,前缀和优化。对于当前点x，答案为 sum[x][r]-sum[x][l-1]-(sum[z][r]-sum[z][l-1]) 整理为 sum[x][r]-sum[z][r]-(sum[x][l-1]-sum[z][l-1]) 我们已知x和sum[x][r],对于z我们枚举，对于sum[x][l-1]-sum[z][l-1]我们需要一...
2017.10.6 QBXT 模拟赛
weixin_30834019的博客
10-10 62
题目链接 T1 Sort 一下与原数组比较 ，若有两个数或者没有数发生位置交换 ，则输出YES ，否则输出NO #include <algorithm> #include <cctype> #include <cstdio> #define N 1005000 int n, cnt1, cnt2, sum; struct node...
2017.10.4 QBXT 模拟赛
weixin_30481087的博客
10-10 81
题目链接 T1 维护一个单调栈 #include <iostream> #include <cstdio> #define N 500000 #define rep(a,b,c) for(int a=b;a<=c;++a) #define Rep(a,b,c) for(int a=b;a>=c;--a) using namespac...
2017.10.2 QBXT 模拟赛
weixin_30376083的博客
10-10 60
题目链接 T1 我们所要求得是(a*b)|x 也就是 使(a*b)的倍数小于x的个数之和 1<=x<=n 我们可以 找一个c使得 (a*b*c)<=x 由于我们所求的是一个三元有序对 即 (1,2,3) 与 (1,3,2) 是两种不同的方案数 所以 我们可以强制规定 a<b<c 最坏的情况是(a*a*a)==x 所以我们就可以确定a的枚举范围...
2017.10.7 QBXT 模拟赛
weixin_30865427的博客
10-10 67
题目链接 T1 容斥原理，根据奇偶性进行加减 #include<iostream> #include<cstdio> using namespace std; typedef long long ll; typedef long double ld; #define rep(i,a,n) for(int i=a;i<=n;i++) ld...
2017.10.3 QBXT 模拟赛
weixin_30670965的博客
10-10 58
题目链接 T1 模拟 #include <cstring> #include <cstdio> #define N 105000 int L,R; char s[N]; int main() { freopen("bracket.in","r",stdin); freopen("bracket.out","w",stdout);...
【模拟赛】qbxt金秋冲刺训练营杯NOIP模拟赛第二场D1
ars4me
09-23 226
这一次是正儿八经参加啦
Qbxt 模拟题 day3(am) T3 选数字 (select)(贪心)
weixin_30668887的博客
11-07 102
选数字 (select Time Limit:3000ms Memory Limit:64MB 题目描述 LYK 找到了一个 n*m 的矩阵，这个矩阵上都填有一些数字，对于第 i 行第 j 列的位置上 的数为 ai,j。 由于它 AK 了 noip2016 的初赛，最近显得非常无聊，便想到了一个方法自娱自乐一番。 ...
NOIP 2018 提高组模拟（QBXT版）
08-15
【标题】"NOIP 2018 提高组模拟（QBXT版）" 指的是全国青少年信息学奥林匹克竞赛（National Olympiad in Informatics, NOIP）2018年提高组的一次模拟赛。这个模拟赛是针对QBXT（可能是一个解题平台或训练团队的代号...
NOIP模拟赛 T3:与众不同(线段树+滑动窗口)
Continue
11-16 2212
NOIP模拟赛 T3:与众不同问题描述:A是某公司的CEO,每个月都会有员工把公司的盈利数据送给A,A是个与众不同的怪人,A不注重盈利还是亏本,而是喜欢研究”完美序列”:连续的互不相同的序列.A想知道区间[L,R]之间最长的完美序列.输入格式:第一行两个整数N,M(1<=N,M<=200000)，N表示连续N个月,编号为0到N-1，M表示询问的次数。 第二行N个整数(绝对值不超过10^6),第i个
2021年 第12届 蓝桥杯 第2次模拟赛真题详解及小结【Java版】
upward
01-13 8410
模拟赛网页截图 一、试题A——答案：168 二、试题B——答案：14 三、试题C——答案：2002 四、试题D——答案：65 解法一：找规律（累加求和） 五、试题E——答案：162 解法一：双重for循环遍历所有组合 六、试题F 解法一：整数相乘 七、试题G 解法一：整数运算 八、试题H 解法一：【求巨佬代码】 九、试题I 解法一：【求巨佬代码】 十、试题J 小结
Java实现 第十一届 蓝桥杯 （本科组）省内模拟赛
热门推荐
南 墙
04-13 1万+
有错误的或者有问题的欢迎评论 计算机存储中有多少字节 合法括号序列 无向连通图最少包含多少条边 字母重新排列 凯撒密码加密 反倍数 正整数的摆动序列 螺旋矩阵 小明植树 户户通电 计算机存储中有多少字节 题目 问题描述 在计算机存储中，12.5MB是多少字节？ 答案提交 这是一道结果填空的题，你只需要算出结果后提交即可。本题的结果为一个整数，在提交答案时只填写这个整数，填写多余的内容将无法得分。 ...
C# winform使用纯opencvsharp部署yolov11-onnx图像分类模型源码
10-06
【测试环境】 vs2019 net framework4.7.2 opencvsharp4.8.0 【视频演示】 https://www.bilibili.com/video/BV1EB1iYXEoi/ 【博文地址】 https://blog.csdn.net/FL1623863129/article/details/142728931
现代密码学课程实验.zip
10-06
项目学习分享。【项目资源】：包含前端、后端、移动开发、操作系统、人工智能、物联网、信息化管理、数据库、硬件开发、大数据、课程资源、音视频、网站开发等各种技术项目的源码。包括STM32、ESP8266、PHP、QT、Linux、iOS、C++、Java、python、web、C#、EDA、proteus、RTOS等项目的源码。 【技术】 Java、Python、Node.js、Spring Boot、Django、Express、MySQL、PostgreSQL、MongoDB、React、Angular、Vue、Bootstrap、Material-UI、Redis、Docker、Kubernetes
2022人工智能技术创新大赛-赛道1-电商关键属性匹配.zip
最新发布
10-06
项目学习分享。【项目资源】：包含前端、后端、移动开发、操作系统、人工智能、物联网、信息化管理、数据库、硬件开发、大数据、课程资源、音视频、网站开发等各种技术项目的源码。包括STM32、ESP8266、PHP、QT、Linux、iOS、C++、Java、python、web、C#、EDA、proteus、RTOS等项目的源码。 【技术】 Java、Python、Node.js、Spring Boot、Django、Express、MySQL、PostgreSQL、MongoDB、React、Angular、Vue、Bootstrap、Material-UI、Redis、Docker、Kubernetes
基于ssm大学图书管理系统.zip
10-06
基于ssm大学图书管理系统 食物管理.. 运动管理 用户管理 Android端 手机登陆 手机注册… 手机功能选择 个人健康信息 食物计算器、 生活记录本
水化放热COMSOL计算仿真 一区 Cement & Concrete Composites 半绝热和蒸汽养护条件下硅酸盐水泥砂
10-06
水化放热COMSOL计算仿真 一区 Cement & Concrete Composites 半绝热和蒸汽养护条件下硅酸盐水泥砂浆水化过程传热传质的数值模拟 水化放热 其他一样的均为直接搬运无源文件 COMSOL6.1 6.2版本
C语言编程：模拟乒乓球比赛
该系统利用8个LED发光管模拟乒乓球在两个球台间的移动，并通过两个按钮模拟发球和接球动作。比赛规则遵循一般的乒乓球比赛，包括规定时间内未完成击球动作的一方会失去一分。 1. 基本要求： - 球每0.5秒滚动过1个...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值