这是我第一次参加CSP,一共得了260分,100,70,70,20,0。这两天试着写一下题解,大家哪里看不懂直接留言问我就好。
当时不是这么写的,当时建了三个数组,再写的时候发现没必要,记住上一个B的值就可以了,这样可以节省空间(虽然没太大必要)
#include <bits/stdc++.h>
#define For(i,n) for(int i=0;i<n;i++)
using namespace std;
int n;
int lastb,b;
int aMax;
int aMin;
int sumMax=0;
int sumMin=0;
int main() {
scanf("%d",&n);
int b;
For(i,n){
scanf("%d",&b);
if(i==0){
aMax = b;
aMin = b;
}
else{
if(b>lastb){
aMax = b;
aMin = b;
}
else if(b==lastb){
aMax = b;
aMin = 0;
}
}
lastb=b;
sumMax+=aMax;
sumMin+=aMin;
}
printf("%d\n%d\n",sumMax,sumMin);
return 0;
}
下面这段这是满分代码,花了我好长时间,基本上是重写了一遍。减小时间复杂度的方法是通过上一次的非零段划分情况来推算这次的情况。
最近感觉优先队列是个挺好用的东西,不过这题的重点不在优先队列。重点是c[i]和lastc[i]这两个数组,它们表示的都是第i个数是否已经变成0。lastc是上一次的c。还有一个比较细节的地方是我从其他博客里看来的。首先,优先队列里存的是元素值与位置的映射,细节就在于多往里存个(0,0)和(0,n+1),这样后续的处理就会更加方便。
大家看不懂优先队列的可以先把优先队列里的元素按顺序打印出来,应该就比较容易看懂了。
这题比较容易出错的点就是ans在什么时候加或减,if里的条件要写全,我一开始样例都对了,但就是爆0,后来自己写了个
5
1 2 1 2 1
的测试样例,才发现自己的if条件还有问题。
#include <bits/stdc++.h>
using namespace std;
int n;
int ans = 0;
int maxans = 0;
typedef pair<int, int> pii;
priority_queue<pii, vector<pii>, greater<pii>> Q;
bool lastc[500005];
bool c[500005];
int main() {
int ai;
scanf("%d", &n);
Q.push(pii(0,0));
for (int i = 1; i <= n; i++) {
scanf("%d", &ai);
Q.push(pii(ai, i));
}
Q.push(pii(0, n+1));
memset(c, false, sizeof c);
memset(lastc, false, sizeof c);
pair<int, int> q = Q.top();
Q.pop();
int left = q.second;
int right = q.second;
while (!Q.empty()) {
while (Q.top().first == q.first && !Q.empty()) {
c[q.second] = true;
if (q.second == Q.top().second - 1)
++right;
else {
if (left > 0 && !c[left - 1] && ((right < n+1 && !lastc[right + 1]) || right==n+1))
++ans;
if (right < n+1 && left > 0 && lastc[right + 1] && lastc[left - 1])
--ans;
left = Q.top().second;
right = Q.top().second;
}
q = Q.top();
Q.pop();
}
c[q.second] = true;
if (left > 0 && !c[left - 1] && ((right < n+1 && !lastc[right + 1]) || right==n+1))
++ans;
if (right < n+1 && left > 0 && lastc[right + 1] && lastc[left - 1])
--ans;
if (maxans < ans)
maxans = ans;
//cout << q.first << " " << ans <<endl;
if (!Q.empty()) {
q = Q.top();
Q.pop();
left = q.second;
right = q.second;
}
memcpy(lastc, c, sizeof c);
}
printf("%d", maxans);
return 0;
}
这个是满分代码,在自己改的基础上参照了其他博主的做法,对Dmax取模来节省空间。
传送门
在这个过程中我遇到一个问题,我按着和其他博主一样的思路写的代码就是超时,后来我把放在神经元结构体里的w数组拿出来就不超时了,
这告诉我们一个道理,数组容量大的时候不要放在结构体里面。。
关于G数组的作用,有人问到了,我解释一下
#include <bits/stdc++.h>
using namespace std;
static unsigned long nex = 1;
/* RAND_MAX assumed to be 32767 */
int myrand() {
nex = nex * 1103515245 + 12345;
return ((unsigned) (nex / 65536) % 32768);
}
int N, S, P, T;
double ct;
double maxv = -DBL_MAX;
double minv = DBL_MAX;
double v, u, a, b, c, d;
int maxtim = -INT_MAX;
int mintim = INT_MAX;
struct SJY {
double v, u, a, b, c, d;
int tim;
} sjy[1005];
double getww[1005][1005];
vector<int> G[2005];
struct TU {
int s, t;
double w;
int D;
};
vector<TU> tu;
int r[1005];
int maxD=0;
int main() {
int sjyn = 0;
int RN;
scanf("%d%d%d%d", &N, &S, &P, &T);
scanf("%lf", &ct);
while (true) {
scanf("%d%lf%lf%lf%lf%lf%lf", &RN, &v, &u, &a, &b, &c, &d);
while (RN--) {
sjy[sjyn].v = v;
sjy[sjyn].u = u;
sjy[sjyn].a = a;
sjy[sjyn].b = b;
sjy[sjyn].c = c;
sjy[sjyn].d = d;
sjy[sjyn].tim = 0;
sjyn++;
}
if (sjyn == N)
break;
}
for (int i = 0; i < P; i++)
scanf("%d", &r[i]);
int s, t, D;
double w;
for (int i = 0; i < S; i++) {
scanf("%d%d%lf%d", &s, &t, &w, &D);
G[s].push_back(i);
tu.push_back({
s, t, w, D});
if(D>maxD)
maxD=D;
}
maxD++;
for (int nowt = 1; nowt <= T; nowt++) {
for (int j = 0; j < N; j++) {
double v1 = sjy[j].v, u1 = sjy[j].u;
sjy[j].v = v1 + ct * (0.04 * v1 * v1 + 5 * v1 + 140 - u1) + getww[nowt%maxD][j];
sjy[j].u = u1 + ct * sjy[j].a * (sjy[j].b * v1 - u1);
getww[nowt%maxD][j] = 0;
if (sjy[j].v >= 30) {
sjy[j].tim++;
sjy[j].v = sjy[j].c;
sjy[j].u += sjy[j].d;
for (int ii = 0; ii < G[j].size(); ii++) {
TU tu1(tu[G[j][ii]]);
getww[(nowt+tu1.D)%maxD][tu1.t] += tu1.w;
}
}
}
for (int j = 0; j < P; j++) {
if (r[j] > myrand()) {
for (int ii = 0; ii < G[j + N].size(); ii++) {
TU tu1(tu[G[j + N][ii]]);
getww[(nowt+tu1.D)%maxD][tu1.t] += tu1.w;
}
}
}
}
for(int j=0;j<N;j++){
if (sjy[j].v > maxv)
maxv = sjy[j].v;
if (sjy[j].v < minv)
minv = sjy[j].v;
if (sjy[j].tim > maxtim)
maxtim = sjy[j].tim;
if (sjy[j].tim < mintim)
mintim = sjy[j].tim;
}
printf("%.3f %.3f\n", minv, maxv);
printf("%d %d", mintim, maxtim);
return 0;
}
参照另一博主(隔壁李叟)的做法写的,传送门
这其实是一道相对来说比较简单的的状压dp,但是自己推dp水平太差了。
#include <bits/stdc++.h>
using namespace std;
int n, k;
double p[17];
double dp[1 << 17][77];
bool d[1 << 17][77];
double dfs(int t, int sum, int val, double l) {
if (val + (sum - val) / k >= n)
return sum;
if (d[t][sum])
return dp[t][sum];
for (int i = 0; i < n; i++) {
if (!((t >> i) & 1)) {
dp[t][sum] += p[i] * dfs(t | (1 << i), sum + 1, val + 1, l + p[i]);
}
}
if (t) {
dp[t][sum] += l * dfs(t, sum + 1, val, l);
}
d[t][sum] = true;
return dp[t][sum];
}
int main() {
scanf("%d%d", &n, &k);
for (int i = 0; i < n; i++) {
scanf("%lf", &p[i]);
}
printf("%.10lf", dfs(0, 0, 0, 0));
return 0;
}
这题依然是看的其他人的题解,主要是参考的CSDN上另一位博主的代码,传送门
用到的知识是树链剖分线段树。
因为没有用LCT(动态树),所以只能拿50分。(异或操作强制在线,树链剖分要知道树的形态)。LCT目前我还没学,我上面发的那个博主写了,b站也有一个题解。大家可以自行搜索。
#include <bits/stdc++.h>
#define mid ((l+r)>>1)
#define lson rt<<1,l,mid
#define rson rt<<1|1,mid+1,r
using namespace std;
typedef long long ll;
const int maxm = 300000 + 10;
struct req {
ll type, x2, s, l, r, y2;
req(ll type, ll x2, ll s, ll l, ll r, ll y2)
: type(type), x2(x2), s(s), l(l), r(r), y2(y2) {
}
};
vector<req> R;
ll m, p, T;
//ll e = 0, beg[maxm], nex[maxm], to[maxm];
ll dep[maxm], siz[maxm], fa[maxm], son[maxm];
ll id[maxm], top[maxm], cnt = 0;
ll sum[maxm << 2], lazy[maxm << 2];
ll now, res, tot = 0,rt;
ll A = 0;
vector<ll> edge[maxm];
vector<pair<ll, ll> > E[maxm];
//inline void add(ll x, ll y) {
// to[++e] = y;
// nex[e] = beg[x];
// beg[x] = e;
//}
inline void dfs1(ll x) {
//x当前节点,f父亲,deep深度
siz[x] = 1;
ll maxson = 0;
for (ll y:edge[x]) {
// ll y = to[i];
if (y == fa[x])
continue;
dfs1(y);
siz[x] += siz[y];
if (siz[y] > maxson) {
son[x] = y;
maxson = siz[y];
}
}
}
inline void dfs2(ll x, ll topf) {
//x当前节点,topf当前链的最顶端的节点
id[x] = ++cnt;
top[x] = topf;
if (!son[x])
return;
dfs2(son[x], topf);
for (ll y:edge[x]) {
// ll y = to[i];
if (y == fa[x] || y == son[x])
continue;
dfs2(y, y);
}
}
inline void pushdown(ll rt) {
if (lazy[rt] != 1) {
lazy[rt << 1] = (lazy[rt] * lazy[rt << 1]) % p;
lazy[rt << 1 | 1] = (lazy[rt] * lazy[rt << 1 | 1]) % p;
sum[rt << 1] = (lazy[rt] * sum[rt << 1]) % p;
sum[rt << 1 | 1] = (lazy[rt] * sum[rt << 1 | 1]) % p;
lazy[rt] = 1;
}
}
inline void insert(ll rt, ll l, ll r, ll v, ll y) {
if (l == r) {
sum[rt] = y;
sum[rt] %= p;
lazy[rt] = 1;
return;
} else {
if (lazy[rt] != 1)
pushdown(rt);
if (v <= mid)
insert(lson, v, y);
if (v > mid)
insert(rson, v, y);
sum[rt] = (sum[rt << 1] + sum[rt << 1 | 1]) % p;
}
}
inline void query(ll rt, ll l, ll r, ll L, ll R) {
if (L <= l && r <= R) {
res += sum[rt];
res %= p;
return;
} else {
if (lazy[rt] != 1)
pushdown(rt);
if (L <= mid)
query(lson, L, R);
if (R > mid)
query(rson, L, R);
}
}
inline void update(ll rt, ll l, ll r, ll L, ll R, ll k) {
if (L <= l && r <= R) {
lazy[rt] = (lazy[rt] * k) % p;
sum[rt] = (sum[rt] * k) % p;
} else {
if (lazy[rt] != 1)pushdown(rt);
if (L <= mid)update(lson, L, R, k);
if (R > mid)update(rson, L, R, k);
sum[rt] = (sum[rt << 1] + sum[rt << 1 | 1]) % p;
}
}
inline ll qRange(ll x, ll y) {
ll ans = 0;
while (top[x] != top[y]) {
if (dep[top[x]] < dep[top[y]])swap(x, y);
res = 0;
query(1, 1, tot, id[top[x]], id[x]);
ans += res;
ans %= p;
x = fa[top[x]];
}
if (dep[x] > dep[y])swap(x, y);
res = 0;
query(1, 1, tot, id[x], id[y]);
ans += res;
return ans % p;
}
inline void updRange(ll x, ll y, ll k) {
k %= p;
while (top[x] != top[y]) {
if (dep[top[x]] < dep[top[y]]) swap(x, y);
update(1, 1, tot, id[top[x]], id[x], k);
x = fa[top[x]];
}
if (dep[x] > dep[y]) swap(x, y);
update(1, 1, tot, id[x], id[y], k);
}
inline ll find(ll x, ll y) {
ll l = 0, r = E[x].size() - 1, ans = 0;
while (l <= r) {
ll mi = (l + r) >> 1;
if (E[x][mi].first <= y) {
ans = E[x][mi].second;
l = mi + 1;
} else
r = mi - 1;
}
return ans;
}
int main() {
ll type;
scanf("%lld%lld%lld", &m, &p, &T);
rt = now = 0;
tot = 0;
for (int i = 1; i <= m << 2; i++) {
lazy[i] = 1;
sum[i] = 0;
}
dep[0] = 0;
for (ll t = 1, x2, y2, s, l, r; t <= m; t++) {
scanf("%lld", &type);
x2 = y2 = s = l = r = 0;
if (type == 1) {
scanf("%lld", &x2);
if (x2 == 0) {
now = fa[now];
} else {
s = ++tot;
fa[tot] = now;
dep[tot] = dep[now] + 1;
edge[now].push_back(tot);
// add(now, tot);
// add(tot, now);
E[dep[tot]].push_back({
t, tot});
now = tot;
}
} else if (type == 2) {
scanf("%lld%lld%lld%lld", &s, &l, &r, &y2);
} else if (type == 3) {
scanf("%lld%lld%lld", &s, &l, &r);
}
R.push_back(req(type, x2, s, l, r, y2));
}
A = 0;
tot++;
dfs1(rt);
dfs2(rt, rt);
for (req re: R) {
if (re.type == 1) {
if (re.x2) {
insert(1, 1, tot, id[re.s], re.x2 ^ A);
}
} else {
ll x = find(re.l, re.s);
ll y = find(re.r, re.s);
if (re.type == 2) {
updRange(x, y, re.y2 ^ A);
} else {
res = qRange(x, y);
printf("%lld\n", res);
if (T)
A = res;
}
}
}
return 0;
}
文章浏览阅读331次。第一部分:准备工作1 安装虚拟机2 安装centos73 安装JDK以上三步是准备工作,至此已经完成一台已安装JDK的主机第二部分:准备3台虚拟机以下所有工作最好都在root权限下操作1 克隆上面已经有一台虚拟机了,现在对master进行克隆,克隆出另外2台子机;1.1 进行克隆21.2 下一步1.3 下一步1.4 下一步1.5 根据子机需要,命名和安装路径1.6 ..._创建一个hadoop项目
文章浏览阅读1.7k次。心脏滴血漏洞HeartBleed CVE-2014-0160 是由heartbeat功能引入的,本文从深入码层面的分析该漏洞产生的原因_heartbleed代码分析
文章浏览阅读1.4k次。前言ofd是国家文档标准,其对标的文档格式是pdf。ofd文档是容器格式文件,ofd其实就是压缩包。将ofd文件后缀改为.zip,解压后可看到文件包含的内容。ofd文件分析工具下载:点我下载。ofd文件解压后,可以看到如下内容: 对于xml文件,可以用文本工具查看。但是对于印章文件(Seal.esl)、签名文件(SignedValue.dat)就无法查看其内容了。本人开发一款ofd内容查看器,..._signedvalue.dat
文章浏览阅读1.8w次,点赞29次,收藏313次。整体系统设计本设计主要是对ADC和DAC的使用,主要实现功能流程为:首先通过串口向FPGA发送控制信号,控制DAC芯片tlv5618进行DA装换,转换的数据存在ROM中,转换开始时读取ROM中数据进行读取转换。其次用按键控制adc128s052进行模数转换100次,模数转换数据存储到FIFO中,再从FIFO中读取数据通过串口输出显示在pc上。其整体系统框图如下:图1:FPGA数据采集系统框图从图中可以看出,该系统主要包括9个模块:串口接收模块、按键消抖模块、按键控制模块、ROM模块、D.._基于fpga的信息采集
文章浏览阅读2.5w次。1.背景错误信息:-- [http-nio-9904-exec-5] o.s.c.n.z.filters.post.SendErrorFilter : Error during filteringcom.netflix.zuul.exception.ZuulException: Forwarding error at org.springframework.cloud..._com.netflix.zuul.exception.zuulexception
文章浏览阅读358次。1.介绍图的相关概念 图是由顶点的有穷非空集和一个描述顶点之间关系-边(或者弧)的集合组成。通常,图中的数据元素被称为顶点,顶点间的关系用边表示,图通常用字母G表示,图的顶点通常用字母V表示,所以图可以定义为: G=(V,E)其中,V(G)是图中顶点的有穷非空集合,E(G)是V(G)中顶点的边的有穷集合1.1 无向图:图中任意两个顶点构成的边是没有方向的1.2 有向图:图中..._给定一个邻接矩阵未必能够造出一个图
文章浏览阅读321次。(十二)、WDS服务器安装通过前面的测试我们会发现,每次安装的时候需要加域光盘映像,这是一个比较麻烦的事情,试想一个上万个的公司,你天天带着一个光盘与光驱去给别人装系统,这将是一个多么痛苦的事情啊,有什么方法可以解决这个问题了?答案是肯定的,下面我们就来简单说一下。WDS服务器,它是Windows自带的一个免费的基于系统本身角色的一个功能,它主要提供一种简单、安全的通过网络快速、远程将Window..._doc server2012上通过wds+mdt无人值守部署win11系统.doc
文章浏览阅读219次。python–xlrd/xlwt/xlutilsxlrd只能读取,不能改,支持 xlsx和xls 格式xlwt只能改,不能读xlwt只能保存为.xls格式xlutils能将xlrd.Book转为xlwt.Workbook,从而得以在现有xls的基础上修改数据,并创建一个新的xls,实现修改xlrd打开文件import xlrdexcel=xlrd.open_workbook('E:/test.xlsx') 返回值为xlrd.book.Book对象,不能修改获取sheett_xlutils模块可以读xlsx吗
文章浏览阅读8.2w次,点赞267次,收藏656次。运行Selenium出现'WebDriver' object has no attribute 'find_element_by_id'或AttributeError: 'WebDriver' object has no attribute 'find_element_by_xpath'等定位元素代码错误,是因为selenium更新到了新的版本,以前的一些语法经过改动。..............._unresolved attribute reference 'find_element_by_id' for class 'webdriver
文章浏览阅读198次。一:模态窗口//父页面JSwindow.showModalDialog(ifrmehref, window, 'dialogWidth:550px;dialogHeight:150px;help:no;resizable:no;status:no');//子页面获取父页面DOM对象//window.showModalDialog的DOM对象var v=parentWin..._jquery获取父window下的dom对象
文章浏览阅读1.7w次,点赞15次,收藏129次。算法(algorithm)是解决一系列问题的清晰指令,也就是,能对一定规范的输入,在有限的时间内获得所要求的输出。 简单来说,算法就是解决一个问题的具体方法和步骤。算法是程序的灵 魂。二、算法的特征1.可行性 算法中执行的任何计算步骤都可以分解为基本可执行的操作步,即每个计算步都可以在有限时间里完成(也称之为有效性) 算法的每一步都要有确切的意义,不能有二义性。例如“增加x的值”,并没有说增加多少,计算机就无法执行明确的运算。 _算法
文章浏览阅读1.5k次,点赞18次,收藏26次。网络安全的标准和规范是网络安全领域的重要组成部分。它们为网络安全提供了技术依据,规定了网络安全的技术要求和操作方式,帮助我们构建安全的网络环境。下面,我们将详细介绍一些主要的网络安全标准和规范,以及它们在实际操作中的应用。_网络安全标准规范