环游衢州
(walk.pas/c/cpp)
【问题描述】
Jason 想要带着他新叫的朋友环游 QZ, 于是他把 QZ 划分成了 n 个地方,已知这 n 个地方由 m 条双向边相连接,并且他把 n 个地方分别赋予了一个 happy 值,他希望按照 happy 值严格递减的方法去游览QZ, 这样他可以把一些最美的地方介绍给他朋友。 然而 Jason 的心态不太好,所以他有 k 次从小向大走的机会,问他最多能够带他的朋友游览几个地方,如果能够游览的地方数超过 10^9+7则输出“ infinity”
【输入格式】
第一行两个数 n,m, k 表示地方数、边数和反向次数
第二行 n 个数表示 n 个地点的 happy 值
然后 m 行每行两个数 a,b 表示 a, b 之间有一条无向边
【输出格式】
一行一个数表示最多能够游览的地方数或者“ infinity”
【输入样例1】
1 3 1
2147483647
1 1
1 1
1 1
【输出样例1】
2
【输入样例2】
10 30 20
720472486 1616090782 1659830133 376600248 1485380712 819748825 1072214931 1957734249 95302927 622052677
2 6
1 6
10 2
8 4
8 2
9 6
1 4
9 7
1 4
4 7
7 1
4 5
10 6
7 7
8 5
7 8
6 2
7 2
3 9
6 8
3 4
9 8
8 8
4 5
8 7
5 6
2 2
9 7
3 1
9 1
【输出样例2】
105
【样例解释】
在 1 号点的自环上走一步即可到达两个点。
【 数据说明】
对于 10%的数据: 1<=n<=10, 1<=m<=30, k=0
对于 30%的数据: 1<=n<=500, 1<=m<=2500, k=0
对于 60%的数据: 1<=n<=1000, m<=3000, k<=20
对于 100%的数据: 1<=n<=10^5, 1<=m<=2*10^5,k<=80, happy 值小于等于 2^31-1
【提示】
一个地点经过两次按两个地点算。
pascal 请勿使用 readln。
【题目分析】
这道题是典型的将求最长单调序列套在图的模型上,那么先将这个模型。
首先为了方便推,a[i]需要排序,但是这样会有一个严重的问题就是原来的图会被打乱,所以在实际操作中不能这样做,正确的做法应该另设一个数组id[i],表示排名为i的点在图中的具体位置,初始时id[i]=i,排序的时候排的标准是a[i],但是调动的是id[i],这样就可以不会改变原来的图了。
接下来就是DP的推法(以最长降为例):
f[i]同在序列中的定义,就是以i为结尾的最长降的路。转移时为了能得到更优的答案所以是从大到小推(这就是为什么需要进行排序),比较简单,代码如下:
for (int i=n;i;i--)
for (int j=lnk[id[i]];j;j=nxt[j])
if (a[id[i]]<a[son[j]]) f[id[i]]=max(f[id[i]],f[son[j]]+1);
回到这道题,不过就是添加了一个可以逆走的次数,那么就把逆走的次数带到转移中,f[T][i]表示使用了T次逆走,目前在i点的距离,转移方程就比较简单,详细看代码。
【复杂度】
时间:O(nT)(实际偏大); 空间:O(nT);
【代码】
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
using namespace std;
#define maxn 100005
#define maxe 400005
#define tt 1000000007
int n,e,k,tot,ans,a[maxn],id[maxn],son[maxe],nxt[maxe],lnk[maxn],f[85][maxn];
inline void readi(int &x){
x=0; char ch=getchar();
while ('0'>ch||ch>'9') ch=getchar();
while ('0'<=ch&&ch<='9') {x=x*10+ch-'0'; ch=getchar();}
}
void _add(int x,int y){
tot++; son[tot]=y; nxt[tot]=lnk[x]; lnk[x]=tot;
}
inline int cmp_id(int i,int j){return a[i]<a[j];}
void _init()
{
freopen("walk.in","r",stdin);
freopen("walk.out","w",stdout);
readi(n); readi(e); readi(k);
for (int i=1;i<=n;i++) {readi(a[i]); id[i]=i;}
tot=0;
memset(son,0,sizeof(son));
memset(nxt,0,sizeof(nxt));
memset(lnk,0,sizeof(lnk));
for (int i=1;i<=e;i++)
{
int x,y; readi(x); readi(y);
_add(x,y); _add(y,x);
}
sort(id+1,id+n+1,cmp_id);
}
void _solve()
{
memset(f,0,sizeof(f));
for (int i=n;i;i--)
for (int j=lnk[id[i]];j;j=nxt[j])
if (a[son[j]]>a[id[i]])
f[0][id[i]]=max(f[0][id[i]],f[0][son[j]]+1);
for (int t=1;t<=k;t++)
for (int i=n;i;i--)
for (int j=lnk[id[i]];j;j=nxt[j])
if (a[id[i]]<a[son[j]]) f[t][id[i]]=max(f[t][id[i]],f[t][son[j]]+1);
else f[t][id[i]]=max(f[t][id[i]],f[t-1][son[j]]+1);
ans=0;
for (int i=1;i<=n;i++)
for (int j=0;j<=k;j++) ans=max(ans,f[j][i]);
printf("%d\n",ans+1);
}
int main()
{
_init();
_solve();
return 0;
}
PS:
如果各位神犇发现本博客有误,请一定及时告诉作者我这个菜鸟噢,新手上路,请多指教。