NOI2014 购票

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<vector>
using namespace std;
typedef long long LL;
inline LL read()
{
    LL x=0;bool f=0;char c=getchar();
    for (;c<'0'||c>'9';c=getchar()) f=c=='-'?1:0;
    for (;c>='0'&&c<='9';c=getchar()) x=x*10+c-'0';
    return f?-x:x;
}
const int N=200010;
int n,ll=1,to[N],next[N],head[N],fa[N];
LL p[N],q[N],lv[N],fl[N],dis[N],dp[N],_ans,_k;
int ct,id[N],ln[N<<2],rn[N<<2],_l,_r;
 
struct na
{
    LL x,y;
    na (LL a,LL b) {x=a;y=b;}
};
 
inline bool cmp(const na &a,const na &b,const na &c)
{
    return ((double)b.y-a.y)/(b.x-a.x)<((double)c.y-a.y)/(c.x-a.x);
}
 
struct pch
{
    na *h,*dh;
    int *pos,*tp,top;
    void init(int len)
    {
        top=0;
        h=(na*)malloc(sizeof(na)*(len+1));
        dh=(na*)malloc(sizeof(na)*(len+1));
        pos=(int*)malloc(sizeof(int)*(len+1));
        tp=(int*)malloc(sizeof(int)*(len+1));
    }
    void ins(const na &a)
    {
        if (top<2||cmp(h[top-2],h[top-1],a))
            *(dh++)=h[top],*(pos++)=*(tp++)=top,h[top++]=a;
        else if (!cmp(h[0],h[1],a)) 
            *(dh++)=h[1],*(pos++)=1,*(tp++)=top,h[1]=a,top=2;
        else
        {
            int l=0,r=top-2,mid;
            while (l+1<r)
            {
                mid=l+r>>1;
                if (cmp(h[mid],h[mid+1],a)) l=mid;
                else r=mid;
            }
            *(dh++)=h[r+1],*(pos++)=r+1,h[r+1]=a;
            *(tp++)=top,top=r+2;
        }
    }
    void del()
    {
        int p=*(--pos);
        h[p]=*(--dh);top=*(--tp);
    }
    inline double calk(int i)
    {
        if (!i) return -1e18;
        if (i==top) return 1e18;
        return ((double)(h+i)->y-(h+i-1)->y)/((double)(h+i)->x-(h+i-1)->x);
    }
    LL cal(LL k)
    {
        int l=0,r=top,mid;
        while (l+1<r)
        {
            if (calk(mid=l+r>>1)<=k) l=mid;
            else r=mid;
        }
        return (h+l)->y-k*(h+l)->x;
    }
}e[N<<2];
 
inline void adde(int u,int v,LL w)
{
    to[++ll]=v;next[ll]=head[u];head[u]=ll;fl[ll]=w;
}
 
int find(LL x)
{
    int l=0,r=ct-1,mid;
    while (l+1<r)
    {
        if (x<=dis[id[mid=l+r>>1]]) r=mid;
        else l=mid;
    }
    return r;
}
void build(int o,int l,int r)
{
    ln[o]=l;rn[o]=r;
    e[o].init(r-l+1);
    if (l==r) return;
    int mid=l+r>>1;
    build(o<<1,l,mid);
    build(o<<1|1,mid+1,r);
}
 
void query(int o)
{
    if (_l<=ln[o]&&rn[o]<=_r)
    {
        _ans=min(_ans,e[o].cal(_k));
        return;
    }
    int mid=ln[o]+rn[o]>>1;
    if (_l<=mid) query(o<<1);
    if (_r>mid) query(o<<1|1);
}
 
void ins(int p,LL x,LL y)
{
    int o=1,l=1,r=n,mid;
    na pt(x,y);
    e[o].ins(pt);
    while (l!=r)
    {
        if (p<=(mid=l+r>>1)) o=o<<1,r=mid;
        else o=o<<1|1,l=mid+1;
        e[o].ins(pt);
    }
}
 
void del(int x)
{
    int o=1,l=1,r=n,mid;
    e[o].del();
    while (l!=r)
    {
        if (x<=(mid=l+r>>1)) o=o<<1,r=mid;
        else o=o<<1|1,l=mid+1;
        e[o].del();
    }
}
 
void dfs(int u)
{
    id[++ct]=u;
    int lim=find(dis[u]-lv[u]);
    _l=lim;_r=ct-1;_ans=(LL)1e18;_k=p[u];
    query(1);dp[u]=_ans+p[u]*dis[u]+q[u];
    ins(ct,dis[u],dp[u]);
    for (int i=head[u];i;i=next[i]) 
        dis[to[i]]=dis[u]+fl[i],dfs(to[i]);
    del(ct--);
}
 
int main()
{
    n=read();read();
    for (int i=2;i<=n;i++)
    {
        fa[i]=read();
        adde(fa[i],i,read());
        p[i]=read();q[i]=read();lv[i]=read();
    }
    build(1,1,n);
    id[ct=1]=1;ins(1,0,0);
    for (int i=head[1];i;i=next[i])
        dis[to[i]]=fl[i],dfs(to[i]);
    for (int i=2;i<=n;i++) printf("%lld\n",dp[i]);
    return 0;
}


P2375 [NOI2014] 动物园是一道经典的动态规划题目,以下是该题的详细题意和解题思路。 【题意描述】 有两个长度为 $n$ 的整数序列 $a$ 和 $b$,你需要从这两个序列中各选出一些数,使得这些数构成一个新的序列 $c$。其中,$c$ 序列中的元素必须在原序列中严格递增。每个元素都有一个价值,你的任务是选出的元素的总价值最大。 【解题思路】 这是一道经典的动态规划题目,可以采用记忆化搜索的方法解决,也可以采用递推的方法解决。 记忆化搜索的代码如下: ```c++ #include <iostream> #include <cstdio> #include <cstring> using namespace std; const int MAXN = 1005; int dp[MAXN][MAXN], a[MAXN], b[MAXN], n; int dfs(int x, int y) { if (dp[x][y] != -1) return dp[x][y]; if (x == n || y == n) return 0; int res = max(dfs(x + 1, y), dfs(x + 1, y + 1)); if (a[x] > b[y]) { res = max(res, dfs(x, y + 1) + b[y]); } return dp[x][y] = res; } int main() { scanf("%d", &n); for (int i = 0; i < n; i++) scanf("%d", &a[i]); for (int i = 0; i < n; i++) scanf("%d", &b[i]); memset(dp, -1, sizeof(dp)); printf("%d\n", dfs(0, 0)); return 0; } ``` 其中,dp[i][j]表示选到a数组中第i个元素和b数组中第j个元素时的最大价值,-1表示未计算过。dfs(x,y)表示选到a数组中第x个元素和b数组中第y个元素时的最大价值,如果dp[x][y]已经计算过,则直接返回dp[x][y]的值。如果x==n或者y==n,表示已经遍历完一个数组,直接返回0。然后就是状态转移方程了,如果a[x] > b[y],则可以尝试选b[y],递归调用dfs(x, y+1)计算以后的最大价值。否则,只能继续遍历数组a,递归调用dfs(x+1, y)计算最大价值。最后,返回dp[0][0]的值即可。 递推的代码如下: ```c++ #include <iostream> #include <cstdio> #include <cstring> using namespace std; const int MAXN = 1005; int dp[MAXN][MAXN], a[MAXN], b[MAXN], n; int main() { scanf("%d", &n); for (int i = 0; i < n; i++) scanf("%d", &a[i]); for (int i = 0; i < n; i++) scanf("%d", &b[i]); for (int i = n - 1; i >= 0; i--) { for (int j = n - 1; j >= 0; j--) { dp[i][j] = max(dp[i + 1][j], dp[i + 1][j + 1]); if (a[i] > b[j]) { dp[i][j] = max(dp[i][j], dp[i][j + 1] + b[j]); } } } printf("%d\n", dp[0][0]); return 0; } ``` 其中,dp[i][j]表示选到a数组中第i个元素和b数组中第j个元素时的最大价值。从后往前遍历数组a和数组b,依次计算dp[i][j]的值。状态转移方程和记忆化搜索的方法是一样的。 【参考链接】 P2375 [NOI2014] 动物园:https://www.luogu.com.cn/problem/P2375
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值