zoj_2278-Fight for Food

本文介绍了一种通过动态规划寻找捕获最多老鼠的最优路径算法。该算法利用一维数组跟踪每个出现老鼠的位置,并考虑时间因素,确保路径可行。通过预处理可达点之间的最短距离并按时间降序排序点,实现高效求解。
//*************************************************************************************
//1.刚开始容易想到的DP状态是三维的即f[行][列][时间]。这样也可以求解,不过由于时间T可能会很大这样显然会超时。实际只要依次处理每一个会出
//现rat的点平且判断点与点之间是否能在要求的时间内到达就可以了,这样上只需要一维就可以了。
// 
//2.很容易想到按时间顺序对所有的点进行排序并且对于起点L不能到达的点要全部去掉。f[i]=max(f[j])+1 0<=j<i。f[i]表示在第i个点能够catch的最
//大数目的rat
//
//3.还有个重要的优化点,就是点与点之间的时间间隔有可能非常大。例如有的点rat在T=20s出现,有的点在T=1000s出现。其实可以对所有点按rat出现时
//间进行降序排列,并记录下从当前i点(位置[i,j])到所有可到达点的最大距离maxDistance[i][j]。如果rat在点j出现的时间减去rat在点i出现的时间
//大于maxDistance[i][j],那么显然小于j(0<=k<j)的所有点的rat从i点出发必然也都可以catch到
//*************************************************************************************

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

#define max(a,b) (((a) > (b)) ? (a) : (b))
#define N 10
#define NUM 30000

struct Rat{
	int row,col;
	int appear_time;
}rats[NUM];


bool bFlag[N][N];
int maxDistance[N][N];
int f[NUM];
int maxRatsNum[NUM];
int rowNum,colNum;
int ratsNum;

struct Position{
	int row,col;
	int distance;
};
int shorestDistance[N][N][N][N];
Position queue[N*N];

void getShorestDistance(int start_x,int start_y)
{
	if(shorestDistance[start_x][start_y][start_x][start_y]!=-1) return;

	int dir[][2]={{-1,0},{0,1},{1,0},{0,-1}};
	
	Position cur,next;
	int head=0,rear=0;
	cur.row=start_x,cur.col=start_y,cur.distance=0;
	queue[rear++] = cur;
	shorestDistance[start_x][start_y][start_x][start_y] = 0;

	while(head<rear)
	{
		cur = queue[head++];

		for(int i = 0;i<4;++i)
		{
			int newRow = cur.row+dir[i][0],newCol = cur.col+dir[i][1];
			if(newRow>=0&&newRow<rowNum&&newCol>=0&&newCol<colNum&&bFlag[newRow][newCol]&&shorestDistance[start_x][start_y][newRow][newCol]==-1)
			{
				next.row=newRow;
				next.col=newCol;
				next.distance = cur.distance+1;
				queue[rear++]=next;

				shorestDistance[start_x][start_y][newRow][newCol] = next.distance;
				maxDistance[start_x][start_y] = next.distance;
			}
		}
	}
}


int cmpRat(const void* a,const void* b)
{
	return (((Rat*)b)->appear_time-((Rat*)a)->appear_time);
}


int calc_max_rats()
{
	qsort(rats,ratsNum,sizeof(Rat),cmpRat);
	int row,col,newRow,newCol,leftTime;

	for(int i=0;i<ratsNum;++i)
	{
		row = rats[i].row,col=rats[i].col;
		f[i]=1;
		getShorestDistance(row,col);

		for(int j=i-1;j>=0;--j)
		{
			leftTime = rats[j].appear_time - rats[i].appear_time;
			if(leftTime>=maxDistance[row][col])
			{
				if(f[i]<maxRatsNum[j]+1) f[i] = maxRatsNum[j]+1;
				break;
			} 
			else
			{
				newRow = rats[j].row,newCol = rats[j].col;
				if(leftTime>=shorestDistance[row][col][newRow][newCol]&&f[i]<f[j]+1) f[i] = f[j]+1;
			}
		}

		maxRatsNum[i]=(i>0?max(maxRatsNum[i-1],f[i]):f[0]);
	}

	return ratsNum>0?maxRatsNum[ratsNum-1]:0;
}


int main()
{
	//freopen("input.txt","r",stdin);
	char ch;
	while(scanf("%d %d",&rowNum,&colNum)==2)
	{
		memset(bFlag,false,sizeof(bFlag));
		memset(shorestDistance,255,sizeof(shorestDistance));

		int start_x,start_y;
		for(int i = 0;i<rowNum;++i)
		{
			getchar();
			for(int j = 0;j<colNum;++j)
			{
				scanf("%c",&ch);
				if(ch=='.') bFlag[i][j] = true;
				else if(ch=='L'){
					start_x=i;
					start_y=j;
					bFlag[i][j] = true;
				}
			}
		}

		getShorestDistance(start_x,start_y);

		scanf("%d",&ratsNum);
		int row,col,appear_time;
		int i=0,distance;

		while(ratsNum>0)
		{
			scanf("%d %d %d",&row,&col,&appear_time);
			distance = shorestDistance[start_x][start_y][--row][--col];
			
			if(distance!=-1&&appear_time>=distance){
				rats[i].row = row;
				rats[i].col = col;
				rats[i].appear_time = appear_time;
				++i;
			}
			--ratsNum;
		}

		ratsNum = i;
		printf("%d\n",calc_max_rats());
	}

	return 0;
}

一、 内容概要 本资源提供了一个完整的“金属板材压弯成型”非线性仿真案例,基于ABAQUS/Explicit或Standard求解器完成。案例精确模拟了模具(凸模、凹模)与金属板材之间的接触、压合过程,直至板材发生塑性弯曲成型。 模型特点:包含完整的模具-工件装配体,定义了刚体约束、通用接触(或面面接触)及摩擦系数。 材料定义:金属板材采用弹塑性材料模型,定义了完整的屈服强度、塑性应变等真实应力-应变数据。 关键结果:提供了成型过程中的板材应力(Mises应力)、塑性应变(PE)、厚度变化​ 云图,以及模具受力(接触力)曲线,完整再现了压弯工艺的力学状态。 二、 适用人群 CAE工程师/工艺工程师:从事钣金冲压、模具设计、金属成型工艺分析与优化的专业人员。 高校师生:学习ABAQUS非线性分析、金属塑性成形理论,或从事相关课题研究的硕士/博士生。 结构设计工程师:需要评估钣金件可制造性(DFM)或预测成型回弹的设计人员。 三、 使用场景及目标 学习目标: 掌握在ABAQUS中设置金属塑性成形仿真的全流程,包括材料定义、复杂接触设置、边界条件与载荷步。 学习如何调试和分析大变形、非线性接触问题的收敛性技巧。 理解如何通过仿真预测成型缺陷(如减薄、破裂、回弹),并与理论或实验进行对比验证。 应用价值:本案例的建模方法与分析思路可直接应用于汽车覆盖件、电器外壳、结构件等钣金产品的冲压工艺开发与模具设计优化,减少试模成本。 四、 其他说明 资源包内包含参数化的INP文件、CAE模型文件、材料数据参考及一份简要的操作要点说明文档。INP文件便于用户直接修改关键参数(如压边力、摩擦系数、行程)进行自主研究。 建议使用ABAQUS 2022或更高版本打开。显式动力学分析(如用Explicit)对计算资源有一定要求。 本案例为教学与工程参考目的提供,用户可基于此框架进行拓展,应用于V型弯曲
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值