【NOI2006】聰明的導遊
【NOI2006】聰明的導遊
題面
洛谷
題目描述
小佳最近迷上了導遊這個工作,一天到晚想著帶遊客參觀各處的景點。正好 M 市在舉行 NOI,來參觀的人特別的多。不少朋友給小佳介紹了需要導遊的人。
M 市有nnn 個著名的景點,小佳將這些景點從111 至nnn 編號。有一些景點之間存在雙向的路。小佳可以讓遊客們在任何一個景點集合,然後帶著他們參觀,最後也可以在任何一個景點結束參觀。不過,來參觀的遊客們都不願去已經參觀過的地方。所以,小佳不能帶遊客們經過同一個景點兩次或兩次以上。
小佳希望你幫助他設計一個方案, 走可行的路線, 帶遊客們參觀盡可能多的地方。
輸入輸出格式
輸入格式:
輸入文件為 guide1.in~guide10.in,第一行為兩個整數n,mn,mn,m ,分別表示景點數和路的條數。接下來mmm 行,每行兩個整數a,ba,ba,b ,表示景點aaa 和景點bbb 之間有一條雙向路。
輸出格式:
你需要將答案輸出到 guide1.out~guide10.out 中,guide?.out 為對應 guide?.in
的答案。輸出的第一行為ppp ,表示你能找到的路徑所經過的景點個數。接下來ppp 行,每行一個整數,按順序表示你所找到的路徑上的每一個景點。
題解
數據什麽的到洛谷上去下載
也到洛谷上提交
回歸題目
看到題答題,一般都是怎麽想著怎麽亂搞
於是我YY了一種亂搞方法:
每次從任意一個點開始
訪問之前計算出來的答案中最大的那個兒子
知道走不了為止
然後把這個鏈的長度設為當前為止的答案
然後重復多次這個過程,
結合1,2兩個點的打爆力,第10個點的樹形\(dp\)
於是拿到了\(86\)
附上我一頓亂搞的代碼(因為代碼很亂搞,所以代碼也很亂)
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define RG register
#define MAX 11000
inline int read()
{
RG int x=0,t=1;RG char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
int n,m;
bool vis[MAX];
struct Line{int v,next;}e[MAX<<2];
int h[MAX],cnt=1;
vector<int> ans,st;
int du[MAX];
int f[MAX],ret=0,dep[MAX],fa[MAX];
inline void Add(int u,int v){e[cnt]=(Line){v,h[u]};h[u]=cnt++;}
void init()
{
memset(h,0,sizeof(h));cnt=1;
memset(f,0,sizeof(f));memset(dep,0,sizeof(dep));
memset(fa,0,sizeof(fa));ret=0;
for(int i=1;i<=n;++i)du[i]=i;
}
int getf(int x){return x==du[x]?x:du[x]=getf(du[x]);}
void dp(int u,int ff)
{
f[u]=1;dep[u]=dep[ff]+1;fa[u]=ff;
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;
if(v==ff)continue;
dp(v,u);
ret=max(ret,f[u]+f[v]);
f[u]=max(f[u],f[v]+1);
}
}
int bfn[MAX];
int bfs(int S)
{
memset(vis,0,sizeof(vis));
queue<int> Q;
Q.push(S);
int ret,tim=0;
while(!Q.empty())
{
int u=Q.front();Q.pop();
ret=u;vis[u]=true;bfn[u]=++tim;
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;
if(vis[v])continue;
vis[v]=true;Q.push(v);
}
}
return ret;
}
vector<int> ans1,ans2;
void jump(int u,int v)
{
ans1.clear();ans2.clear();
if(dep[u]<dep[v])swap(u,v);
while(dep[u]!=dep[v])ans1.push_back(u),u=fa[u];
while(u!=v)
{
ans1.push_back(u);ans2.push_back(v);
u=fa[u];v=fa[v];
}
ans1.push_back(v);
for(int i=ans2.size()-1;i>=0;--i)ans1.push_back(ans2[i]);
if(ans1.size()>ans.size())ans=ans1;
}
struct edge{int u,v;}E[MAX<<1];
void Build()
{
for(int i=1;i<=m;++i)
{
int u=E[i].u,v=E[i].v;
if(getf(u)==getf(v))continue;
Add(u,v);Add(v,u);
du[getf(u)]=getf(v);
}
}
int son[MAX],sum,Ans[MAX];
bool cmp(int a,int b){if(Ans[a]!=Ans[b])return Ans[a]<Ans[b];else return rand()%2;}
int main()
{
srand(time(NULL));
freopen("guide4.in","r",stdin);
freopen("guide4.out","w",stdout);
n=read();m=read();
for(int i=1;i<=m;++i)
{
int u=read(),v=read();
Add(u,v);Add(v,u);
}
for(int T=1;T<=1000;++T)
for(int i=1;i<=n;++i)
{
int u=i;ans1.clear();
memset(vis,0,sizeof(vis));
while(233)
{
sum=0;int nt=0;vis[u]=true;ans1.push_back(u);
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;if(vis[v])continue;
if(Ans[v]>Ans[nt])nt=v;
}
if(!nt)break;
u=nt;
}
if(ans1.size()>ans.size())ans=ans1;
if(rand()%20==0)
Ans[i]=ans1.size();
}
printf("%d\n",(int)(ans.size()));
for(int i=0;i<(int)(ans.size());++i)printf("%d\n",ans[i]);
return 0;
}
這樣跑完,
\(1,2,5,7,8,9,10\)都是正確的解了
但是因為我這份代碼不涉及隨機化(盡管我後來加了)
所以答案沒有辦法更優了
於是直接搜別人的題解。。。
和我一開始的想法比較類似:隨機找一棵生成樹,求直徑
但是我就是真正的隨機一棵生成樹。。。。
那麽如果找生成樹呢?
貪心的想,如果能夠把樹直接拉成一條鏈,一定更優
一條鏈有什麽特征?
每個點的入度/出度很小
所以,我們按照度數來隨機。
大致的過程就是:
首先隨機一個點作為根節點
然後檢查所有沒有被訪問過兒子
找到其中出度最小的兒子
然後繼續遞歸處理這個兒子
回溯到這個兒子之後,如果還有兒子沒有訪問過,直接訪問
構建出樹之後在上面跑\(bfs\)求直徑即可
代碼
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define RG register
#define MAX 11000
inline int read()
{
RG int x=0,t=1;RG char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
int n,m;
bool vis[MAX];
struct Line{int v,next;}e[MAX<<2];
int h[MAX],cnt=1;
vector<int> ans;
int du[MAX],dg[MAX];
int f[MAX],ret=0,dep[MAX],fa[MAX];
inline void Add(int u,int v){e[cnt]=(Line){v,h[u]};h[u]=cnt++;}
void init()
{
memset(h,0,sizeof(h));cnt=1;
memset(f,0,sizeof(f));memset(dep,0,sizeof(dep));
memset(fa,0,sizeof(fa));ret=0;
for(int i=1;i<=n;++i)dg[i]=du[i];
}
int bfs(int S)
{
memset(vis,0,sizeof(vis));
queue<int> Q;
Q.push(S);
int ret;
while(!Q.empty())
{
int u=Q.front();Q.pop();
ret=u;vis[u]=true;
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;
if(vis[v])continue;
vis[v]=true;Q.push(v);
}
}
return ret;
}
vector<int> ans1,ans2;
void jump(int u,int v)
{
ans1.clear();ans2.clear();
if(dep[u]<dep[v])swap(u,v);
while(dep[u]!=dep[v])ans1.push_back(u),u=fa[u];
while(u!=v)
{
ans1.push_back(u);ans2.push_back(v);
u=fa[u];v=fa[v];
}
ans1.push_back(v);
for(int i=ans2.size()-1;i>=0;--i)ans1.push_back(ans2[i]);
if(ans1.size()>ans.size())ans=ans1;
}
vector<int> Ed[MAX];
void Build(int u)
{
vis[u]=true;
for(int i=0;i<Ed[u].size();i++)--dg[Ed[u][i]];
int nt=0;dg[0]=1e9;
for(int i=0;i<Ed[u].size();i++)
{
int v=Ed[u][i];
if(vis[v])continue;
if(dg[nt]>dg[v])nt=v;
else if(dg[nt]==dg[v]&&rand()%2)nt=v;
}
if(nt)Add(u,nt),Add(nt,u),Build(nt);
for(int i=0;i<Ed[u].size();i++)
if(!vis[Ed[u][i]])Add(u,Ed[u][i]),Add(Ed[u][i],u),Build(Ed[u][i]);
}
void dfs(int u,int ff)
{
dep[u]=dep[ff]+1;fa[u]=ff;
for(int i=h[u];i;i=e[i].next)
if(e[i].v!=ff)dfs(e[i].v,u);
}
int main()
{
srand(time(NULL));
freopen("guide6.in","r",stdin);
freopen("guide6.out","w",stdout);
n=read();m=read();
for(int i=1;i<=m;++i)
{
int u=read(),v=read();
Ed[u].push_back(v);
Ed[v].push_back(u);
du[u]++;du[v]++;
}
int T=1000;
while(T--)
{
init();
int rt=rand()%n+1;
memset(vis,0,sizeof(vis));
Build(rt);
dfs(rt,0);
int u=bfs(rt),v=bfs(u);
jump(u,v);
}
printf("%d\n",(int)(ans.size()));
for(int i=0;i<(int)(ans.size());++i)printf("%d\n",ans[i]);
return 0;
}
提供一下我的答案
戳我下載
【NOI2006】聰明的導遊