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

传送门

文章目录

  • 题意:
  • 思路:

题意:

你现在有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-11

2≤n,p,m,M≤4e52\le n,p,m,M\le 4e52n,p,m,M4e5

思路:

对于n,mn,mn,m这几个条件,显然是2−sat2-sat2sat板子了,以下设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,l1],[1,r]

我们把f∈[1,li−1],f∈[1,ri]f\in [1,l_i-1],f\in [1,r_i]f[1,li1],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,li1] 0 0 1 1f[1,ri]0101i010

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

(1)(1)(1)f∈[1,li−1]f\in [1,l_i-1]f[1,li1],则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 ifi,否则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 ifi,那么f≤i+1f\le i+1fi+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-sat2sat的小技巧了。

跑完之后,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;
}

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

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

相关文章

C. Orac and LCM(数论lcm, gcd)

C. Orac and LCM 思路 题目非常简单&#xff0c;就是求gcd(lcm(i,j))foriinrange(n),forjinrange(n),i<jgcd(lcm_(i,\ j))\ for\ i\ in\ range(n),\ for\ j\ in\ range(n),\ i\ <\ jgcd(lcm(​i, j)) for i in range(n), for j in range(n), i < j 对于包含a1a_1a1…

.net core 实现基于 cron 表达式的任务调度

.net core 实现基于 cron 表达式的任务调度Intro上次我们实现了一个简单的基于 Timer 的定时任务&#xff0c;详细信息可以看这篇文章 。但是使用过程中慢慢发现这种方式可能并不太合适&#xff0c;有些任务可能只希望在某个时间段内执行&#xff0c;只使用 timer 就显得不是那…

P6378 [PA2010] Riddle 2-sat + 前缀和优化建图

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你nnn个点mmm调变的无向图被分成kkk个部分&#xff0c;每个部分包含若干点&#xff0c;请选择一些关键点&#xff0c;使得每个部分恰好有一个关键点&#xff0c;且每条边至少有一个是关键点。 1≤k,w≤n≤…

E:Tree Queries(假树链剖分写法)

博客园地址 E&#xff1a;Tree Queries 思路 当我写完A完这道题后&#xff0c;百度了一下&#xff0c;发现好像没有人是用类树链剖分来写的&#xff0c;都是LCALCALCA&#xff0c;于是我就来水一篇树链剖分题解了。 第一步&#xff1a;贪心取点 我们可以发现&#xff0c;要…

ASP.NET Core Identity自定义数据库结构和完全使用Dapper而非EntityFramework Core

前言原本本节内容是不存在的&#xff0c;出于有几个人问到了我&#xff1a;我想使用ASP.NET Core Identity&#xff0c;但是我又不想使用默认生成的数据库表&#xff0c;想自定义一套&#xff0c;我想要使用ASP. NE Core Identity又不想使用EntityFramework Core。真难伺候&…

Educational Codeforces Round 17 E. Radio stations cdq分治 + 树状数组

传送门 文章目录题意思路&#xff1a;题意 有nnn个电台&#xff0c;对于每个电台iii有三个参数xi,ri,fix_i,r_i,f_ixi​,ri​,fi​&#xff0c;分别指他们的坐标、作用半径、频率。如果两个电台频率差值在kkk以内&#xff0c;并且他们的作用范围都能覆盖到彼此&#xff0c;那么…

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

Educational Codeforces Round 89 (Rated for Div. 2) A. Shovels and Swords 思路 题意非常简单&#xff0c;就是得到最多的物品嘛&#xff0c;我们假定a,ba, ba,b中aaa是最小的一个&#xff0c;分两种情况。 如果2∗a<b2 * a < b2∗a<b&#xff0c;那么我们只需…

什么是微服务?为什么你要用微服务?

前言最近几年微服务很火&#xff0c;大家都在建设微服务&#xff0c;仿佛不谈点微服务相关的技术&#xff0c;都显得不是那么主流了。近几年见识到身边朋友的很多公司和团队都在尝试进行微服务的改变&#xff0c;但很多团队并没有实际微服务踩坑经验&#xff0c;很多团队甚至强…

P5367 【模板】康托展开

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 存个板子 // Problem: P5367 【模板】康托展开 // Contest: Luogu // URL: https://www.luogu.com.cn/problem/P5367 // Memory Limit: 64 MB // Time Limit: 1200 ms // // Powered by …

微软开源基于.NET Core的量子开发工具包 QDK

微软最近开源了量子开发工具包&#xff08;Quantum Development Kit&#xff0c;QDK&#xff09;&#xff0c;旨在使“量子计算和算法开发对开发人员来说更容易、更透明”。微软 QDK 包括 Q#编译器、量子库和量子模拟器。微软在 2017 年底的 Ignite 大会上发布了量子开发工具包…

UVA11525 Permutation 逆康托展开

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 逆康托展开板子 // Problem: UVA11525 Permutation // Contest: Luogu // URL: https://www.luogu.com.cn/problem/UVA11525 // Memory Limit: 0 MB // Time Limit: 3000 ms // // Power…

Azure 上使用 Windows Server Core 运行 ASP.NET Core 网站

点击上方蓝字关注“汪宇杰博客”导语微软智慧云 Azure 上虽然早就有 App Service 这种完全托管的 PaaS 服务可以让我们分分钟建网站。但是不自己配一下环境&#xff0c;就不能体现技术含量&#xff0c;容易被说微软的人都只会点鼠标。年轻的时候不敲命令&#xff0c;什么时候可…

Codeforces Round #285 (Div. 2) D. Misha and Permutations Summation 康托展开 + 线段树

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先肯定不能模n!n!n!&#xff0c;所以考虑先将a,ba,ba,b做一个逆康托展开&#xff0c;得到a′,b′a,ba′,b′数组&#xff0c;以及a′b′sumabsuma′b′sum数组&#xff0c;让后我们可以通…

「PowerBI」使用TabularEditor进行PowerBIDeskTop模型开发最佳实践

前面系列文章介绍的场景&#xff0c;设定的工具使用对象是Sqlserver和Azure 的SSAS数据模型开发&#xff0c;其实TabularEditor亦可以有限度地使用在PowerBIDeskTop的模型开发上&#xff0c;本文简单介绍下其最佳的使用场景。PowerBIDeskTop模型不同于Sqlserver的SSAS模型虽然大…

E:Modular Stability(组合数)

Modular Stability 思路 (((xmoda1)moda2)……modak−1)modak(((xmodp1)modp2)……modpk−1)modpk(((x \mod a_1) \mod a_2) …… \mod a_{k - 1}) \mod a_{k} (((x \mod p_1) \mod p_2) …… \mod p_{k - 1}) \mod p_{k}(((xmoda1​)moda2​)……modak−1​)modak​(((xmodp…

E:K-periodic Garland(DP)

思路 每个点我们有两种决策&#xff0c;其值为0或1&#xff1a; 如果点我们放置0的话&#xff0c;我们有其前一位数字是零&#xff0c;或者其前一位数字是一。 如果这个点我们放置1的话&#xff0c;我们有其前面是按照每k个数字都出现一次1的排列&#xff0c;也有可能其前面…

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

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 考虑最终的序列是什么鸭子的&#xff0c;首先序列肯定单调不降&#xff0c;也就是a1≤a2≤a3≤...≤ana_1\le a_2\le a_3\le ...\le a_na1​≤a2​≤a3​≤...≤an​&#xff0c;显然不可…

沃尔玛招聘.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) 对于有不止一条路径能到这个点的点编号为…