这是我第一次参加CSP，一共得了260分，100,70,70,20,0。这两天试着写一下题解，大家哪里看不懂直接留言问我就好。

第一题：数组推导（100分）

当时不是这么写的，当时建了三个数组，再写的时候发现没必要，记住上一个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;
}

第二题：非零段划分（100分）

下面这段这是满分代码，花了我好长时间，基本上是重写了一遍。减小时间复杂度的方法是通过上一次的非零段划分情况来推算这次的情况。
最近感觉优先队列是个挺好用的东西，不过这题的重点不在优先队列。重点是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;
}

第三题脉冲神经网络（100分）

这个是满分代码，在自己改的基础上参照了其他博主的做法，对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;
}

第四题收集卡牌（100分）

参照另一博主（隔壁李叟）的做法写的，传送门
这其实是一道相对来说比较简单的的状压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;
}

第五题箱根山岳险天下（50分）

这题依然是看的其他人的题解，主要是参考的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,rusing 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;} elser = 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;
}

第23次CSP认证题解相关推荐

第23次CSP认证(202109)
第23次CSP认证---202109 第一题 :数组推导 (贪心) 第二题 :非零段划分 (差分) 第三题 :脉冲神经网络(模拟) 第四题 :收集卡牌( DP ) 第一题 :数组推导 ...
CSP认证题解第一题
文章目录 202012-1 期末预测之安全指数 202009-1 称检测点查询 202006-1 线性分类器 201912-1 报数 201909-1 小明种苹果 201903-1 小中大 20181 ...
2019年3月CCF CSP 认证题解第一题
第一题:201903-1 [题目背景] 在数据分析中,最小值最大值以及中位数是常用的统计信息. [题目描述] 老师给了你n个整数组成的测量数据,保证有序(可能为升序或降序),可能存在重复的数据.请统计 ...
CSP认证题解第二题
文章目录 202012-2 期末预测之最佳阈值 202009-2 风险人群 201909-2 小明种苹果(续) 201903-2 二十四点 201812-2 小明放学 201809-2 买菜 2018 ...
2022年9月CSP认证题解如此编码（k进制），何以包邮？（背包问题），吉祥物投票（珂朵莉树、懒标记、并查集）
T1 如此编码思路由公式和前缀乘积定义得m=b1+a1×b2+⋅⋅⋅+a1×a2×⋅⋅⋅×an−1×bnm=b_1+a_1\times b_2+···+a_1\times a_2\times· ...
CCF CSP认证题解：201703-3 Markdown（Java语言原创）
问题描述 Markdown 是一种很流行的轻量级标记语言(lightweight markup language),广泛用于撰写带格式的文档.例如以下这段文本就是用 Markdown 的语法写成的: ...
第23次CSP认证第4题收集卡牌（记忆化搜索，状压）
链接:http://118.190.20.162/view.page?gpid=T132 思路:先写一个暴力搜索,然后把中间结点的状态存起来,就是记忆化搜索,状态要存搜索到的层数和选择的卡牌数,0代表 ...
CSP认证201509-2 日期计算[C++题解]：枚举、模拟
题目分析来源:acwing 分析: 主要是判断2月多少天.然后就是遍历每个月,依次减去当前月的天数,发现days还剩多少天就是当月几号. 其实这里是背的模板,日期题模板更难一点的题目是 CSP认证 ...
CSP认证 202009-4 星际旅行 Python 题解思路清晰
CSP认证 202009-4 星际旅行 Python 题解思路清晰本题找对数学关系就很容易求了,分三种情况: P1,P2的连线与圆心的距离 >=r,即这条直线在园外,P1和P2的最短距离就 ...
CCF-CSP认证题解
算法设计老师作业要求,看看能补多少补多少吧我去年入围选拔赛坐校车去南海校区第一次参加csp 100+70+70=240 计院第一全校第四,拿了300块奖金那一场是第23次,排名5% 第25次CC ...

第23次CSP认证题解

目录

第一题：数组推导（100分）

第二题：非零段划分（100分）

第三题脉冲神经网络（100分）

第四题收集卡牌（100分）

第五题箱根山岳险天下（50分）

第23次CSP认证题解相关推荐

最新文章

热门文章

第23次CSP认证题解

目录

第一题：数组推导（100分）

第二题：非零段划分（100分）

第三题 脉冲神经网络（100分）

第四题 收集卡牌（100分）

第五题 箱根山岳险天下（50分）

第23次CSP认证题解相关推荐

最新文章

热门文章

第三题脉冲神经网络（100分）

第四题收集卡牌（100分）

第五题箱根山岳险天下（50分）