[JZOJ4648] 【NOIP2016提高A组模拟7.17】锦标赛
题目
描述
题目大意
有nnn个人,你要确定一个出场序列。每次新上台的人就会和擂主打一架,胜利的人继续当擂主。题目给出两两之间打架胜利(失败)的概率。
问111选手坚持到最后的最大概率。
思考历程
看这数据范围这么小,立即想到状压DP!
自然而然地想到状态:设fs,if_{s,i}fs,i表示上了台的人的状态为sss,当前擂主为iii的最大概率。
于是我很快就发现了这个方法的bug。每次转移都会有输或硬两种状态,但万一它们以后要合并到一起呢?
假设状态AAA转移到状态BBB和CCC,接下来状态BBB和CCC又要转移到DDD。按理来说,从AAA到DDD的这两个不同状态是可以合并的,但怎么合并呢?
没有办法合并啊!它们会被看做是不同的状态,然后取个maxmaxmax。
所以我就迷茫了,不知道该怎么做。
开始往什么贪心。结论题方面想……
正解
实际上这题的正解就是状压DP。
我很震惊,难道我找出来的这个东西对答案没有任何影响?
哼,一定是数据错了
然后我就发现了区别:反着和正着的区别。
我一开始并不理解。反着推和正着推有什么不一样吗?状态之间的转移都是一个有向无环图啊!
后来看看题解才发现了本质上的区别:正着推时,每个状态都会被分裂为许多个状态,有些状态在后面还需要合并,有些就不用,这样就很难处理;反着推时,由于我们已经预知了未来,我们只需要从后往前合并状态,并没有什么分裂状态的鬼玩意儿。这样我们就不用想之前的那些情况了。
然后还有另一种神方法,由LYL发明的O(n2)O(n^2)O(n2)贪心做法。
我们思考一下111号选手放在哪里最优。
111号选手要和排在它后面的选手决斗一次,并且必须要胜利,所以概率为这些选手的胜率之积。
胜率总是不超过111的,所以越乘越小。所以将111放在最后面是最优的。
接着我们就用和DP做法一样的思路,从后往前推。设个fif_ifi表示这一个为iii时的最大概率,方程式一样的。
我们贪心地选择其中最大的那个。
接下来就将这个当做新的111号,重复这样的过程。
%%%LYL
代码
状压DP
using namespace std;
#include <cstdio>
#include <cstring>
#include <algorithm>
inline void update(double &a,double b){a<b?a=b:0;}
int n;
double a[18][18];
double f[1<<18][18];
int main(){
scanf("%d",&n);
for (int i=0;i<n;++i)
for (int j=0;j<n;++j)
scanf("%lf",&a[i][j]);
f[(1<<n)-1][0]=1;
for (int i=(1<<n)-2;i>=0;--i)
for (int j=0;j<n;++j)
if (i>>j&1)
for (int k=0;k<n;++k)
if (!(i>>k&1))
update(f[i][j],f[i|1<<k][j]*a[j][k]+f[i|1<<k][k]*a[k][j]);
double ans=0;
for (int i=0;i<n;++i)
update(ans,f[1<<i][i]);
printf("%.7lf\n",ans);
return 0;
}
LYL神贪心
using namespace std;
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 5010
int n;
double a[N][N];
int q[N];
bool vis[N];
double f[N];
int main(){
scanf("%d",&n);
for (int i=1;i<=n;++i)
for (int j=1;j<=n;++j)
scanf("%lf",&a[i][j]);
vis[q[0]=1]=1;
f[0]=-1,f[1]=1;
for (int i=1;i<n;++i){
int mx=0;
for (int j=1;j<=n;++j)
if (!vis[j]){
f[j]=f[j]*a[j][q[i-1]]+f[q[i-1]]*a[q[i-1]][j];
if (f[j]>f[mx])
mx=j;
}
vis[q[i]=mx]=1;
}
printf("%.7lf\n",f[q[n-1]]);
return 0;
}
总结
有时候正着推和反着推真的很不一样……
[JZOJ4648] 【NOIP2016提高A组模拟7.17】锦标赛的更多相关文章
- 【JZOJ4787】【NOIP2016提高A组模拟9.17】数格子
题目描述 输入 输出 样例输入 1 10000 3 10000 5 10000 0 0 样例输出 1 11 95 数据范围 每个测试点数据组数不超过10组 解法 状态压缩动态规划. 设f[i][j]表 ...
- 【NOIP2016提高A组模拟9.17】序列
题目 分析 首先用\(a_i\)表示达到目标的步数\(B_i-A_i(mod 4)\) 根据粉刷栅栏,先不管mod 4的情况,答案就是\(\sum\max(a_i-a_{i+1},0)\) 那我们刚才 ...
- NOIP2016提高A组模拟9.17总结
第一题,典型的隔板问题, 但是我忘记隔板问题怎么打,一开始在花了1小时,还是没想出来,果断弃疗, 最后的40分钟,我打完了第二题,接着又用了20分钟推敲出一种极其猥琐的式子来代替,可惜预处理的阶乘忘记 ...
- 【NOIP2016提高A组模拟9.17】数格子
题目 分析 设表示每一行的状态,用一个4位的二进制来表示,当前这一行中的每一个位数对下一位有没有影响. 设\(f_{i,s}\)表示,做完了的i行,其状态为s,的方案数. 两个状态之间是否可以转移就留 ...
- 【NOIP2016提高A组模拟9.17】小a的强迫症
题目 分析 题目要求第i种颜色的最后一个珠子要在第i+1种颜色的最后一个珠子之前, 那么我们从小到大枚举做到第i种,把第i种的最后一颗珠子取出,将剩下的\(num(i)-1\)个珠子插入已排好的前i- ...
- 【NOIP2016提高A组模拟8.17】(雅礼联考day1)总结
考的还ok,暴力分很多,但有点意外的错误. 第一题找规律的题目,推了好久.100分 第二题dp,没想到. 第三题树状数组.比赛上打了个分段,准备拿60分,因为时间不够,没有对拍,其中有分段的20分莫名 ...
- 【NOIP2016提高A组模拟8.17】(雅礼联考day1)Binary
题目 分析 首先每个数对\(2^i\)取模.也就是把每个数的第i位以后删去. 把它们放进树状数组里面. 那么当查询操作, 答案就位于区间\([2^i-x,2^{i-1}-1-x]\)中,直接查询就可以 ...
- 【NOIP2016提高A组模拟8.17】(雅礼联考day1)Value
题目 分析 易证,最优的答案一定是按\(w_i\)从小到大放. 我们考虑dp, 先将w从小到大排个序,再设\(f_{i,j}\)表示当前做到第i个物品,已选择了j个物品的最大值.转移就是\[f_{i, ...
- 【NOIP2016提高A组模拟8.17】(雅礼联考day1)Matrix
题目 分析 假设,我们从\(F_{i,2}\)出发,那么对\(F_{n,n}\)的贡献就是\(某个系数乘以a^{n-i}b^{n-1}r_i\): 同理,如果从\(F_{2,i}\)出发,那么对\(F ...
随机推荐
- 《DSP using MATLAB》Problem 8.45
代码: %% ------------------------------------------------------------------------ %% Output Info about ...
- sklearn算法中的顶层设计
sklearn监督学习的各个模块 neighbors近邻算法,svm支持向量机,kernal_ridge核岭回归,discriminant_analysis判别分析,linear_model广义线性模 ...
- Python匹马行天下之python之父
龟叔和他的python 经过了漫长的旅程,终于要看到主角Python了.Python是现在非常非常流行的编程语言,在我们能看到的大部分编程语言排行榜中,Python都能在前三甲中拥有一席之地 ,并且发 ...
- UEditor 编辑模板
读取模板,放到ueditor中进行编辑 @model WeiXin_Shop.Models.WX_GoodsDetails @Html.Partial("_MasterPage") ...
- PHP算法之字符串转换整数 (atoi)
请你来实现一个 atoi 函数,使其能将字符串转换成整数. 首先,该函数会根据需要丢弃无用的开头空格字符,直到寻找到第一个非空格的字符为止. 当我们寻找到的第一个非空字符为正或者负号时,则将该符号与之 ...
- CF528E Triangles3000
题意:给你一个不存在三线共交点的一次函数组a[i]x+b[i]y+c[i]=0. 问等概率选取三条直线,围成三角形的面积的期望. n<=3000. 标程: #include<bits/st ...
- Kmeans算法实现
下面的demo是根据kmeans算法原理实现的demo,使用到的数据是kmeans.txt 1.658985 4.285136 -3.453687 3.424321 4.838138 -1.15153 ...
- 11 Python Libraries You Might Not Know
11 Python Libraries You Might Not Know by Greg | January 20, 2015 There are tons of Python packages ...
- django2 rest_framework + vue.js + mysql5.6 实现增删改查
1.安装pymysql,mysqlclient,创建项目django-admin startproject django3 2.在Mysql中创建一个数据库叫django3db,打开项目,修改一下数据 ...
- 多线程--GIL锁
GIL 即全局解释器锁,是一个互斥锁,防止多个线程在同一时间执行python代码,因为在一个python进程中,不仅有主线程而且还有该主线程开启的子线程,还有解释器开启的垃圾回收机等解释器级别的线程. ...