做题打卡第24天

本文介绍了一种通过树状结构和排序算法解决砂时計問題的方法,该問題涉及砂时計翻转后砂粒数量的变化计算。文章提供了一个具体的编程实现方案。

又到了开心的周三
老规矩
日语练习开始
music~
パーツAとパーツBからなる砂時計があります。これらのパーツにはいくらかの砂が入っています。 砂時計を置くときはパーツAとパーツBのどちらかが上になり、そうでないほうが下になります。

11 秒間に 11 [g] の砂が上にあるパーツから下にあるパーツに落ちます。 ただし、上のパーツにもう砂が残っていない場合は砂は落ちません。

はじめ時刻 00 にパーツAが上にあり、 aa [g] の砂がパーツAに入っていて、 X-aX−a [g] の砂がパーツBに入っています(すなわち、合計 XX [g] の砂が入っています)。

時刻 r_1,r_2,\ …,\ r_Kr
1

,r
2

, …, r
K

に砂時計をひっくり返します。この操作は瞬間的に行われ、時間はかからないものとします。なお、時刻 tt とは時刻 00 の tt 秒後を指します。

クエリが QQ 個与えられます。 各クエリは (t_i,a_i)(t
i

,a
i

) の形をしています。 各クエリに対し、 a=a_ia=a
i

だとして、時刻 t_it
i

にパーツAに入っている砂の量が何gか答えてください。

输入格式
The input is given from Standard Input in the following format:

$ X $
$ K $
$ r_1 $ $ r_2 $ … $ r_K $
$ Q $
$ t_1 $ $ a_1 $
$ t_2 $ $ a_2 $
$ : $
$ t_Q $ $ a_Q $
输出格式
For each query, print the answer in its own line.
原题链接:https://www.luogu.com.cn/problem/AT3527
基本上考的是排序和树
下面是代码
(才不给翻译ne)度娘会告诉你的(:

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;

int n,m,X,tot=0,root,r[100005];
int lch[8000005],rch[8000005];
int minv[8000005],maxv[8000005];
int set[8000005],add[8000005];
struct Question{
	int t,a,id;
}q[100005];
int ans[100005];

bool comp(Question a,Question b){return a.t<b.t;}

int newnode(int &o,int x,int y){
	o=++tot,minv[o]=x,maxv[o]=y,set[o]=-1;
}

void pushup(int o){
	minv[o]=min(minv[lch[o]],minv[rch[o]]);
	maxv[o]=max(maxv[lch[o]],maxv[rch[o]]);
}

void pushdown(int o,int l,int r){
	int mid=l+r>>1;
	int &lc=lch[o],&rc=rch[o];
	if(!lc)newnode(lc,l,mid-1);
	if(!rc)newnode(rc,mid,r-1);
	if(~set[o]){
		minv[lc]=maxv[lc]=set[lc]=set[o];
		minv[rc]=maxv[rc]=set[rc]=set[o];
		add[lc]=add[rc]=0,set[o]=-1;
	}
	if(add[o]){
		minv[lc]+=add[o],maxv[lc]+=add[o],add[lc]+=add[o];
		minv[rc]+=add[o],maxv[rc]+=add[o],add[rc]+=add[o];
		add[o]=0;
	}
}

void Sub(int &o,int l,int r,int x){
	if(o==0)newnode(o,l,r-1);
	if(maxv[o]<=x)
		maxv[o]=minv[o]=add[o]=set[o]=0;
	else if(minv[o]>=x)
		maxv[o]-=x,minv[o]-=x,add[o]-=x;
	else{
		pushdown(o,l,r);
		int mid=l+r>>1;
		Sub(lch[o],l,mid,x),Sub(rch[o],mid,r,x);
		pushup(o);
	}
}

void Add(int &o,int l,int r,int x){
	if(o==0)newnode(o,l,r-1);
	if(minv[o]>=X-x)
		maxv[o]=minv[o]=set[o]=X,add[o]=0;
	else if(maxv[o]<=X-x)
		maxv[o]+=x,minv[o]+=x,add[o]+=x;
	else{
		pushdown(o,l,r);
		int mid=l+r>>1;
		Add(lch[o],l,mid,x),Add(rch[o],mid,r,x);
		pushup(o);
	}
}

int query(int o,int l,int r,int x){
	if(o==0)newnode(o,l,r-1);
	if(r-l==1)return minv[o];
	pushdown(o,l,r);
	int mid=l+r>>1;
	if(x<mid)return query(lch[o],l,mid,x);
	else return query(rch[o],mid,r,x);
}

int main(){
	scanf("%d%d",&X,&n);
	for(int i=1;i<=n;i++)scanf("%d",r+i);
	r[n+1]=1e9;
	scanf("%d",&m);
	for(int i=1;i<=m;i++){
		scanf("%d%d",&q[i].t,&q[i].a);
		q[i].id=i;
	}
	sort(q+1,q+m+1,comp);
	for(int i=0,j=1,rev=0;i<=n;i++,rev^=1){
		while(j<=m&&q[j].t<=r[i+1]){
			int tmp=query(root,1,X+1,q[j].a);
			if(rev)tmp=min(X,tmp+q[j].t-r[i]);
			else tmp=max(0,tmp-q[j].t+r[i]);
			ans[q[j++].id]=tmp;
		}
		if(j>m)break;
		int tmp=r[i+1]-r[i];
		if(rev)Add(root,1,X+1,tmp);
		else Sub(root,1,X+1,tmp);
	}
	for(int i=1;i<=m;i++)printf("%d\n",ans[i]);
} 
一、 内容概要 本资源提供了一个完整的“金属板材压弯成型”非线性仿真案例,基于ABAQUS/Explicit或Standard求解器完成。案例精确模拟了模具(凸模、凹模)与金属板材之间的接触、压合过程,直至板材发生塑性弯曲成型。 模型特点:包含完整的模具-工件装配体,定义了刚体约束、通用接触(或面面接触)及摩擦系数。 材料定义:金属板材采用弹塑性材料模型,定义了完整的屈服强度、塑性应变等真实应力-应变数据。 关键结果:提供了成型过程中的板材应力(Mises应力)、塑性应变(PE)、厚度变化​ 云图,以及模具受力(接触力)曲线,完整再现了压弯工艺的力学状态。 二、 适用人群 CAE工程师/工艺工程师:从事钣金冲压、模具设计、金属成型工艺分析与优化的专业人员。 高校师生:学习ABAQUS非线性分析、金属塑性成形理论,或从事相关课题研究的硕士/博士生。 结构设计工程师:需要评估钣金件可制造性(DFM)或预测成型回弹的设计人员。 三、 使用场景及目标 学习目标: 掌握在ABAQUS中设置金属塑性成形仿真的全流程,包括材料定义、复杂接触设置、边界条件与载荷步。 学习如何调试和分析大变形、非线性接触问题的收敛性技巧。 理解如何通过仿真预测成型缺陷(如减薄、破裂、回弹),并与理论或实验进行对比验证。 应用价值:本案例的建模方法与分析思路可直接应用于汽车覆盖件、电器外壳、结构件等钣金产品的冲压工艺开发与模具设计优化,减少试模成本。 四、 其他说明 资源包内包含参数化的INP文件、CAE模型文件、材料数据参考及一份简要的操作要点说明文档。INP文件便于用户直接修改关键参数(如压边力、摩擦系数、行程)进行自主研究。 建议使用ABAQUS 2022或更高版本打开。显式动力学分析(如用Explicit)对计算资源有一定要求。 本案例为教学与工程参考目的提供,用户可基于此框架进行拓展,应用于V型弯曲
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值