这是我第一次参加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;
}
文章浏览阅读1w次。开源埋点工具分析查找和比较_开源埋点系统
文章浏览阅读63次。答:选择To check syntax, find problems, and enforce code style(检查语法、发现问题并强制执行代码风格)1.在vue3中,规则vue/name-property--casing已经被vue/component-definition-name-casing代替,在使用前者会出错。a).安装postcss和postcss-preset-env插件:npm install postcss postcss-preset-env -D。(你想遵循哪一种风格指南?
文章浏览阅读1.4w次,点赞6次,收藏11次。使用IDM下载百度云盘文件提示“此站点/路径要求授权”解决方法问题描述解决方法问题描述利用IDM下载百度云盘文件时,未下载完成时提示“此站点/路径要求授权”,要求输入账户和密码,然而输入啥都没用。解决方法停止下载,右键下载项目——属性,可以发现下载的URL改变了。所以要把URL改回来。在浏览器再次下载文件,IDM会弹出“下载文件信息”将图中的URL信息放到之前下载项目右键——属性的..._idm此站点路径要求授权
文章浏览阅读595次,点赞7次,收藏6次。_idea将todo添加到侧边栏
文章浏览阅读208次。【IT168报道】为适应网络强国和数字中国的发展战略,奠定互联网+发展的网络基础,满足垂直行业对网络的多样化需求,迎接万物互联的网络新时代,传统网络向开放化、软件化、智能化的新一代信息网络演进成为产业界的广泛共识。国内三大运..._sdn2018中国
文章浏览阅读1.1k次,点赞18次,收藏26次。校园点餐系统小程序采用B/S模式、采用JAVA语言、springboot框架、mysql数据库、小程序框架uniapp等开工具,促进了校园点餐系统小程序的业务发展。与传统线下点餐相比,校园点餐系统小程序不但节省了卖家的线下运营的成本和门店曝光等难题,对学生而言提供了足不出户送餐上门的服务,针对现代校园而言,出门用餐路程远,排队麻烦等问题,特开发校园点餐系统小程序。本系统使用动态网页开发技术JSP,Java作为系统的开发语言,采用springboot框架、mysql作为后台数据库。_微信点餐小程序源代码
文章浏览阅读1k次。前言接下来就是数据结构的第二大板块:树之前的第一板块:线性表中的顺序表和链表已经更完。但由于队列和栈是特殊的线性表,既然特殊,实际上就普通的顺序表和链表也能实现,加上队列和栈可直接STL实现。故由此,先暂时跳过队列和栈的实现(实际上队列和栈的实现也是仿STL的)笔记(1)为顺序表仿STL的实现,笔记(2)为链表的实现,笔记(3)为队列的仿STL实现,笔记(4)为栈的仿STL实现,而本篇笔记(5)则介绍树形结构树的定义树是一种数据结构,它是由n(n≥0)个有限节点组成一个具有层次关系_三叉树
文章浏览阅读54次。struts简介Struts是Apache软件基金会(ASF)赞助的一个开源项目。它最初是jakarta项目中的一个子项目,并在2004年3月成为ASF的顶级项目。它通过采用JavaServlet/JSP技术,实现了基于JavaEEWeb应用的MVC设计模式的应用框架,是MVC经典设计模式中的一个经典产品。Struts发展历史Struts是作为ApacheJakarta项目的组成部分,...
文章浏览阅读391次。Javaweb分页技术实现分页技术就是通过SQL语句(如下)来获取数据,具体实现看下面代码//分页查询语句select * from 表名 where limit page , count;和//获取表中的总数据,确定页数select count(*) from 表名;不说废话直接上代码前端代码:<!DOCTYPE html>..._javaweb之分页查询 post请求写
文章浏览阅读2.7k次。稳定性范畴, 参考5.x kernel。kernel OopsOops指的就是内核的不正确行为,比如对驱动来说:static int i82092aa_pci_probe(struct pci_dev *dev,const struct pci_device_id *id){unsigned char configbyte;int i, ret;ret = pci_enable_device(de..._kernel warn
文章浏览阅读266次。_forward 和 redirect 区别?
文章浏览阅读5.6k次,点赞3次,收藏7次。前几天我发了个cad三维线模导入Tekla快速建模的文章,有不少朋友让我做个具体的操作视频教程,说实在的,我这把年纪你让我录个视频教程,还要语音的,那真是要我老命了,普通话都讲不好,想来想去还是做个图文教程吧。言归正传首先要快速导入之前,需要下载两个小插件(都是完全免费的)1、SSBIM for CAD下载地址及安装方法:http://zt.ggditu.com/Proscenium/w..._cad三维建模导图插件