标签：轮廓 int 存档 long ch 记忆 du define
很好的一道期望题，暑假时学长讲的，不过当时用的是 WerKeyTom_FTD 给的第二种解法，现在用的是 skyh 学长的决策单调性（分治写法）
首先不难得到期望转移式
我们令 \(a_{i,j}\) 表示存档点在 \(i\) 点 走到 \(j\) 点的期望步数
\(g_{i}\) 表示从该错误节点走到存档点的期望步数
则:
\[a_{i,j+1}=a_{i,j}+\dfrac{1}{du_j}+\left( \dfrac{du_j-1}{du_j}+\dfrac{1}{du_j}\sum_{k\epsilon son(j)}g(k) \right) \]\[a_{i,J+1}=du_j*a_{i,j}+1+S(j) \]其中 \(S(i)=\sum_{k\epsilon son(i)}g(k)\)
接着我们会得到一个类似背包的一个东西：
令 \(f_{i,j}\) 表示在 \(i\) 点存档 共存档 \(j\) 次的最小到达期望步数
则 \(f_{i,j}=min\{ f_{k,j-1}+a_{k,i} \}\)
这样转移是 \(n^3\) 的，我们不妨来观察一下该 \(a\) 式子
\[a_{i,j+1}-a_{i,j} \]\[(du_j-1)*a_{i,j}+du_j+S(j) \]\[a_{i+1,j+1}-a_{i+1,j} \]\[(du_j-1)*a_{i+1,j}+du_j+S(j) \]显然有 \(a_{i,j}>=a_{i+1,j}\)
所以:
\[a_{i,j+1}-a_{i,j}>=a_{i+1,j+1}-a_{i+1,j} \]\[a_{i,j+1}-a_{i+1,j+1}>=a_{i,j}-a_{i+1,j} \]也就是说对于任意 \(l>j\)，\(a_{x,l}\) 的增量大于 \(a_{x,j}\) 的增量
形式化的，我们令函数 \(S(j)=f_{x,p-1}+a_{x,j},S(l)=f_{x,p-1}+a_{x,l}\)
那么两个函数有且仅有一个交点，所以我们其实就是要维护一个类似下凸包的东西
其实也就是对于 \(l\) 点和 \(j\) 点 \((l>j)\) ，\(l\) 的最优决策点一定大于等于 \(j\) 的最优决策点
我们可以用分治来解决此问题
Code
/* memory */
#include <bits/stdc++.h>
#define re register
// #define int long long
#define lls long long
#define fr first
#define sc second
#define pb push_back
#define db double
using namespace std;
const int maxn=1e5+10;
const int mol=998244353;
const db INF=1e99;
inline int read() {
int w=1,s=0; char ch=getchar();
while(ch<'0'||ch>'9') { if(ch=='-') w=-1; ch=getchar(); }
while(ch>='0'&&ch<='9') { s=s*10+ch-'0'; ch=getchar(); }
return s*w;
}
int n,m,t,p,du[2020];
struct EDGE { int var,nxt; } edge[maxn<<1];
int head[maxn],cnt;
inline void add(int a,int b) { edge[++cnt]=(EDGE){ b,head[a] }; head[a]=cnt; }
db g[2020],a[2020][2020],f[2020][2020];
inline void clear() {
for(re int i=1;i<=m;i++) head[i]=-1,g[i]=0,du[i]=0; cnt=0;
}
inline db dfs(int now,int fa) {
bool o=0; db as=0;
for(re int i=head[now],to;~i;i=edge[i].nxt) if((to=edge[i].var)!=fa) {
as+=(dfs(to,now)+1.0)/(db)(du[now]-1.0); o=1;
}
if(!o) return 1.0; return as;
}
inline void wor(int p,int l,int r,int ll,int rr) {
if(l>r) return ;
int fg=ll,mid=(l+r)>>1; db tmp;
f[mid][p]=INF;
for(re int i=ll,lim=min(mid-1,rr);i<=lim;i++) {
tmp=f[i][p-1]+a[i][mid];
if(tmp<f[mid][p]) { f[mid][p]=tmp; fg=i; }
}
wor(p,l,mid-1,ll,fg); wor(p,mid+1,r,fg,rr);
}
inline void slove() {
n=read(); m=read(); p=read(); clear();
for(re int i=1,a,b;i<=m-n;i++) { a=read(); b=read(); add(a,b); add(b,a); du[a]++; du[b]++; }
for(re int now=1;now<=n;now++) { for(re int i=head[now];~i;i=edge[i].nxt) g[now]+=dfs(edge[i].var,now); }
for(re int i=1;i<=n;i++) for(re int j=i+1;j<=n;j++) a[i][j]=a[i][j-1]*(db)(du[j-1]+1.0)+(db)(du[j-1]+1.0)+g[j-1];
f[1][1]=0; for(re int i=2;i<=n;i++) f[i][1]=INF;
for(re int i=2;i<=p;i++) wor(i,1,n,1,n);
printf("%.4lf\n",f[n][p]);
}
signed main(void) {
t=read(); while(t--) { slove(); }
}
标签：轮廓,int,存档,long,ch,记忆,du,define 来源： https://www.cnblogs.com/zjxlm/p/15575378.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。