题意:给定一个长为n的序列,每个位置可以选择取或不取,要求构造方案使得:

1.对于前M1个约束,区间【L,R】内取的数量必须严格不少于K

2.对于后M2个约束,区间【L,R】外取的数量必须严格不少于K

3.满足所有M1+M2个约束的前提下使得取的个数最少,输出最少取的个数

n,M1,M2<=3e3

思路:做法一:

特殊的SPFA判负环的技巧见https://www.cnblogs.com/myx12345/p/6212893.html

大致说来就是用栈和初始置0两个地方

 #include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned int uint;
typedef unsigned long long ull;
typedef long double ld;
typedef pair<int,int> PII;
typedef pair<ll,ll> Pll;
typedef vector<int> VI;
typedef vector<PII> VII;
typedef pair<ll,ll>P;
#define N 500010
#define M 1000000
#define INF 1e9
#define fi first
#define se second
#define MP make_pair
#define pb push_back
#define pi acos(-1)
#define mem(a,b) memset(a,b,sizeof(a))
#define rep(i,a,b) for(int i=(int)a;i<=(int)b;i++)
#define per(i,a,b) for(int i=(int)a;i>=(int)b;i--)
#define lowbit(x) x&(-x)
#define Rand (rand()*(1<<16)+rand())
#define id(x) ((x)<=B?(x):m-n/(x)+1)
#define ls p<<1
#define rs p<<1|1
#define fors(i) for(auto i:e[x]) if(i!=p) const int MOD=1e9+,inv2=(MOD+)/;
double eps=1e-;
int dx[]={-,,,};
int dy[]={,,-,}; struct edge
{
int x,y,z;
}a[N],b[N]; int head[N],vet[N],nxt[N],inq[N],dis[N],stk[N],len[N],t[N],
n,m1,m2,tot; int read()
{
int v=,f=;
char c=getchar();
while(c<||<c) {if(c=='-') f=-; c=getchar();}
while(<=c&&c<=) v=(v<<)+v+v+c-,c=getchar();
return v*f;
} ll readll()
{
ll v=,f=;
char c=getchar();
while(c<||<c) {if(c=='-') f=-; c=getchar();}
while(<=c&&c<=) v=(v<<)+v+v+c-,c=getchar();
return v*f;
} void add(int a,int b,int c)
{
nxt[++tot]=head[a];
vet[tot]=b;
len[tot]=c;
head[a]=tot;
} void build(int k)
{
//printf("k=%d\n",k);
tot=;
rep(i,,n) head[i]=;
rep(i,,n-)
{
add(i+,i,);
add(i,i+,);
}
rep(i,,m1) add(a[i].y,a[i].x-,-a[i].z);
rep(i,,m2) add(b[i].x-,b[i].y,k-b[i].z);
add(,n,k);
add(n,,-k);
} int isok()
{
rep(i,,n) inq[i]=dis[i]=t[i]=;
int top=;
rep(i,,n)
{
top++;
stk[top]=i;
inq[i]=;
t[]=i;
}
while(top)
{
int u=stk[top];
top--;
inq[u]=;
int e=head[u];
while(e)
{
int v=vet[e];
if(dis[u]+len[e]<dis[v])
{
dis[v]=dis[u]+len[e];
if(!inq[v])
{
top++;
stk[top]=v;
inq[v]=;
t[v]++;
if(t[v]==n+) return ;
}
}
e=nxt[e];
}
}
return ;
} void solve()
{
n=read(),m1=read(),m2=read();
rep(i,,m1) a[i].x=read(),a[i].y=read(),a[i].z=read();
rep(i,,m2) b[i].x=read(),b[i].y=read(),b[i].z=read();
int l=,r=n,last=n;
while(l<=r)
{
int mid=(l+r)>>;
build(mid);
if(isok()){last=mid; r=mid-;}
else l=mid+;
}
printf("%d\n",last);
} int main()
{
//freopen("1.in","r",stdin);
int cas=read();
while(cas--) solve();
return ;
}

最新文章

  1. 四HttpServletResponse常见应用——生成验证码
  2. where, group by, having
  3. boost在linux下的编译和使用
  4. 【Linux 命令】Linux系统下强制用户下线——who,pkill
  5. 398. Random Pick Index
  6. [iOS]リソースファイルの取得方法
  7. sql中在查询语句中加判断,控制输出的内容
  8. HDU 2079-课程时间(生成函数)
  9. VBS 选择文件夹框
  10. 201521123034《Java程序设计》第八周学习总结
  11. [10] 过滤器 Filter
  12. bugku web 变量1
  13. Python中常见的序列及其函数
  14. &lt;Android基础&gt;(二) Activity Part 1
  15. HDU6341 Let Sudoku Rotate (杭电多校4J)
  16. python之迭代器、生成器与面向过程编程
  17. C memset
  18. AI大牛阿里VP贾扬清
  19. C# 中对COOKIES的操作
  20. 第五周加分题--mybash的实现

热门文章

  1. Dubbo从入门到精通
  2. [转帖]2018年全球ERP软件行业市场规模与发展趋势分析 云ERP将兴起【组图】
  3. 小记---------网页采集之Jsoup
  4. java集群技术(转)
  5. springboot -- 2.0版本自定义ReidsCacheManager的改变
  6. php Excel 导入
  7. 企业面试题|最常问的MySQL面试题集合(三)
  8. 1-Kubernetes基本概念
  9. RabbitMQ从安装到使用
  10. linux centos中安装flash player