Berland Army CodeForces - 883B (贪心,拓扑排序)
大意: n个点, 点$i$的等级为$r_i$, 只给出部分点的$r$值, $r_i$的范围为[1,k], 且[1,k]都至少有一个. 给定m条有向边, (x,y)表示$r[x]>r[y]$, 求一个合法的$r$序列.
我们记条件(1): [1,k]都至少有一个. 条件(2): 所有有向边(x,y), 满足$r[x]>r[y]$.
正反各拓扑一次, 这样可以求出满足条件(2)时每个点的可用等级范围[L,R], 我们考虑在不违背条件(2)的情况下, 使(1)成立.
所有点按$L$排序后, 将$1,...,k$从小到大分配给各点, 假设遍历到$i$, $[1,i-1]$都已经分配完毕, 用一个集合$S$维护所有未分配的且$L\le i$的点. 对于$S$中所有$R\le i$的点, 显然全分配为$i$会更优, 这是因为如果分配为$j<i$的话, 可能会与之前的点违背条件(2). 而对于$S$中所有$R> i$的点, 我们尽量先不分配一定不会影响条件(2), 这是由于之前拓扑排序保证的. 需要注意如果不存在$R \le i$的点的话, 那么必须从$R>i$中选择一个, 显然选择一个$R$小的会使条件(1)更优
#include <iostream>
#include <sstream>
#include <algorithm>
#include <cstdio>
#include <math.h>
#include <set>
#include <map>
#include <queue>
#include <string>
#include <string.h>
#include <bitset>
#define REP(i,a,n) for(int i=a;i<=n;++i)
#define PER(i,a,n) for(int i=n;i>=a;--i)
#define hr putchar(10)
#define pb push_back
#define lc (o<<1)
#define rc (lc|1)
#define mid ((l+r)>>1)
#define ls lc,l,mid
#define rs rc,mid+1,r
#define x first
#define y second
#define io std::ios::sync_with_stdio(false)
#define endl '\n'
#define DB(a) ({REP(__i,1,n) cout<<a[__i]<<' ';hr;})
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
const int P = 1e9+7, INF = 0x3f3f3f3f;
ll gcd(ll a,ll b) {return b?gcd(b,a%b):a;}
ll qpow(ll a,ll n) {ll r=1%P;for (a%=P;n;a=a*a%P,n>>=1)if(n&1)r=r*a%P;return r;}
ll inv(ll x){return x<=1?1:inv(P%x)*(P-P/x)%P;}
inline int rd() {int x=0;char p=getchar();while(p<'0'||p>'9')p=getchar();while(p>='0'&&p<='9')x=x*10+p-'0',p=getchar();return x;}
//head #ifdef ONLINE_JUDGE
const int N = 1e6+10;
#else
const int N = 111;
#endif int n, m, k, a[N], deg1[N], deg2[N];
int L[N], R[N], vis[N];
int x[N], y[N], b[N];
vector<int> g1[N], g2[N]; int main() {
scanf("%d%d%d", &n, &m, &k);
REP(i,1,n) {
scanf("%d", a+i);
if (a[i]) L[i]=R[i]=a[i];
else L[i]=1,R[i]=k;
}
REP(i,1,m) {
scanf("%d%d",x+i,y+i);
g2[x[i]].pb(y[i]),++deg2[y[i]];
g1[y[i]].pb(x[i]),++deg1[x[i]];
}
queue<int> q;
REP(i,1,n) if (!deg1[i]) q.push(i);
int cnt = 0;
while (q.size()) {
int u = q.front(); q.pop();
if (L[u]>R[u]) return puts("-1"),0;
++cnt;
for (auto &v:g1[u]) {
L[v] = max(L[v], L[u]+1);
if (!--deg1[v]) q.push(v);
}
}
if (cnt!=n) return puts("-1"),0;
REP(i,1,n) if (!deg2[i]) q.push(i);
while (q.size()) {
int u = q.front(); q.pop();
if (L[u]>R[u]) return puts("-1"),0;
for (auto &v:g2[u]) {
R[v] = min(R[v], R[u]-1);
if (!--deg2[v]) q.push(v);
}
}
REP(i,1,n) b[i] = i;
sort(b+1,b+1+n,[](int a,int b){return L[a]<L[b];});
int now = 1;
set<pii> s;
REP(i,1,k) {
while (now<=n&&L[b[now]]<=i) {
s.insert(pii(R[b[now]],b[now]));
++now;
}
int ok = 0;
while (s.size()) {
int id = s.begin()->y;
if (ok&&R[id]>i) break;
ok = 1;
s.erase(s.begin());
a[id] = i;
}
if (!ok) return puts("-1"),0;
}
REP(i,1,n) printf("%d ", a[i]);hr;
}
Berland Army CodeForces - 883B (贪心,拓扑排序)的更多相关文章
- [POI2004] SZP (贪心+拓扑排序)
[问题描述] Byteotian 中央情报局(BIA) 雇佣了许多特工. 他们每个人的工作就是监视 另一名特工. Byteasar 国王需要进行一次秘密行动,所以他要挑选尽量多的信得过的特工. 但 是 ...
- hdu-5695 Gym Class(贪心+拓扑排序)
题目链接: Gym Class Time Limit: 6000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) ...
- Aizu 2456 Usoperanto 贪心 拓扑排序
Usoperanto Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://www.bnuoj.com/v3/contest_show.php?cid= ...
- 2016"百度之星" - 初赛(Astar Round2A)1006 Gym Class(HDU5695)——贪心+拓扑排序
分析:首先,利用贪心可知,如果要所有人的分数和最高,需要把序号大的优先放在前面.其次,对于a的前面不能为b,那么只能a在b前面了,那么就建立一条从a到b的边,并且b的入度加1.然后就是拓扑排序了.要分 ...
- Codeforces 510C (拓扑排序)
原题:http://codeforces.com/problemset/problem/510/C C. Fox And Names time limit per test:2 seconds mem ...
- E - E CodeForces - 1100E(拓扑排序 + 二分)
E - E CodeForces - 1100E 一个n个节点的有向图,节点标号从1到n,存在m条单向边.每条单向边有一个权值,代表翻转其方向所需的代价.求使图变成无环图,其中翻转的最大边权值最小的方 ...
- CodeForces 721C Journey(拓扑排序+DP)
<题目链接> 题目大意:一个DAG图有n个点,m条边,走过每条边都会花费一定的时间,问你在不超过T时间的条件下,从1到n点最多能够经过几个节点. 解题分析:对这个有向图,我们进行拓扑排序, ...
- poj 3687 Labeling Balls - 贪心 - 拓扑排序
Windy has N balls of distinct weights from 1 unit to N units. Now he tries to label them with 1 to N ...
- CodeForces - 919D Substring (拓扑排序+dp)
题意:将一个字符串上的n个字符视作点,给出m条有向边,求图中路径上最长出现的相同字母数. 分析:首先如果这张图中有环,则可以取无限大的字符数,在求拓扑排序的同时可以确定是否存在环. 之后在拓扑排序的结 ...
随机推荐
- AtCoder AGC017C Snuke and Spells
题目链接 https://atcoder.jp/contests/agc017/tasks/agc017_c 题解 很久前不会做看了题解,现在又看了一下,只想说,这种智商题真的杀我... 转化成如果现 ...
- H5 video全屏与取消全屏兼容
H5 video全屏与取消全屏各浏览器兼容, requestFullscreen()全屏方法,exitFullscreen()退出全屏方法.兼容各个浏览器与css3兼容一样加个前缀即可. // 全屏 ...
- [CSP-S模拟测试]:硬币(博弈论+DP+拓展域并查集)
题目传送门(内部题135) 输入格式 第一行包含一个整数$T$,表示数据组数. 对于每组数据,第一行两个整数$h,w$,表示棋盘大小. 接下来$h$行,每行一个长度为$w$的字符串,每个位置由为$o, ...
- Dynamic Web Module 2.3 change to 3.0 maven工程解决办法
pom.xml <plugins> <plugin> <groupId>org.apache.maven.plugins</groupId> <a ...
- ztree复选框
var setting = { check: { enable: true // chkboxType : { "Y" : "", "N" ...
- iOS13即将到来,iOS推送Device Token适配方案详解
关于提前适配iOS13 苹果推送DeviceToken的通知 随着苹果iOS13系统即将发布,个推提前推出DeviceToken适配方案,以确保新版本的兼容与APP推送服务的正常使用.iOS13的一个 ...
- POI 生成excel
POI生成原生Excel-工具类 https://www.jianshu.com/p/2dfe7fe7d02e JAVA poi 帮助类 https://www.cnblogs.com/Ca ...
- springboot的优点
2013年12月12日,spring发布了4.0版本.这个本来只是作为Java平台上的控制反转容器的库,经过将近10年的发展已经成为了一个巨无霸产品.不过其依靠良好的分层设计,每个功能模块都能保持较好 ...
- Fatal error: Allowed memory size of 134217728 bytes exhausted (tried to allocate 1099 bytes) in
解释是可用内存已耗尽,这关系到PHP的memory_limit的设置问题. 我在网上看到,有两种方法解决 1.修改php.ini memory_limit = 128 这种方法需要重启服务器,很显然, ...
- leetcode 分割回文串
这个方法有问题,这是计算所有子串组成的所有回文子串:而不是所有分割的回文子串: class Solution { public: vector<vector<string>> ...