BJ 集训8 测试 B 新访问计划

本文针对APIO2010巡逻问题提供了一种改进算法,通过增加可修建道路数量至2e4,并引入跳跃成本的概念,利用动态规划优化路径选择策略。详细介绍了算法实现过程及核心代码。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

http://www.elijahqi.win/archives/2757
题意:参考apio2010巡逻 只不过这次将可以修建的道路增加到2e4

那么考虑增加枚举每次跳跃的花费 那么可以知道随着跳跃的花费增加我的跳跃次数是递减的

所以那么考虑显然如果我跳跃的答案更小的话 当我跳跃次数<=k的时候我一定尽量跳跃更多的即可 那么设dp数组dp[x][0/1]分别表示当前我处在x节点 我没有走回/走回当前节点 我的最优答案是多少 在满足这个最优答案的时候我需要跳跃几次才可以满足条件 注意当结尾的时候我可能存在一些问题 比如我跳跃的时候代价可能正好等于我走边的边权的时候 那么我可能用跳跃来代替了我这个走边权的机会那么因为我跳跃固定在c 所以我就会多增加一些费用 所以应该恰好将k条边多增加的代价卡下来

#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 22000
#define pa pair<int,int>
#define mp make_pair
#define fi first
#define se second
#define inf 0x3f3f3f3f
using namespace std;
inline int read(){
    int x=0;char ch=getchar();
    while(ch<'0'||ch>'9') ch=getchar();
    while(ch<='9'&&ch>='0') x=x*10+ch-'0',ch=getchar();
    return x;
}
struct node{
    int y,next,z;
}data[N<<1];
pa dp[N][2];
int h[N],n,k,c,cost,num;
inline void dfs(int x,int fa){
    dp[x][0].fi=dp[x][1].fi=dp[x][0].se=dp[x][1].se=0;
    for (int i=h[x];i;i=data[i].next){
        int y=data[i].y,z=data[i].z;if (y==fa) continue;pa tmp0,tmp1;dfs(y,x);
        tmp0=min(mp(dp[x][1].fi+dp[y][0].fi,dp[x][1].se+dp[y][0].se),mp(dp[x][0].fi+dp[y][1].fi+z,dp[x][0].se+dp[y][1].se));
        tmp1=min(mp(dp[x][1].fi+dp[y][1].fi+z,dp[x][1].se+dp[y][1].se),mp(dp[x][0].fi+dp[y][0].fi+cost,dp[x][0].se+dp[y][0].se+1));
        dp[x][0]=tmp0;dp[x][1]=tmp1;
    }
}
inline bool check(int x){
    cost=x;dfs(1,1);return dp[1][1].se<=k;
}
int main(){
    freopen("t2.in","r",stdin);
    while(~scanf("%d",&n)){
        k=read();c=read();
        memset(h,0,sizeof(int)*(n+10));num=0;int sum=0;
        for (int i=1;i<n;++i) {
            int x=read()+1,y=read()+1,z=read();sum+=z;
            data[++num].y=y;data[num].next=h[x];data[num].z=z;h[x]=num;
            data[++num].y=x;data[num].next=h[y];data[num].z=z;h[y]=num;
        }cost=c;dfs(1,1);if (dp[1][1].se<=k) {printf("%d\n",dp[1][1].fi+sum);continue;}
        int l=c+1,r=1e9,mid;
        while(l<=r){
            mid=l+r>>1;
            if (check(mid)) r=mid-1;else l=mid+1;
        }check(l);printf("%d\n",sum+dp[1][1].fi-(cost-c)*k);
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值