NOIP2017提高組DAY2T2
阿新 • • 發佈:2018-12-19
Chating
說實話,NOIP的暴力分都挺好拿的
比如這道題,前40%的資料,隨便亂搞都可以拿到手(畢竟如此之小)
然後思考一下正解,感覺可以狀壓,但是並不知道怎麼搞。。
後來手賤翻了題解,暴擊++++
很很(hao)簡(li)單(jie)的狀壓啊
打算練習一下狀壓(但是ldw老師說這個dp嘛,練了提升也不大。。我有點方)
Analysis
可以列舉起點i,設立陣列dp[S]表示以i為起點且到達狀態S所需要的最小費用,因此答案就是
然後從起點開始搜尋,每一次列舉已經挖掘的所有點(設當前點為j,當前狀態為S),再列舉這一次要打通的點k。記錄dep[j]表示從j到起點經過的點的數量,也就是題目中的K。則
這樣的複雜度並不是N!,因為dp陣列是有條件更新的(有點類似記搜) (但仍然不知道複雜度怎麼算)
Code(100pts)
#include<bits/stdc++.h>
#define in read()
#define N 30
#define M 2009
#define ll long long
using namespace std;
inline int read(){
int data=0;int w=1; char ch=0;
ch=getchar();
while(ch!='-' && (ch<'0' || ch>'9')) ch=getchar();
if(ch=='-') w=-1,ch=getchar();
while(ch>='0' && ch<='9') data=(data<<3)+(data<<1)+ch-'0',ch=getchar();
return data*w;
}
int n,m;
int f[20][20],inf;
ll dp[ 5000],dep[20],ans=(1ll<<62)-1;
inline void dfs(int sta){
for(int i=1;i<=n;++i){//列舉在狀態S中的點
if(!(sta&(1<<i-1))) continue;
for(int j=1;j<=n;++j){//列舉下一個要更新的點
if(sta&(1<<j-1)) continue;
if(f[i][j]==inf) continue;
if(dp[sta]+dep[i]*1ll*f[i][j]<dp[sta|(1<<j-1)]){
dp[sta|(1<<j-1)]=dp[sta]+dep[i]*1ll*f[i][j];
dep[j]=dep[i]+1;
dfs(sta|(1<<j-1));
}
}
}
}
int main(){
n=in;m=in;
memset(f,127/3,sizeof(f));
inf=f[0][0];
int i,j,u,v,w;
for(i=1;i<=m;++i){
u=in;v=in;w=in;
if(w<f[u][v]) f[u][v]=f[v][u]=w;
}
for(i=1;i<=n;++i){
memset(dep,0,sizeof(dep));
memset(dp,127/3,sizeof(dp));
dep[i]=1;
dp[(1<<i-1)]=0;
dfs(1<<i-1);
ans=min(ans,dp[(1<<n)-1]);
}
cout<<ans;
return 0;
}
Code(40pts)
#include<bits/stdc++.h>
#define in read()
#define N 30
#define M 2009
#define ll long long
using namespace std;
inline int read(){
int data=0;int w=1; char ch=0;
ch=getchar();
while(ch!='-' && (ch<'0' || ch>'9')) ch=getchar();
if(ch=='-') w=-1,ch=getchar();
while(ch>='0' && ch<='9') data=(data<<3)+(data<<1)+ch-'0',ch=getchar();
return data*w;
}
int n,m;
int nxt[M],to[M],head[N],w[M],ecnt=0;
int f[20][20],inf;
inline void add(int x,int y,int z){
nxt[++ecnt]=head[x];head[x]=ecnt;to[ecnt]=y;w[ecnt]=z;
}
int main(){
n=in;m=in;
memset(f,127/3,sizeof(f));//果然啊。。不能開太大,相加容易爆int
inf=f[0][0];
int i,j,k,u,v,w,flag=0,ww;
for(i=1;i<=m;++i){
u=in;v=in;w=in;
add(u,v,w);add(v,u,w);
f[u][v]=f[v][u]=1;
if(i!=1&&w!=ww) flag=1;
ww=w;
}
for(i=1;i<=n;++i) f[i][i]=0;
if(!flag){
for(k=1;k<=n;++k)
for(i=1;i<=n;++i)
for(j=1;j<=n;++j)
if(f[i][k]+f[k][j]<f[i][j]) f[i][j]=f[i][k]+f[k][j];
ll ans=(1ll<<62)-1;
for(i=1;i<=n;++i){//列舉每一個點,作為贊助商幫忙開發的寶藏屋
ll tmp=0;
for(j=1;j<=n;++j){
if(i==j) continue;
if(f[i][j]==inf) break;
tmp+=f[i][j]*1ll*w;
}
if(j<=n) continue;
ans=min(tmp,ans);
}
cout<<ans;
}
return 0;
}
後記
NOIP Day1,2 T1必須要拿下 T2爭取拿下 T3暴力拿完