NOIP2017提高组DAY2T2 - 宝藏

传送门

Chating

说实话，NOIP的暴力分都挺好拿的
比如这道题，前40%的数据，随便乱搞都可以拿到手（毕竟$n$如此之小）
然后思考一下正解，感觉可以状压,但是并不知道怎么搞。。
后来手贱翻了题解，暴击++++
很很（hao）简（li）单（jie）的状压啊
打算练习一下状压（但是ldw老师说这个dp嘛，练了提升也不大。。我有点方）

Analysis

可以枚举起点i，设立数组dp[S]表示以i为起点且到达状态S所需要的最小费用，因此答案就是$Min(dp[2^n-1])$

然后从起点开始搜索，每一次枚举已经挖掘的所有点（设当前点为j，当前状态为S），再枚举这一次要打通的点k。记录dep[j]表示从j到起点经过的点的数量，也就是题目中的K。则
$dp[S|k]=Min(dp[S]+dep[j]\ast cost(j\to k))$

这样的复杂度并不是N!，因为dp数组是有条件更新的（有点类似记搜）
（但仍然不知道复杂度怎么算，应该是$O（能过）$吧）

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暴力拿完