【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了,不过三角形的部分分妥妥的.. ...
随机推荐
- Java容器之Collections
Collections 类来源于 java.util.Collections,从 java.lang.object继承. 此类完全由在 collection 上进行操作或返回 collection 的 ...
- lintcode-191-乘积最大子序列
191-乘积最大子序列 找出一个序列中乘积最大的连续子序列(至少包含一个数). 样例 比如, 序列 [2,3,-2,4] 中乘积最大的子序列为 [2,3] ,其乘积为6. 标签 子数组 领英 动态规划 ...
- OSG数学基础:坐标系变换
三维实体对象需要经过一系列的坐标变换才能正确.真实地显示在屏幕上.在一个场景中,当读者对场景中的物体进行各种变换及相关操作时,坐标系变换是非常频繁的. 坐标系变换通常包括:世界坐标系-物体坐标系变换. ...
- mysql授权远程连接
查一下你的MYSQL用户表里, 是否允许远程连接 1.授权 mysql>grant all privileges on *.* to 'root'@'%' identified by ...
- printf以及各种变种
int printf(char *format, [argument]); 其向终端(显示器.控制台等)输出字符 int fprintf(FILE*stream, const char*format, ...
- C# 使用this的形参
示例1: public static RectangleF TransformRect(this Matrix mat, RectangleF rect) 是向Matrix类扩展带有Rectangle ...
- JVM(二)GC
GC简介 Java堆内存 在运行时,java的实例被存放在堆内存区域.当一个对象不在被引用,满足条件就会从堆内存移除并且内存空间被回收.堆内存由三个主要区域 1.新生代 Eden空间(任何实例 ...
- [STL] vector基本用法
vector的数据安排以及操作方式,与array非常相似.两者的唯一区别在于空间的运用的灵活性.array是静态空间,一旦配置了就不能改变.vector是动态空间,随着元素的加入,它的内部机制会自行扩 ...
- 在select中,载入时默认select为空白,选项内不显示空白项
有两种办法,一种纯css实现,一种借助js实现. html: <body onload="load()"> <select id="abc" ...
- WPF对某控件添加右键属性
代码创建右键属性 ContextMenu cm = new ContextMenu(); MenuItem mi = new MenuItem(); mi.Header = "打开此文件所有 ...