CSP 202309-3 梯度求解 详解（60分debug到100分）

已于 2024-08-18 19:40:16 修改
阅读量525 收藏 3
点赞数 10
分类专栏： CSP 文章标签： 算法 c++
于 2024-08-18 19:31:15 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/2301_79020141/article/details/141303270
版权

题目浅析

        题目主要需要解决，如何将逆波兰式还原为表达式并进行运用，以及如何计算整个表达式的导数，两大问题。

解决逆波兰式
要点

        根据题目提到的逆波兰式联系通过二叉树的后序遍历，通过建二叉树树将原表达式存储。

详解

        定义二叉树结点时，定义id来识别该结点，并用于进行左右孩子的连接。同时，将每个结点对应的字符串存储下来。此外为了便于后续计算导数还定义了d_v表示以该节点为根的子树的表达式的导数的值，pre_v表示以该节点为根的子树的表达式， 直接带入各个数值的原始值。（这两个属性的用途后续会再解释）。

        定义结点类型的栈用于将逆波兰式还原为表达式的过程。

        开始遍历逆波兰式中以空格为分隔的每一个字符串，根据逆波兰式是表达式二叉树后续遍历的结果可得，叶子结点都是自变量或常数，运算符号一定都有左右孩子。

        据上，当遍历到的字符串是“+”或“-”或“*”某运算符号时，首先建立一个树的结点赋值id和字符串po，接着从栈中取出两个结点分别作为右孩子和左孩子，注意此处先出栈的元素作为右孩子，顺序不可交换，否则"-"时会出现错误。将左右孩子设定好之后，将该结点入栈。

        当遍历到其他字符串即自变量或常数时，建立结点，对相应属性赋值，并将结点放入栈中。

        直到最后一个字符串建立好结点，建树完成。

计算导数
要点

        计算以每个结点为根节点的子树对应的表达式的导数和原始值，根节点对应的表达式的导数即为整个表达式的结果。

 详解

        以样例1的表达式x1 *（x1 * x1 + x2）为例进行分析 计算整体的导数最终都会落实到计算局部的导数和局部的表达式的值，对原式求导得到{x1}' *(x1*x1+x2)+x1*{(x1*x1+x2)}'，其他项此时可以得到，继续对{(x1*x1+x2)}'计算，得到{(x1*x1)}'+{x2}'，继续对{(x1*x1)}'部分进行求导，得到{x1}'*x1+x1*{x1}'，此时带入各值即可得到对原表达式求导的值。

        根据上述计算过程可以考虑在代码中实现以下过程：将每个结点处，以该节点为根的子树形成表达式的导数和原始值先计算出来，即，由于计算整体的导数最终都会落实到，计算局部的导数和局部的表达式的值，则先将局部导数和局部的表达式的值计算出来，再逐层往上最终得到整体的导数值。

        通过对表达式二叉树的深度遍历，计算出每个结点的局部表达式的导数和原始值，最终根节点处的导数值即为整体表达式的导数值。

 60分代码及debug

​
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int MOD=1e9+7;
int n,m,cnt;
string st, str;
struct nodes
{
	int id=-1;
	int l=-1;//该节点左孩子的id 为-1时则表示没有左孩子 
	int r=-1;//该节点右孩子的id 为-1时则表示没有右孩子 
	string po;//该节点自身表示的字符串 
	int d_v;//以该节点为根的子树的表达式的导数的值 
	int pre_v;//以该节点为根的子树的表达式 直接带入各个数值的原始值 
}node[240];
stack<nodes> s;
int x[110];
int xnum;
int qmi(int a,int b)
{
	a=a%MOD;
	b=b%MOD;
	int res=0;
	res=(a*b)%MOD;
	return res;
}
void dfs(int cnt)
{
	if(node[cnt].l==-1&&node[cnt].r==-1)//没有左右孩子是叶子结点 
	{//叶子结点一定是某自变量或常数 

		string temp=node[cnt].po;
		if(temp[0]=='x')//是自变量 
		{
	        //算出xi的i 
			int num=temp[1]-'0';
			

			if(num==xnum)//是设定的需要求偏导的那个自变量 
			{
				
				node[cnt].d_v=1;
				node[cnt].pre_v=x[num]; 
			}
			else//其他自变量 当做常数处理 
			{
			
				node[cnt].d_v=0;
				node[cnt].pre_v=x[num];
			}
		}
		else//是常数本身 
		{
				node[cnt].d_v=0;
				node[cnt].pre_v=atoi(temp.c_str());
		}
		return;
	}
	else
	{//不是叶子结点一定是某个符号 

		string temp=node[cnt].po;
		dfs(node[cnt].l);//遍历左右，用于下方计算出该结点子树形成的表达式的导数值和原始值 
		dfs(node[cnt].r);

		if(temp=="+")
		{
			node[cnt].pre_v=(node[node[cnt].l].pre_v+node[node[cnt].r].pre_v)%MOD;
			node[cnt].d_v=(node[node[cnt].l].d_v+node[node[cnt].r].d_v)%MOD;
		}
		else if(temp=="-")
		{
			node[cnt].pre_v=node[node[cnt].l].pre_v-node[node[cnt].r].pre_v;
			node[cnt].d_v=node[node[cnt].l].d_v-node[node[cnt].r].d_v;
		}else if(temp=="*")
		{
			int a=node[node[cnt].l].d_v;
			int b=node[node[cnt].l].pre_v;
			int c=node[node[cnt].r].d_v;
			int d=node[node[cnt].r].pre_v;
			node[cnt].pre_v=qmi(b,d);
			node[cnt].d_v=qmi(a,d)+qmi(b,c);
		}
	}
}
signed main()
{
    cin >> n >> m;
    getchar();
    getline(cin, st);
    stringstream ss(st);
    while (getline(ss, str, ' ')) {

		if(str=="+"||str=="-"||str=="*")
		{
			cnt++;//给树的结点编号从1开始 
			node[cnt].po=str;
			node[cnt].id=cnt;
			
			//从栈中取出两个元素 设为其左右孩子
			nodes temp1;
			temp1=s.top();
			s.pop();
			node[cnt].r=temp1.id;

			nodes temp2;
			temp2=s.top();
			s.pop();
			node[cnt].l=temp2.id;

			s.push(node[cnt]);

		}
		else
		{
			cnt++;
			node[cnt].po=str;
			node[cnt].id=cnt;	
			s.push(node[cnt]);		
		}
	}
	//root的id=cnt

	int root=cnt;
	while(m--)
	{
		xnum=0;
		memset(x,0,sizeof(x));
		cin>>xnum;
		for(int i=1;i<=n;i++) cin>>x[i];
		
		dfs(root);//以编号为cnt的结点为根节点遍历每个结点计算导数的值 
		int res=node[root].d_v;
		
		res=res%(MOD);
		if(res>=0)
		{
			cout<<res<<endl;			
		}
		else
		{
			res+=MOD;
			cout<<res<<endl;
		}
	}

	return 0;
}

​

        在判定某自变量是第几个自变量时，由于自变量个数会超过9个，因此直接将temp[1]-‘0’的值当做自变量xi的i值是错误的，下方代码进行了改正。

100分代码

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int MOD=1e9+7;
int n,m,cnt;
string st, str;
struct nodes
{
	int id=-1;
	int l=-1;//该节点左孩子的id 为-1时则表示没有左孩子 
	int r=-1;//该节点右孩子的id 为-1时则表示没有右孩子 
	string po;//该节点自身表示的字符串 
	int d_v;//以该节点为根的子树的表达式的导数的值 
	int pre_v;//以该节点为根的子树的表达式 直接带入各个数值的原始值 
}node[240];
stack<nodes> s;
int x[110];
int xnum;
int qmi(int a,int b)
{
	a=a%MOD;
	b=b%MOD;
	int res=0;
	res=(a*b)%MOD;
	return res;
}
void dfs(int cnt)
{
	if(node[cnt].l==-1&&node[cnt].r==-1)//没有左右孩子是叶子结点 
	{//叶子结点一定是某自变量或常数 

		string temp=node[cnt].po;
		if(temp[0]=='x')//是自变量 
		{
			//int num=temp[1]-'0';
			
			//算出xi的i 
			int num=0;
			for(int i=1;temp[i]!='\0';i++)
			{
				num=num*10+(temp[i]-'0');
			}
			
			if(num==xnum)//是设定的需要求偏导的那个自变量 
			{
				
				node[cnt].d_v=1;
				node[cnt].pre_v=x[num]; 
			}
			else//其他自变量 当做常数处理 
			{
			
				node[cnt].d_v=0;
				node[cnt].pre_v=x[num];
			}
		}
		else//是常数本身 
		{
				node[cnt].d_v=0;
				node[cnt].pre_v=atoi(temp.c_str());
		}
		return;
	}
	else
	{//不是叶子结点一定是某个符号 

		string temp=node[cnt].po;
		dfs(node[cnt].l);//遍历左右，用于下方计算出该结点子树形成的表达式的导数值和原始值 
		dfs(node[cnt].r);

		if(temp=="+")
		{
			node[cnt].pre_v=(node[node[cnt].l].pre_v+node[node[cnt].r].pre_v)%MOD;
			node[cnt].d_v=(node[node[cnt].l].d_v+node[node[cnt].r].d_v)%MOD;
		}
		else if(temp=="-")
		{
			node[cnt].pre_v=node[node[cnt].l].pre_v-node[node[cnt].r].pre_v;
			node[cnt].d_v=node[node[cnt].l].d_v-node[node[cnt].r].d_v;
		}else if(temp=="*")
		{
			int a=node[node[cnt].l].d_v;
			int b=node[node[cnt].l].pre_v;
			int c=node[node[cnt].r].d_v;
			int d=node[node[cnt].r].pre_v;
			node[cnt].pre_v=qmi(b,d);
			node[cnt].d_v=qmi(a,d)+qmi(b,c);
		}
	}
}
signed main()
{
    cin >> n >> m;
    getchar();
    getline(cin, st);
    stringstream ss(st);
    while (getline(ss, str, ' ')) {

		if(str=="+"||str=="-"||str=="*")
		{
			cnt++;//给树的结点编号从1开始 
			node[cnt].po=str;
			node[cnt].id=cnt;
			
			//从栈中取出两个元素 设为其左右孩子
			nodes temp1;
			temp1=s.top();
			s.pop();
			node[cnt].r=temp1.id;
            //注意此处要先设为右孩子，后出栈的设为左孩子
			nodes temp2;
			temp2=s.top();
			s.pop();
			node[cnt].l=temp2.id;

			s.push(node[cnt]);

		}
		else
		{
			cnt++;
			node[cnt].po=str;
			node[cnt].id=cnt;	
			s.push(node[cnt]);		
		}
	}
	//root的id=cnt

	int root=cnt;
	while(m--)
	{
		xnum=0;
		memset(x,0,sizeof(x));
		cin>>xnum;
		for(int i=1;i<=n;i++) cin>>x[i];
		
		dfs(root);//以编号为cnt的结点为根节点遍历每个结点计算导数的值 
		int res=node[root].d_v;
		
		res=res%(MOD);
		if(res>=0)
		{
			cout<<res<<endl;			
		}
		else
		{
			res+=MOD;
			cout<<res<<endl;
		}
	}

	return 0;
}
其他注意点：

        在计算各个结点处对应的导数和原始值的过程中就一边进行对题目要求的MOD进行取模，防止数字过大。

"daydayup"
关注
  • 10
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
csp13-21java100分.zip
12-08
CSP 考试是由中国计算机学会组织的一项旨在评测个人计算机编程和算法应用能力的考试，它分为 CSP-J（初级）和 CSP-S（提高级）两个级别。这个压缩包可能包含了从 2013 年至 2021 年历年 CSP 考试中一、二题的题目、...
CSP-S2022估分器
09-20
CSP-S2022估分器是一款专为计算机软件能力认证（CSP）初赛设计的评分工具。CSP，全称为“Certified Software Professional”，是由中国计算机学会（CCF）主办的一项针对计算机编程技能的评级考试。该考试旨在评估...
CSP 202309-3 梯度求解
weixin_44340652的博客
09-26 3552
比较简单的一道题，只需要把其他变量对应的数代入进去，在栈里存表达式的要求的偏导的系数和对应的指数就可以，最后就变成处理关于偏导变量的一元多次方程的导数，就比较简单了。比赛的时候没有调出来，好亏。这种解法对于变量无关的求解是可行的，实际上利用二叉树的方法更好。
CCF-CSP 202309-3 梯度求解 大模拟 100分
qq_62571980的博客
10-21 4339
有了上面的基础，多项式可以以下面的方式进行存储，vector中每一项都是一个乘积项，且默认这些乘积项都由加号连接，那么遇到减号该如何处理呢？选择合适的数据结构是解决大模拟问题的关键，一旦选定了合适的数据结构，题目也就迎刃而解了，通常使用的数据结构有栈（STL模板库），map,vector,结构体及其复杂嵌套结构，具体应该如何嵌套，具体情况具体分析。到此为止，最关键的数据结构模型已经建立好了，随后的多项式加法，减法，乘法，乃至求偏导，都是非常基础的模拟过程，具体请详见下方完整代码。
CCF-CSP真题《202309-3 梯度求解》思路+python，c++满分题解
weixin_53919192的博客
11-05 8546
本文介绍了CCF-CSP真题《202309-3 梯度求解》思路+python，c++满分题解，希望对你有帮助~
CSP-202309-3-梯度求解
fzy2003的博客
02-22 1033
并在每次迭代中处理这些子字符串。这种处理方式对于解析和执行基于逆波兰表示法（RPN）的算术表达式非常有效，因为它允许程序按照操作的顺序（从左到右）逐步计算表达式的结果。：从标准输入读取一个字符。这里它的作用可能是用来“吃掉”（消耗）前一个输入后留下的。是一个流类，可以像输入/输出流一样操作字符串。，然后将读取到的文本（不包括换行符）存储到之前声明的字符串变量中。（如空格、制表符等）来拆分原始字符串。，要符合导数运算的规则。循环的作用是逐个读取。由于求的是导数，这里的。常数，类似于非变量的。
2023.09 CCF-CSP 认证 第三题 《梯度求解》 80分
qq_29817727的博客
09-26 1947
刚开始完全没想法，索性用导数定义求解，将自变量和自变量减0.000001带入求函数值，用两个函数值之差除以0.000001，最后结果再四舍五入。后来想到最适合求导的就是多项式形式，一个数组就是一个多项式，下标就是一项的指数，值就是一项的系数。后缀表达式求值的方法也可以用来算多项式。不过还是没AC，只有80分，不知道问题在哪。
CCF-CSP真题《202309-4 阴阳龙》思路+python，c++满分题解
weixin_53919192的博客
11-05 5629
本文介绍了CCF-CSP真题《202309-4 阴阳龙》思路+python，c++满分题解，希望对你有帮助~~
CCF CSP题解：坐标变换（其二）（202309-2）
曹无悔的博客
10-12 6092
使用前缀和的思想，求解CCF CSP认证的试题《坐标变换（其二）（202309-2）》。
CCF-CSP真题《202309-2 坐标变换（其二）》思路+python，c++满分题解
weixin_53919192的博客
11-05 1万+
本文介绍了CCF-CSP真题《202309-2 坐标变换（其二）》思路+python，c++满分题解，希望对你有帮助~~
CCF-CSP真题《202309-5 阻击》思路+ c++题解
weixin_53919192的博客
12-17 3151
本文用c++题解了CCF-CSP真题《202309-5 阻击》，希望对你有帮助~~
CCF-CSP真题《202309-1 坐标变换（其一）》思路+python，c++，java满分题解
热门推荐
weixin_53919192的博客
09-25 1万+
本文为大家题解了CCF-CSP真题《202309-1 坐标变换（其一）》，希望对你有帮助~~
信息学奥赛入门组复赛CSP 2021-J2题目详解.pdf
07-13
【信息学奥赛入门组复赛CSP 2021-J2题目详解】涉及到的是信息学竞赛中的算法问题，主要包括两道题目：分糖果（candy）和插入排序（sort）。以下是这两道题目的详细解析及解决方案： **第一题 分糖果（candy）** ...
信息学奥赛入门组复赛CSP 2019-J2题目详解.pdf
07-13
本篇内容主要涵盖了2019年CCF非专业级软件能力认证CSP-J/S的第二轮入门组CSP 2019-J2的三道题目，包括T1数字游戏、T2交通换乘和T3纪念品的解题策略和代码实现。 **T1 数字游戏** 这是一个关于字符串处理的问题。...
信息学奥赛入门组复赛CSP 2020-J2题目详解.pdf
07-13
本资料主要涉及的信息学竞赛是CCF非专业级软件能力认证（CSP-J/S），特别是2020年第二轮认证的入门组（CSP 2020-J2）。这个竞赛旨在培养青少年的计算机科学能力，通过解决一系列编程问题来评估他们的逻辑思维和技术...
算法的学习笔记—删除链表中重复的结点（牛客JZ76）
最新发布
业精于勤，荒于嬉
08-14 1201
本文介绍了一种删除排序链表中重复节点的高效算法。通过递归的方式，算法能够有效地删除所有重复节点，并保留链表中仅出现一次的节点。此方法的实现简单且易于理解，同时在时间和空间复杂度方面表现优异，适用于处理大规模数据链表中的重复节点问题。😁热门专栏推荐想学习vue的可以看看这个java基础合集数据库合集redis合集nginx合集linux合集手写机制微服务组件spring_尘觉springMVCmybits🤔欢迎大家加入我的社区尘觉社区。
算法：魔法字典
Garson的博客
08-12 568
设计一个使用单词列表进行初始化的数据结构，单词列表中的单词。如果给出一个单词，请判定能否只将这个单词中字母换成另一个字母，使得所形成的新单词存在于你构建的字典中。实现dictionarydictionarysearchWordtruefalse// 返回 False// 将第二个 'h' 替换为 'e' 可以匹配 "hello" ，所以返回 True// 返回 False// 返回 FalsedictionarybuildDictsearch100search。
备战秋招60天算法挑战，Day17
溜达虎的博客
08-13 396
给定两个字符串 和 ，编写一个函数来判断 是否是 的字母异位词。注意： 若 和 中每个字符出现的次数都相同，则称 和 互为字母异位词。举个例子： 视频题解 有效的字母异位词 思路来源 思路来源 思路解析 方法一 排序 分别对和进行排序，然后比较排序后的字符串是否相等。 java代码 python代码 复杂度分析 时间复杂度： 快排的时间复杂度是O(nlogn)，其中是字符串的长度。空间复杂度： 快排的空间复杂度是O(logn)。假设， 。 java代码 python代码 复杂度分析 时
【LeetCode每日一题】——623.在二叉树中增加一行
IronmanJay的博客
08-14 1133
给定一个二叉树的根 root 和两个整数 val 和 depth ，在给定的深度 depth 处添加一个值为 val 的节点行。 注意，根节点 root 位于深度 1 。 加法规则如下: 给定整数 depth，对于深度为 depth - 1 的每个非空树节点 cur ，创建两个值为 val 的树节点作为 cur 的左子树根和右子树根。 cur 原来的左子树应该是新的左子树根的左子树。 cur 原来的右子树应该是新的右子树根的右子树。 如果 depth == 1 意味着 depth - 1 根本没有深度
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值