【Foreign】Bumb [模拟退火]
Bumb
Time Limit: 20 Sec Memory Limit: 512 MB
Description
Input
Output
Sample Input
4
1 5 1 4
Sample Output
5
HINT
Solution
首先,我们对于一个已知的k,可以O(n)得到Ans,这样就有60%了。
那么怎么做90%呢?老老实实写O(nlogn)是不可能的!模拟退火美滋滋!
传授一点人生的经验吧:写个对拍用于调参,由于这种题不能很单峰,显然温度变化设大一点比较容易正确,无限running + 卡时即可。
是的,没有错!BearChild就这样拿到了90%!
我们来考虑100%怎么做,显然复杂度O(n)。我们把( a[i - 1], a[i] ]看作若干个区间挂在数轴上。
然后考虑k:1->m对于答案的变化。显然可以记录cnt表示k包含的区间个数,以及sum和,这样就可以做完啦!
Code
#include<iostream>
#include<string>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<vector>
#include<ctime>
using namespace std;
typedef long long s64; const int ONE = ;
const s64 INF = 1e18; int n, m;
int a[ONE];
int Now, A;
s64 Ans = INF; int get()
{
int res=,Q=;char c;
while( (c=getchar())< || c> )
if(c=='-')Q=-;
res=c-;
while( (c=getchar())>= && c<= )
res=res*+c-;
return res*Q;
} void TimeE()
{
if((double)clock() / CLOCKS_PER_SEC > 0.97)
{
printf("%lld", Ans);
exit();
}
} s64 Get(int x, int y)
{
TimeE();
if(x <= y) return y - x;
return y + m - x;
} s64 Judge(int k)
{
s64 res = ;
for(int i = ; i <= n; i++)
res += min(Get(a[i - ], a[i]), + Get(k, a[i]));
Ans = min(Ans, res);
return res;
} double Random() {return rand() / (double)RAND_MAX;}
void SA(double T)
{
Now = m / ;
while(T >= )
{
A = Now + (int)(T * (Random() * - ));
if(A < || A > m) A = T * Random();
s64 dE = Judge(A) - Judge(Now);
if(dE < ) Now = A;
T *= 0.92;
}
} int main()
{
n = get(); m = get();
for(int i = ; i <= n; i++)
a[i] = get();
if(n > ) {for(;;)SA(m); return ;} for(int i = ; i <= m; i++)
Ans = min(Ans, Judge(i)); printf("%lld", Ans);
}
模拟退火 90%
#include<iostream>
#include<string>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<vector>
using namespace std;
typedef long long s64; const int ONE = ;
const s64 INF = 1e18; int n, m;
int a[ONE];
int L[ONE];
vector <int> R[ONE];
s64 Ans; int get()
{
int res=,Q=;char c;
while( (c=getchar())< || c> )
if(c=='-')Q=-;
res=c-;
while( (c=getchar())>= && c<= )
res=res*+c-;
return res*Q;
} int Get(int x, int y)
{
if(x <= y) return y - x;
return y + m - x;
} int main()
{
n = get(); m = get();
for(int i = ; i <= n; i++)
a[i] = get(); for(int i = ; i <= n; i++)
{
Ans += Get(a[i - ], a[i]);
L[a[i - ]]++, R[a[i]].push_back(a[i - ]);
} s64 sum = Ans, cnt = ; for(int i = ; i <= n; i++)
if(a[i - ] > a[i])
sum -= m - a[i - ], cnt++; for(int k = ; k <= m; k++)
{
int len = R[k - ].size();
cnt = cnt + L[k - ] - len;
sum -= cnt;
for(int i = ; i < len; i++)
sum += Get(R[k - ][i], k - ) - ;
Ans = min(Ans, sum);
} printf("%lld", Ans);
}
正解 100%
【Foreign】Bumb [模拟退火]的更多相关文章
- 【Foreign】咏叹 [模拟退火]
咏叹 Time Limit: 100 Sec Memory Limit: 256 MB Description 有n根木棍,第i根长度为ai.你要贴着墙围出一个矩形区域,木棍围成的矩形边缘必须平行或 ...
- bzoj3680模拟退火
看题意就是一道数学物理题,带权费马点 --这怎么是数学了,这也是物理的 所以要用物理方法,比如FFF 国际著名oi选手miaom曾说 模拟退火初温可以低,但是最好烧个几千次 国际著名物理课代表+1 ...
- Constraint6:更新外键约束(Foreign Key Constraint)的引用列
在SQL Server中,表之间存在引用关系,引用关系通过创建外键约束(Foreign Key Constraint)实现.如果一个Table中的column被其他Table引用,那么该表是参考表,或 ...
- MySQL主从复制中断,报“Error on master: message (format)='Cannot delete or update a parent row: a foreign key constraint fails' error code=1217” 错误
前几天,发现从库挂了,具体报错信息如下: 分析思路 1. 因为我采用的是选择性复制,只针对以下几个库进行复制: card,upay,deal,monitor,collect.所以,不太可能出现对于sa ...
- ORA-02266: unique/primary keys in table referenced by enabled foreign keys
在数据库里面使用TRUNCATE命令截断一个表的数据时,遇到如下错误 SQL >TRUNCATE TABLE ESCMOWNER.SUBX_ITEM ORA-02266: unique/prim ...
- SQL Server 2008 R2——TRUNCATE TABLE 无法截断表 该表正由 FOREIGN KEY 约束引用
=================================版权声明================================= 版权声明:原创文章 禁止转载 请通过右侧公告中的“联系邮 ...
- 【MySQL】Create table 以及 foreign key 删表顺序考究。
1.以下是直接从数据库导出的建表语句. 1 -- ---------------------------- 2 -- Table structure for files 3 -- ---------- ...
- MYSQL外键(Foreign Key)的使用
在MySQL 3.23.44版本后,InnoDB引擎类型的表支持了外键约束.外键的使用条件:1.两个表必须是InnoDB表,MyISAM表暂时不支持外键(据说以后的版本有可能支持,但至少目前不支持): ...
- 无题的题 & 模拟退火...
题意: 给你不超过8条一端在圆心的半径,求他们组成的凸包的最大面积. SOL: 正解怎么搞啊不会啊...然后昨天毛爷爷刚讲过模拟退火...那么就打一个吧... 然后就T了,不过三角形的部分分妥妥的.. ...
随机推荐
- Android开发随笔5
昨天: 对界面的进一步设计补充 可以在界面之间的跳转 研究了对图标等的操作 今天: 实现对库的相关操作 学习视视频内容‘ 复习java的一些知识.
- mysql 复杂查询
1.同一个表下多次查询: sql语句: select b.* ,(select name from exh_common.medicine_type a where b.p_id = a.id) as ...
- lintcode-172-删除元素
172-删除元素 给定一个数组和一个值,在原地删除与值相同的数字,返回新数组的长度. 元素的顺序可以改变,并且对新的数组不会有影响. 样例 给出一个数组 [0,4,4,0,0,2,4,4],和值 4 ...
- LintCode-211.字符串置换
字符串置换 给定两个字符串,请设计一个方法来判定其中一个字符串是否为另一个字符串的置换. 置换的意思是,通过改变顺序可以使得两个字符串相等. 样例 "abc" 为 "cb ...
- 利用闭包判断Dom元素和滚动条的方向
本文收集整理自网上. 一,判断滚动条的方向,利用闭包首先保存滚动条的位置,然后当滚动时候不断更新滚动初始值,然后通过差指判断方向 function scroll(fn) { //利用闭包判断滚动条滚动 ...
- 微信小程序 功能函数 计时器
let lovetime = setInterval(function () { let str = '(' + n + ')' + '重新获取' that.setData({ getText2: s ...
- react-event-pooling
react-event-pooling 事件池 https://codesandbox.io/s/3xp4y9zp7q https://reactjs.org/docs/events.html#eve ...
- 2019 front end jobs collection
2019 front end jobs collection Alibaba https://ant.design/docs/spec/work-with-us-cn https://www.yuqu ...
- P3986 斐波那契数列
题目描述 定义一个数列: f(0)=a,f(1)=b,f(n)=f(n−1)+f(n−2) 其中 a,b均为正整数,n≥2 . 问有多少种 (a,b),使得 k 出现在这个数列里,且不是前两项. 由于 ...
- [ZJOI2010]数字计数 数位DP
最近在写DP,今天把最近写的都放上来好了,,, 题意:给定两个正整数a和b,求在[a,b]中的所有整数中,每个数码(digit)各出现了多少次. 首先询问的是一个区间,显然是要分别求出1 ~ r ,1 ...