传送门

文章目录

  • 题意:
  • 思路:

题意:

你现在有ppp种电台,有nnn对关系(x,y)(x,y)(x,y)代表xxx电台或yyy电台中至少有一个,mmm对关系(x,y)(x,y)(x,y)代表xxx电台或yyy电台中最多有一个,每个电台有两个参数li,ril_i,r_ili​,ri​,你需要在[1,M][1,M][1,M]中选择一个主频fff,如果f∈[li,ri]f\in [l_i,r_i]f∈[li​,ri​],那么第iii个电台你可以选择是否启用,否则无法启用。

请给出fff并构造电台启用方案,无解输出−1-1−1。

2≤n,p,m,M≤4e52\le n,p,m,M\le 4e52≤n,p,m,M≤4e5

思路:

对于n,mn,mn,m这几个条件,显然是2−sat2-sat2−sat板子了,以下设uuu代表选uuu这个点,u′u'u′代表不选uuu这个点。

对于nnn对关系,我们建x′−>y,y′−>xx'->y,y'->xx′−>y,y′−>x。

对于mmm对关系,我们建x−>y′,y−>x′x->y',y->x'x−>y′,y−>x′。

热身完毕,现在开始正菜。

考虑如何制定fff呢?我们运用前缀和的思想,尝试将[l,r][l,r][l,r]拆成[1,l−1],[1,r][1,l-1],[1,r][1,l−1],[1,r]。

我们把f∈[1,li−1],f∈[1,ri]f\in [1,l_i-1],f\in [1,r_i]f∈[1,li​−1],f∈[1,ri​]的真值表写出来。

f∈[1,li−1]f∈[1,ri]i能否启用00001110不存在110\begin{array}{ccc} f\in [1,l_i-1] & f\in [1,r_i] & i能否启用 \\ \hline \ 0&0&0\\ \ 0&1&1\\ \ 1&0&不存在\\ \ 1&1&0\\\\\end{array}f∈[1,li​−1] 0 0 1 1​f∈[1,ri​]0101​i能否启用01不存在0​​

我们可以发现总结一下就是以下三个限制:

(1)(1)(1)若f∈[1,li−1]f\in [1,l_i-1]f∈[1,li​−1],则iii不能启用。

(2)(2)(2)若f∈[1,ri]f\in [1,r_i]f∈[1,ri​]不满足,则iii不能启用。

(3)(3)(3) 若iii启用,则f∈[li,ri]f\in [l_i,r_i]f∈[li​,ri​]。

所以我们拿出来[n+1,n+M+1][n+1,n+M+1][n+1,n+M+1]来代表[0,M][0,M][0,M],我们用这个来继续限制。

我们设满足n+i+1n+i+1n+i+1的时候就是f≤if\le if≤i,否则f>if>if>i。

首先就是n+i+1−>n+i+2,(n+i+2)′−>(n+i+1)′n+i+1->n+i+2,(n+i+2)'->(n+i+1)'n+i+1−>n+i+2,(n+i+2)′−>(n+i+1)′,这个比较显然,因为f≤if\le if≤i,那么f≤i+1f\le i+1f≤i+1,如果f>i+1f>i+1f>i+1,那么f>if>if>i。

让后就是限制一下上面的条件

(1)(1)(1) n+li−>i′,i−>(n+li)′n+l_i->i',i->(n+l_i)'n+li​−>i′,i−>(n+li​)′。

(2)(2)(2)(n+ri+1)′−>i′,i−>n+ri+1(n+r_i+1)'->i',i->n+r_i+1(n+ri​+1)′−>i′,i−>n+ri​+1。

还有就是fff不能为000,这个用n+1−>(n+1)′n+1->(n+1)'n+1−>(n+1)′限制一下即可,算是2−sat2-sat2−sat的小技巧了。

跑完之后,n+f+1n+f+1n+f+1成立而n+fn+fn+f不成立的时候的分界点就是fff。

// Problem: F. Radio Stations
// Contest: Codeforces - Codeforces Round #585 (Div. 2)
// URL: https://codeforces.com/contest/1215/problem/F
// Memory Limit: 256 MB
// Time Limit: 7000 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=2000010,M=N*10,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int p,n,q,m;
int e[M],ne[M],h[N],idx;
int dfn[N],low[N],id[N],tot,cnt;
int stk[N],top;
bool in[N];void add(int a,int b)
{e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}void tarjan(int u)
{dfn[u]=low[u]=++tot;stk[++top]=u; in[u]=true;for(int i=h[u];~i;i=ne[i]){int ver=e[i];if(!dfn[ver]){tarjan(ver);low[u]=min(low[u],low[ver]);}else if(in[ver]) low[u]=min(low[u],dfn[ver]);}if(dfn[u]==low[u]){int y; ++cnt;do{y=stk[top--];in[y]=false; id[y]=cnt;}while(y!=u);}
}int yes(int x) { return x<<1; }int no(int x) { return x<<1|1; }void link(int x,int y) {add(x,y); add(y^1,x^1);
}bool check()
{for(int i=0;i<=n+m;i++) if(id[yes(i)]==id[no(i)]) return false;return true;
}int main()
{//  ios::sync_with_stdio(false);
//  cin.tie(0);scanf("%d%d%d%d",&p,&n,&m,&q);idx=0;memset(h,-1,sizeof(h));for(int i=1;i<=p;i++) {int a,b; scanf("%d%d",&a,&b);a--; b--;// add(2*a,2*b+1); // add(2*b,2*a+1);link(no(a),yes(b));}for(int i=0;i<m;i++) {// add(n*2+i*2+1,n*2+(i+1)*2+1),add(n*2+(i+1)*2,n*2+i*2);link(yes(n+i),yes(n+i+1));}// link(yes(n),no(n));// add(n*2,n*2+1);// add(n*2+1,n*2);for(int i=0;i<n;i++) {int l,r; scanf("%d%d",&l,&r);l--; link(yes(i),no(n+l));link(yes(i),yes(n+r));// add(n*2+l*2+1,i*2); add(i*2+1,n*2+l*2);// add(n*2+r*2,i*2); add(i*2+1,n*2+r*2+1);}for(int i=1;i<=q;i++) {int a,b; scanf("%d%d",&a,&b);a--; b--;link(yes(a),no(b));// add(a*2+1,b*2);// add(b*2+1,a*2);}//n*2+(m+1)*2for(int i=0;i<2*n+(m+1)*2;i++) if(!dfn[i]) tarjan(i);if(check()) {vector<int>ans;for(int i=0;i<n;i++) if(id[yes(i)]<id[no(i)]) ans.pb(i+1);for(int i=1;i<=m;i++) if(id[yes(n+i)]<id[no(n+i)]) {printf("%d %d\n",ans.size(),i);break;}for(auto x:ans) printf("%d ",x); puts("");} else puts("-1");return 0;
}

Codeforces Round #585 (Div. 2) F. Radio Stations 2-sat + 神仙建模相关推荐

  1. Codeforces Round #585 (Div. 2) E. Marbles 状压dp + 逆序对

    传送门 文章目录 题意: 思路: 题意: 思路: 考虑数列最终的状态一定是相同颜色在一起,所以我们发现他的颜色是有顺序的!显然可以用状压dpdpdp来枚举颜色的顺序,但是又有问题了,你怎么确定当前这个 ...

  2. Codeforces Round #699 (Div. 2) F - AB Tree(贪心、树上DP)超级清晰,良心题解,看不懂来打我 ~

    整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 Codeforces Round #699 (Div. 2) F - AB Tree Problem ...

  3. Codeforces Round #644 (Div. 3) F.Spy-string

    Codeforces Round #644 (Div. 3) F.Spy-string 题目链接 You are given n strings a1,a2,-,an: all of them hav ...

  4. Codeforces Round #849 (Div. 4) F. Range Update Point Query

    Codeforces Round #849 (Div. 4) F. Range Update Point Query 题目大意: 给一串数字,有两个操作: 操作1:将 l − r l-r l−r 的数 ...

  5. Codeforces Round #538 (Div. 2) F. Please, another Queries on Array? 线段树 + 欧拉函数

    传送门 文章目录 题意: 思路: 题意: 给你一个序列aaa,你需要实现两种操作: (1)(1)(1) 将[l,r][l,r][l,r]的aia_iai​都乘rrr. (2)(2)(2) 求ϕ(∏i= ...

  6. Codeforces Round #742 (Div. 2) F. One-Four Overload 构造 + 二分图染色

    传送门 文章目录 题意: 思路: 题意: 给你一个n∗mn*mn∗m的矩形,包含...和XXX,你有两种颜色,你需要给...染色使得每个XXX上下左右相邻的...其两种颜色个数相同,输出一种合法方案. ...

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

    传送门 文章目录 题意: 思路: 题意: 思路: 考虑最终的序列是什么鸭子的,首先序列肯定单调不降,也就是a1≤a2≤a3≤...≤ana_1\le a_2\le a_3\le ...\le a_na ...

  8. Codeforces Round #675 (Div. 2) F. Boring Queries 区间lcm + 主席树

    传送门 文章目录 题意: 思路: 题意: 给你一个长度为nnn的序列aaa,qqq个询问,每次询问[l,r][l,r][l,r]内的lcmlcmlcm是多少,对1e9+71e9+71e9+7取模. n ...

  9. Codeforces Round #592 (Div. 2) F. Chips 构造 + 细节

    传送门 文章目录 题意: 思路: 题意: 思路: 恶心的构造题,思路很简单但是代码细节很多,搞了半天. 根据题目的性质不难发现,如果有两个相同颜色的球相邻,那么他们的颜色永远不会改变. 根据这个性质, ...

最新文章

  1. 大致了解写java聊天器所需要的技术
  2. jsp中未登录用户也可以浏览页面的功能实现代码
  3. Spring-学习笔记03【Spring的IOC和DI】
  4. C语言 十进制和八进制相互转换 - C语言零基础入门教程
  5. Go 的 Contex 是线程安全的吗?
  6. UI设计师用的最多的素材网站都收纳在一个导航
  7. Matlab取整函数
  8. 中国五千年统一与分裂
  9. matlab动刚度仿真,基于ANSYS的转子支承动刚度计算分析.doc
  10. 可视化设计之迷失扁平化风潮
  11. Matlab超声影像智能硬件
  12. Python(1-10章)知识点
  13. photoswipe 移动端PC端兼容图片相册预览
  14. python列表(list)的遍历
  15. PSIM仿真+C语言(c block)实现PID控制和最小拍控制
  16. c++入门必学算法 质数筛
  17. GOM引擎登录器配置教程
  18. unity保存游戏对象到mysql_非常好用的游戏数据保存类:ScriptableObject
  19. 为什么这么多人转行产品经理?产品经理发展前景如何?
  20. 习题4.10 设有一长为3000m的绳子,每天减去一半,问需几天时间,绳子的长度会短于5m。

热门文章

  1. python123可以复制粘贴吗_使用win32com和python复制和粘贴隔离
  2. 第7章 C控制语句:分支和跳转
  3. 一元二次方程求根公式的花样变换,你看懂了吗?
  4. 千万别让爸妈帮你P图......
  5. 颜宁问4对科研夫妻:男科学家怎样平衡事业家庭?
  6. 当女朋友学会「监视」男朋友......
  7. 假如,绿茶婊的目标变成女生......
  8. 2018年最后一个月最值得关注的13个优质公号
  9. android 拍照换头像,Android调用相机拍照,裁剪及更换头像功能的实现
  10. qt4.7 mysql_详解Qt 4.7编译和访问Mysql驱动