源代码见文末 部分源代码: % 符号变量 两种表达方式 a=sym('a'); class(a); syms b; b; % 符号常量 c=sym('); c; % 符号表达式 三种表达方式 f1='3*x+4'; % char类型,但是可以参与运算.所以这种形式是可以直接用的,比较方便 syms x f2=3*x+4; % sym类型 f3=sym('3*y+4'); % sym类型 % 符号四则运算 fadd1 = f1 + f2; fmu1 = f1*f2; % 符号表达式的化解 syms…

P1689 方程求解 题目描述 给一个方程,形如X+Y=Z或X-Y=Z.给出了其中两个未知数,请求出第三个数.未知数用‘?’表示,等式中也许会出现一些多余的空格. 输入输出格式 输入格式: 一行,方程. 输出格式: ‘?’代表的值 输入输出样例 输入样例#1: 复制 样例输入1 1+2=? 样例输入2 3 +? = 2 输出样例#1: 复制 样例输出1 3 样例输出1 -1 说明 0<=X,Y,Z<1,000,000,000 模拟 #include<cstdio> #include…

P1689 方程求解 题目描述 给一个方程,形如X+Y=Z或X-Y=Z.给出了其中两个未知数,请求出第三个数.未知数用‘?’表示,等式中也许会出现一些多余的空格. 输入输出格式 输入格式: 一行,方程. 输出格式: ‘?’代表的值 输入输出样例 输入样例#1: 复制 样例输入1 1+2=? 样例输入2 3 +? = 2 输出样例#1: 复制 样例输出1 3 样例输出1 -1 说明 0<=X,Y,Z<1,000,000,000 思路:模拟 #include<cstdio> #incl…

FESTUNG模型介绍-1.对流方程求解 1. 控制方程 对流问题中,控制方程表达式为 \[\partial_t C + \partial_x (u^1 C) + \partial_y (u^2 C) = 0, \quad \mathrm{in} \; \Omega \] 其中边界处包含 Dirichlet 和 Neumann 边界条件,分别为 \[\begin{array}{cl} C = C_D & \mathrm{on} \; \partial \Omega_D, \\ - \nabla…

Time Limit: 1 second Memory Limit: 128 MB [问题描述] 要求Xi(i = 1,2,3,4)是一个[-T..T]中的整数,满足方程AX1 + BX2 + CX3 + DX4 = P的解有多少组? [输入格式] 六个整数分别表示P.T.A.B.C.D. [输出格式] 一个整数表示方程的解的个数. |P| < 109 |A| < 104 |B| < 104 |C| < 104 |D| < 104 0 <= T <= 500 Sa…

话说2015.11.06 ,北京下了第一场雪.16年的今天没下雪,但是雾霾还是不小的,帮媳妇整理她工作时,出现了下面的需求,便想到使用PHP来写程序来进行求解. [需求] 1. 给出一个平均值X,反过来求出来,得到这个平均值X的三个数X1 ,X2, X3,最大值与最小值的差值要小于0.4(X1-X3都是保留1位小数的数) 2. 这三个数X1, X2, X3代表了三组数.满足下面的公式: X1 = [(m1 - m2)/(m1 - m0) ] * 100 (@1); m0, m1, m2三个数的边…

题目:http://acm.timus.ru/problem.aspx?space=1&num=1132 题意:就是给出方程,p为素数,求在区间内的解. 这个思路很简单,详见:http://algo.ftiasch.com/tag/number-theory/ 一开始TLE,原因是我用了二分加法,以后记住:二分加法是适合很大数的,比较小的数就直接乘,不然数据多了可能TLE. #include <stdio.h> #include <stdlib.h> #include &l…

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3329 题意:现在有三个骰子,分别有k1,k2和k3面,面上的点就是1~ki.每次扔骰子,如果这三个骰子的值分别对应为a,b,c,那么将值初始化为0,否则就将三个骰子的点值和相加.求大于等于n的扔骰子次数期望. 思路: 这道题目主要在于推公式,看着别人的题解想了好久. 先设$E(i)$为此时和为i时还需要的期望,易得(1):$E(i)=\sum (E(i+k)P(k))+E…

任意门:http://acm.hdu.edu.cn/showproblem.php?pid=3292 No more tricks, Mr Nanguo Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)Total Submission(s): 587    Accepted Submission(s): 400 Problem Description Now Sailormoon…

//穷举法! /* ====================================================== 题目:求佩尔方程x*x-73*y*y=1的解. ====================================================== */ #include<stdio.h> #include<math.h> int main(void) { int x,y; double t; for(y=1;y<=10000000;y+…