• [JOISC 2020 Day4] 治疗计划


    一、题目

    点此看题

    二、解法

    我们把所有治疗方案按右端点排序,从左到右转移,设 (dp[i]) 表示把 ([1,r_i]) 都治好的最小代价。注意这个状态并没有对时间有特殊限制,我们只需要保证最后所有人能被治好就行了。

    我再进一步地解释这个状态,你可能会说这个状态很奇怪,我不一定要先治疗一个前缀,可以先从中间治起。那么你的说法是对的,所以我再三强调本题的状态对时间是没有限制的,我们可能会被思维定式所局限觉得一定要按时间顺序考虑,但是这个状态的本质是按人去考虑,治前缀只是一种表象,重要的是最后我们会把所有的人治好

    明确了上面的观点以后就可以写转移了,我们考虑右端点在 (i) 之前的一种治疗方案 (j)

    [dp[i]leftarrow dp[j]+c[i] r_j-l_i+1geq |t_i-t_j| ]

    也就是在 (i,j) 时间的空档期,中间那些新感染的人需要被治好。

    考虑优化转移,用数据结构暴力维护是不行的,因为把绝对值拆开以后会有两重偏序关系,数据结构维护不了。观察转移特点,新增代价只和终点有关,可以看成一个 ( t dijkstra) 模型,每次选取最小的还未松弛的 (f[j]),去找能访问到的 (i),不难发现从 (j) 转移到 (i) 是最优的,所以每个点只会被访问一次。

    那么以时间排序就能把绝对值拆掉,然后搞个线段树维护即可,势能分析可知时间复杂度 (O(nlog n))

    三、总结

    当拿到这道题的时候,可能潜意识地就以时间为 (dp) 主体,这是难以发现的思维定式。

    所以选定 (dp) 主体很重要,然后我在此声明,实践是检验 (dp) 的唯一标准,我们只需要看他能不能考虑到所有情况,能不能转移,而不用管他是否奇怪。

    此外本题提供了一种新的优化思路,当转移新增权值之和终点有关时,可以用 ( t dijkstra) 的贪心来做,这样就有每个点只会被更新一次的美妙性质。

    #include <cstdio>
    #include <iostream>
    #include <algorithm>
    #include <queue>
    using namespace std;
    const int M = 100005;
    #define int long long
    const int inf = 1e18;
    int read()
    {
    	int x=0,f=1;char c;
    	while((c=getchar())<'0' || c>'9') {if(c=='-') f=-1;}
    	while(c>='0' && c<='9') {x=(x<<3)+(x<<1)+(c^48);c=getchar();}
    	return x*f;
    }
    int n,m,ans,f[M],mn[4*M][2];
    struct node
    {
    	int t,l,r,c;
    	bool operator < (const node &b) const
    	{
    		return t<b.t;
    	}
    }a[M];
    struct data
    {
    	int u,c;
    	bool operator < (const data &b) const
    	{
    		return c>b.c;
    	}
    };priority_queue<data> q;
    void up(int i)
    {
    	mn[i][0]=min(mn[i<<1][0],mn[i<<1|1][0]);
    	mn[i][1]=min(mn[i<<1][1],mn[i<<1|1][1]);
    }
    void ins(int i,int l,int r,int id,int x,int v)
    {
    	if(l==r)
    	{
    		mn[i][0]=x-v;
    		mn[i][1]=x+v;
    		return ;
    	}
    	int mid=(l+r)>>1;
    	if(mid>=id) ins(i<<1,l,mid,id,x,v);
    	else ins(i<<1|1,mid+1,r,id,x,v);
    	up(i);
    }
    void ask(int i,int l,int r,int L,int R,int op,int v,int w)
    {
    	if(l>R || L>r || mn[i][op]>v) return ;
    	if(l==r)
    	{
    		f[l]=w+a[l].c;q.push(data{l,f[l]});
    		mn[i][0]=mn[i][1]=inf;
    		return ;
    	}
    	int mid=(l+r)>>1;
    	ask(i<<1,l,mid,L,R,op,v,w);
    	ask(i<<1|1,mid+1,r,L,R,op,v,w);
    	up(i);
    }
    signed main()
    {
    	n=read();m=read();ans=inf;
    	for(int i=1;i<=m;i++)
    	{
    		f[i]=inf;
    		a[i].t=read();a[i].l=read();
    		a[i].r=read();a[i].c=read();
    	}
    	sort(a+1,a+1+m);
    	for(int i=1;i<=m;i++)
    	{
    		if(a[i].l==1)
    		{
    			f[i]=a[i].c;
    			q.push(data{i,f[i]});
    			ins(1,1,m,i,inf,0);
    		}
    		else ins(1,1,m,i,a[i].l,a[i].t);
    	}
    	while(!q.empty())
    	{
    		int u=q.top().u;q.pop();
    		ask(1,1,m,1,u-1,0,a[u].r-a[u].t+1,f[u]);
    		ask(1,1,m,u+1,m,1,a[u].r+a[u].t+1,f[u]);
    	}
    	for(int i=1;i<=m;i++)
    		if(a[i].r==n)
    			ans=min(ans,f[i]);
    	if(ans==inf) puts("-1");
    	else printf("%lld
    ",ans);
    }
    
  • 相关阅读:
    专业实训项目需求分析
    2015年秋季个人阅读计划
    场景调研
    二维数组最大连通子数组
    单元测试
    《大道至简——软件工程实践者的思想》阅读笔记之三
    《大道至简——软件工程实践者的思想》阅读笔记之二
    人机交互-输入法使用评价
    第一阶段个人总结10
    第一阶段个人总结09
  • 原文地址:https://www.cnblogs.com/C202044zxy/p/15013140.html
Copyright © 2020-2023  润新知