洛谷P3588 - [POI2015]Pustynia

Portal

Description

给定一个长度为\(n(n\leq10^5)\)的正整数序列\(\{a_n\}\),每个数都在\([1,10^9]\)范围内,告诉你其中\(s\)个数,并给出\(m(m\leq2\times10^5)\)条信息。每条信息包含三个数\(L,R,k(Σk\leq 3\times10^5)\)以及\(k\)个正整数\(\{x_k\}\),表示\(a_L..a_R\)中,任意一个\(x\)均比剩下的\(R-L+1-k\)个数大(严格大于,即没有等号)。请任意构造出一组满足条件的方案,或者判断无解。

Solution

拓扑排序+线段树优化建图。
定义点权为\(val\);存在一条边权为\(w\)的边\((u,v)\)表示\(val[v]\geq val[u]+w\)。
首先考虑朴素的做法。建立\(n\)个点,\(val[i]\)表示\(a_i\)的值。对于每一条信息,新建一个点\(p\),\(val[p]\)表示\(min\{x_k\}\);剩下的数分别向\(p\)连一条边权为\(1\)的边(\(min\{x_k\}\)大于剩下的数),\(p\)向\(x_1..x_k\)分别连一条边权为\(0\)的边(\(x_i\)大于等于\(min\{x_k\}\))。初始时入度为\(0\)的点若没有值则令其\(val=1\),然后进行拓扑排序,如果成环或与初值冲突则无解。这样共有\(O(nm)\)条边。
考虑到边权为\(1\)的边的起点相当于\(k+1\)个区间,我们可以用线段树来优化建图。举例:\(n=8,a_3=7,a_5=4,a_7=2\),\([1,4]\)中最大的是\(\{2,3\}\),\([4,8]\)中是\(\{6\}\),\([1,8]\)中是\(\{2\}\)。

其中虚线边的权值为\(0\),实线边的权值为\(1\)。对于蓝线以上的点(线段树上的点),其\(val\)表示区间中的最大值;对于蓝线以下的点(条件所代表的点),其\(val\)表示\(min\{x_k\}\)。意义还是很明确的:例如边\(([3,4],\{2\})\)的权值为\(1\),表示\(min\{a_2\}>max\{3,4\}\)。同样机型拓扑排序并判断无解即可。至于边数...线段树上有\(2n\)条边,\(m\)条信息总共划分了\(m+Σk\)个区间,每个区间对应\(O(logn)\)条边,\(\{x_k\}\)共对应\(Σk\)条边;共计约\(2n+Σk+(m+Σk)logn\)条边。当然实际上要小很多,因为一条信息中所有区间的和是\([1,n]\),每个区间对应的边数远不足\(logn\)。算这么多干嘛直接开vector能过就行啊

Code

//[POI2015]Pustynia
#include <cstdio>
#include <queue>
#include <vector>
using namespace std;
inline char gc()
{
    static char now[1<<16],*s,*t;
    if(s==t) {t=(s=now)+fread(now,1,1<<16,stdin); if(s==t) return EOF;}
    return *s++;
}
inline int read()
{
    int x=0; char ch=gc();
    while(ch<'0'||'9'<ch) ch=gc();
    while('0'<=ch&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x;
}
inline int max(int x,int y) {return x>y?x:y;}
const int N=1e5+10;
int n,n1,m;
int cnt,rt,chL[N<<1],chR[N<<1];
const int N1=4e5+10;
int edCnt;
int val[N1],a[N1],in[N1],id[N];
vector< pair<int,bool> > son[N1];
inline void edAdd(int u,int v,bool w) {edCnt++; in[v]++; son[u].push_back(make_pair(v,w));}
void bldTr(int &p,int L0,int R0)
{
    if(!p) p=++cnt;
    if(L0==R0) {id[L0]=p; return;}
    int mid=L0+R0>>1;
    bldTr(chL[p],L0,mid),bldTr(chR[p],mid+1,R0);
    edAdd(chL[p],p,0),edAdd(chR[p],p,0);
}
int optL,optR;
void trEdAdd(int p,int L0,int R0)
{
    if(optL<=L0&&R0<=optR) {edAdd(p,cnt,1); return;}
    int mid=L0+R0>>1;
    if(optL<=mid) trEdAdd(chL[p],L0,mid);
    if(mid<optR) trEdAdd(chR[p],mid+1,R0);
}
queue<int> Q;
int main()
{
    n=read(),n1=read(),m=read();
    bldTr(rt,1,n);
    for(int i=1;i<=n1;i++) {int u=id[read()]; val[u]=a[u]=read();}
    for(int i=1;i<=m;i++)
    {
        int L=read(),R=read(),k0=read();
        cnt++; int pre=L;
        while(k0--)
        {
            int x=read();
            optL=pre,optR=x-1; if(optL<=optR) trEdAdd(rt,1,n);
            pre=x+1; edAdd(cnt,id[x],0);
        }
        optL=pre,optR=R; if(optL<=optR) trEdAdd(rt,1,n);
    }
    for(int u=1;u<=cnt;u++) if(!in[u]) val[u]=max(1,a[u]),Q.push(u);
    bool noAns=false;
    while(!noAns&&!Q.empty())
    {
        int u=Q.front(); Q.pop();
        if(a[u]&&val[u]>a[u]) {noAns=true; break;}
        for(int i=0;i<son[u].size();i++)
        {
            int v=son[u][i].first,w=son[u][i].second;
            val[v]=max(val[v],val[u]+w);
            if(--in[v]==0) Q.push(v);
        }
    }
    for(int u=1;u<=cnt&&!noAns;u++) if(in[u]||val[u]>1e9) noAns=true;
    if(noAns) {puts("NIE"); return 0;}
    puts("TAK");
    for(int i=1;i<=n;i++) printf("%d ",val[id[i]]);
    puts("");
    return 0;
}

P.S.

每个数都在\([1,10^9]\)范围内! 意思是说如果你推出来有的数大于\(10^9\)就是无解,坑我半天...

原文地址:https://www.cnblogs.com/VisJiao/p/LgP3588.html

时间: 2024-10-03 07:25:35

洛谷P3588 - [POI2015]Pustynia的相关文章

BZOJ 4385 洛谷3594 POI2015 WIL-Wilcze do?y

[题解] 手残写错调了好久QAQ...... 洛谷的数据似乎比较水.. n个正整数!!这很重要 这道题是个类似two pointer的思想,外加一个单调队列维护当前区间内长度为d的子序列中元素之和的最大值. 枚举右端点,如果左端点到右端点的元素和减去区间内长为d的子序列中元素和的最大值,大于给定的P,那么就把左端点向右挪. #include<cstdio> #include<algorithm> #define N 2000010 #define rg register #defi

洛谷P3585 [POI2015]PIE

传送门 题目大意:有个n*m的格子图,要求'x'点要被染成黑色 有个a*b的印章,'x'是可以染色的印章上的点. 要求用印章去染色格子 (1)印章不可以旋转. (2)不能把墨水印到纸外面. (3)纸上的同一个格子不可以印多次. 题解:模拟 从题目中可以看出,一定要让印章的左上角对应目前n*m方 格中未染色的左上角.因为要求不能重复染色,可以每染完 一个格子就把它赋值为0.(待染色为1). 开始纯模拟,没有任何优化的代码. 加了个读入优化还是T了两个点,3000ms+ #include<iostr

洛谷 P2709 BZOJ 3781 小B的询问

题目描述 小B有一个序列,包含N个1~K之间的整数.他一共有M个询问,每个询问给定一个区间[L..R],求Sigma(c(i)^2)的值,其中i的值从1到K,其中c(i)表示数字i在[L..R]中的重复次数.小B请你帮助他回答询问. 输入输出格式 输入格式: 第一行,三个整数N.M.K. 第二行,N个整数,表示小B的序列. 接下来的M行,每行两个整数L.R. 输出格式: M行,每行一个整数,其中第i行的整数表示第i个询问的答案. 输入输出样例 输入样例#1: 6 4 3 1 3 2 1 1 3

洛谷1231 教辅的组成

洛谷1231 教辅的组成 https://www.luogu.org/problem/show?pid=1231 题目背景 滚粗了的HansBug在收拾旧语文书,然而他发现了什么奇妙的东西. 题目描述 蒟蒻HansBug在一本语文书里面发现了一本答案,然而他却明明记得这书应该还包含一份练习题.然而出现在他眼前的书多得数不胜数,其中有书,有答案,有练习册.已知一个完整的书册均应该包含且仅包含一本书.一本练习册和一份答案,然而现在全都乱做了一团.许多书上面的字迹都已经模糊了,然而HansBug还是可

洛谷教主花园dp

洛谷-教主的花园-动态规划 题目描述 教主有着一个环形的花园,他想在花园周围均匀地种上n棵树,但是教主花园的土壤很特别,每个位置适合种的树都不一样,一些树可能会因为不适合这个位置的土壤而损失观赏价值. 教主最喜欢3种树,这3种树的高度分别为10,20,30.教主希望这一圈树种得有层次感,所以任何一个位置的树要比它相邻的两棵树的高度都高或者都低,并且在此条件下,教主想要你设计出一套方案,使得观赏价值之和最高. 输入输出格式 输入格式: 输入文件garden.in的第1行为一个正整数n,表示需要种的

洛谷 P2801 教主的魔法 题解

此文为博主原创题解,转载时请通知博主,并把原文链接放在正文醒目位置. 题目链接:https://www.luogu.org/problem/show?pid=2801 题目描述 教主最近学会了一种神奇的魔法,能够使人长高.于是他准备演示给XMYZ信息组每个英雄看.于是N个英雄们又一次聚集在了一起,这次他们排成了一列,被编号为1.2.…….N. 每个人的身高一开始都是不超过1000的正整数.教主的魔法每次可以把闭区间[L, R](1≤L≤R≤N)内的英雄的身高全部加上一个整数W.(虽然L=R时并不

洛谷P1466 集合 Subset Sums

洛谷P1466 集合 Subset Sums这题可以看成是背包问题 用空间为 1--n 的物品恰好填充总空间一半的空间 有几种方案 01 背包问题 1.注意因为两个交换一下算同一种方案,所以最终 要 f [ v ] / 2 2.要开 long long 1 #include <cstdio> 2 #include <cstdlib> 3 #include <cmath> 4 #include <cstring> 5 #include <string&g

洛谷P1160 队列安排 链表

洛谷P1160 队列安排   链表 1 #include <cstdio> 2 #include <cstring> 3 #include <cmath> 4 #include <cstdlib> 5 #include <string> 6 #include <algorithm> 7 #include <iomanip> 8 #include <iostream> 9 using namespace std

洛谷 P3367 并查集模板

#include<cstdio> using namespace std; int n,m,p; int father[2000001]; int find(int x) { if(father[x]!=x) father[x]=find(father[x]); return father[x]; } void unionn(int i,int j) { father[j]=i; } int main() { scanf("%d%d",&n,&m); for