【BZOJ-1097】旅游景点atr SPFA + 状压DP

1097: [POI2007]旅游景点atr

Time Limit: 30 Sec  Memory Limit: 357 MB
Submit: 1531  Solved: 352
[Submit][Status][Discuss]

Description

  FGD想从成都去上海旅游。在旅途中他希望经过一些城市并在那里欣赏风景,品尝风味小吃或者做其他的有趣
的事情。经过这些城市的顺序不是完全随意的,比如说FGD不希望在刚吃过一顿大餐之后立刻去下一个城市登山,
而是希望去另外什么地方喝下午茶。幸运的是,FGD的旅程不是既定的,他可以在某些旅行方案之间进行选择。由于
FGD非常讨厌乘车的颠簸,他希望在满足他的要求的情况下,旅行的距离尽量短,这样他就有足够的精力来欣赏风
景或者是泡MM了^_^.整个城市交通网络包含N个城市以及城市与城市之间的双向道路M条。城市自1至N依次编号,道
路亦然。没有从某个城市直接到它自己的道路,两个城市之间最多只有一条道路直接相连,但可以有多条连接两个
城市的路径。任意两条道路如果相遇,则相遇点也必然是这N个城市之一,在中途,由于修建了立交桥和下穿隧道
,道路是不会相交的。每条道路都有一个固定长度。在中途,FGD想要经过K(K<=N-2)个城市。成都编号为1,上海
编号为N,而FGD想要经过的N个城市编号依次为2,3,…,K+1.举例来说,假设交通网络如下图。FGD想要经过城市2,3,
4,5,并且在2停留的时候在3之前,而在4,5停留的时候在3之后。那么最短的旅行方案是1-2-4-3-4-5-8,总长度为1
9。注意FGD为了从城市2到城市4可以路过城市3,但不在城市3停留。这样就不违反FGD的要求了。并且由于FGD想要
走最短的路径,因此这个方案正是FGD需要的。

Input

  第一行包含3个整数N(2<=N<=20000),M(1<=M<=200000),K(0<=K<=20),意义如上所述。

Output

  只包含一行,包含一个整数,表示最短的旅行距离。

Sample Input

8 15 4
1 2 3
1 3 4
1 4 4
1 6 2
1 7 3
2 3 6
2 4 2
2 5 2
3 4 3
3 6 3
3 8 6
4 5 2
4 8 6
5 7 4
5 8 6
3
2 3
3 4
3 5

Sample Output

19

HINT

 上面对应于题目中给出的例子。

Source

Solution

SPFA预处理,状压DP

SPFA预处理出K+1个到各点的距离,状压DP,$dp[i][j]$表示当前状态为$i$(访问过的点的集合),现在在$j$点

转移很好转移..

PS:卡常数...我都以为我要T了...

Code

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
int read()
{
    int x=0,f=1; char ch=getchar();
    while (ch<‘0‘ || ch>‘9‘) {if (ch==‘-‘) f=-1; ch=getchar();}
    while (ch>=‘0‘ && ch<=‘9‘) {x=x*10+ch-‘0‘; ch=getchar();}
    return x*f;
}
#define maxn 20010
#define maxm 400010
int N,M,K,C,dis[maxn],d[30][maxn],ans; bool visit[maxn];
struct EdgeNode{int next,to,len,from;}edge[maxm];
int head[maxn],cnt;
void add(int u,int v,int w)
{
    cnt++;
    edge[cnt].to=v;edge[cnt].from=u;edge[cnt].len=w;edge[cnt].next=head[u];head[u]=cnt;
}
void insert(int u,int v,int w) {add(u,v,w); add(v,u,w);}

#define inf 0x7fffffff
void SPFA(int S)
{
    queue<int>q; memset(visit,0,sizeof(visit));
    for (int i=1; i<=N; i++) dis[i]=inf;
    q.push(S); visit[S]=1; dis[S]=0;
    while (!q.empty())
        {
            int now=q.front(); q.pop(); visit[now]=0;
            for (int i=head[now]; i; i=edge[i].next)
                if (dis[edge[i].to]>dis[now]+edge[i].len)
                    {
                        dis[edge[i].to]=dis[now]+edge[i].len;
                        if (!visit[edge[i].to])
                            q.push(edge[i].to),visit[edge[i].to]=1;
                    }
        }
    for (int i=1; i<=N; i++) d[S][i]=dis[i];
}
int Min(int a,int b)
{
    if (a==-1) return b;
    if (a>b) return b; else return a;
}
int dp[1<<21][30],zt[maxn];
void DP()
{
    for (int i=0; i<=(1<<K)-1; i++)
        for (int j=1; j<=K+1; j++)
            if (dp[i][j]!=-1)
                for (int k=2; k<=K+1; k++)
                    if ((i&zt[k])==zt[k])
                        dp[i|(1<<(k-2))][k]=Min(dp[i|(1<<(k-2))][k],dp[i][j]+d[j][k]);
}
int main()
{
    N=read(),M=read(),K=read();
    for (int u,v,w,i=1; i<=M; i++)
        u=read(),v=read(),w=read(),insert(u,v,w);
    for (int i=1; i<=K+1; i++) SPFA(i);
    C=read();
    for (int x,y,i=1; i<=C; i++)
        x=read(),y=read(),zt[y]|=(1<<(x-2));
    memset(dp,-1,sizeof(dp)); dp[0][1]=0;
    DP(); ans=0x7fffffff;
    for (int i=1; i<=K+1; i++)
        if (dp[(1<<K)-1][i]!=-1) ans=Min(dp[(1<<K)-1][i]+d[i][N],ans);
    printf("%d",ans);
    return 0;
}

开了两个DeBug,第一次忘删了OLE,删了一个还剩一个又WA了...然后卡时..

时间: 2024-12-28 17:59:43

【BZOJ-1097】旅游景点atr SPFA + 状压DP的相关文章

BZOJ 1097 POI2007 旅游景点atr SPFA+状压DP

题目大意:给定一张图,要求从第一个点出发,按照某个拓扑序遍历2~k+1的所有节点,然后到达n,求最短路径 首先将所有关键点之间的最短路用SPFA求出来 然后状压DP 令f[state][p]表示已经走过的点集为state,将要走到p点的最短路 记忆化搜索就行了- - 标准卡时过- - #include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #define M

[AC自动机+spfa+状压dp] hdu 3247 Resource Archiver

题意: 给n个本源串,m个病毒串 求一个最多的长度的单词包含所有的本源串并不包含任意一个病毒串 串均为01串 思路: 只有10个本源串 一开始想的是直接建立完trie图 然后在图上直接spfa 结果发现 dis[60005][1030] 超内存了 这个时候就要想到 其实只有节点的mark值大于0的节点是我们需要用的 就是那些含有状压权值的节点 那么我们先记录下这些节点是哪些 然后发现其实这些不到100个节点 所以跑100遍spfa 求出两两之间的最短路 然后用这个距离 去状压dp 数组就成了 d

[BZOJ 1072] [SCOI2007] 排列perm 【状压DP】

题目链接:BZOJ 1072 这道题使用 C++ STL 的 next_permutation() 函数直接暴力就可以AC .(使用 Set 判断是否重复) 代码如下: #include <iostream> #include <cstdio> #include <cstdlib> #include <cstring> #include <algorithm> #include <cmath> #include <set>

BZOJ 1087: [SCOI2005]互不侵犯King( 状压dp )

简单的状压dp... dp( x , h , s ) 表示当前第 x 行 , 用了 h 个 king , 当前行的状态为 s . 考虑转移 : dp( x , h , s ) = ∑ dp( x - 1 , h - cnt_1( s ) , s' ) ( s and s' 两行不冲突 , cnt_1( s ) 表示 s 状态用了多少个 king ) 我有各种预处理所以 code 的方程和这有点不一样 ------------------------------------------------

BZOJ 3446: [Usaco2014 Feb]Cow Decathlon( 状压dp )

水状压dp. dp(x, s) = max{ dp( x - 1, s - {h} ) } + 奖励(假如拿到的) (h∈s). 时间复杂度O(n * 2^n) ---------------------------------------------------------------------------------- #include<bits/stdc++.h> #define rep(i, n) for(int i = 0; i < n; ++i) #define clr(x

BZOJ 1879 [Sdoi2009]Bill的挑战 ——状压DP

本来打算好好写写SDOI的DP题目,但是忒难了, 太难了,就写的这三道题仿佛是可做的. 生在弱省真是兴奋. 这题目直接状压,f[i][j]表示匹配到i,状态集合为j的方案数,然后递推即可. #include <cstdio> #include <cstring> #include <iostream> #include <algorithm> using namespace std; #define F(i,j,k) for (int i=j;i<=k

BZOJ 1076 SCOI2008 奖励关 期望状压DP

题目大意:给定k次弹出宝物的机会,每次随机弹出n种宝物的机会,如果吃过这种宝物的所有前提宝物就可以吃这种宝物,求最优策略的期望得分 看到数据范围果断状压DP- - 不看数据范围害死人- - 至于吃还是不吃 这是个问题 对于这种最优策略的期望DP 我们一般都是从后往前推 枚举每次出现宝物 枚举此时的状态 枚举宝物是哪种 如果当前的宝物可以吃 就在吃与不吃的后继状态中选择最大值加到当前状态上 如果当前的宝物不能吃 只能选择不吃的后继状态加到当前状态上 最后输出f[1][0]就是答案 #include

BZOJ 4197 NOI 2015 寿司晚宴 状压DP

4197: [Noi2015]寿司晚宴 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 694  Solved: 440[Submit][Status][Discuss] Description 为了庆祝 NOI 的成功开幕,主办方为大家准备了一场寿司晚宴.小 G 和小 W 作为参加 NOI 的选手,也被邀请参加了寿司晚宴. 在晚宴上,主办方为大家提供了 n?1 种不同的寿司,编号 1,2,3,-,n?1,其中第 i 种寿司的美味度为 i+1 (即

bzoj 4197: [Noi2015]寿司晚宴【状压dp】

一个数内可能多个的质因数只有小于根号n的,500内这样的数只有8个,所以考虑状压 把2~n的数处理出小于根号500的质因数集压成s,以及大质数p(没有就是1),然后按p排序 根据题目要求,拥有一个质因数的只能给一个人,所以排序后能给一个人的大质数就是一个区间 然后设f[s1][s2]为一人选s1,另一人选s2的方案数,注意这里的s只压了小于根号500的八个质数 设g[0/1][s1][s2]为一人选s1,另一人选s2的,当前枚举的大质数给小G/小W的方案数 正常转移即可 然后注意把g转到f上时应