12097 - Pie
Time limit: 3.000 seconds
My birthday is coming up and traditionally I'm serving pie. Not just one pie, no, I have a number&N&of them, of various tastes and of various sizes.&F&of my friends are coming to my party and each
of them gets a piece of pie. This should be one piece of one pie, not several small pieces since that looks messy. This piece can be one whole pie though.
My friends are very annoying and if one of them gets a bigger piece than the others, they start complaining. Therefore all of them should get equally sized (but not necessarily equally shaped) pieces, even if this leads
to some pie getting spoiled (which is better than spoiling the party). Of course, I want a piece of pie for myself too, and that piece should also be of the same size.
What is the largest possible piece size all of us can get? All the pies are cylindrical in shape and they all have the same height 1, but the radii of the pies can be different.
One line with a positive integer: the number of test cases. Then for each test case:
One line with two integers&N&and&F&with 1 ≤&N, F&≤ 10000: the number of pies and the number of friends.One line with&N&integers&ri&with 1 ≤&ri&≤ 10000: the radii of the pies.
For each test case, output one line with the largest possible volume&V&such that me and my friends can all get
a pie piece of size&V. The answer should be given as a floating point number with an absolute error of at most 10-3.
Sample Input
1 4 2 3 4 5 6 5 4 2
Sample Output
思路：注意这题是让我们求小数形式的最优解而不是即约分数(最简分数)形式
那有更好的方法求吗？——回想高中解超越方程近似根所用的方法（二分法），我们就有了如下的代码：
/*0.029s*/
#include&cstdio&
#include&cmath&
#include&algorithm&
const double PI = acos(-1.0);
const int maxn = 10005;
double A[maxn];
inline bool ok(double area)
int sum = 0;
for (int i = 0; i & i++)
sum += floor(A[i] / area);
return sum &= f + 1;
int main(void)
scanf(&%d&, &T);
while (T--)
scanf(&%d%d&, &n, &f);
double maxa = -1.0;
for (int i = 0; i & i++)
scanf(&%d&, &r);
A[i] = PI * r *
maxa = max(maxa, A[i]);
double L = 0.0, R =
while (R - L & 1e-5)//二分查找计算，达到要求精度时停止查找
double M = (L + R) / 2;
if (ok(M))
printf(&%f\n&, L);
* 以上用户言论只代表其个人观点，不代表CSDN网站的观点或立场
访问：252824次
积分：9816
积分：9816
排名：第482名
原创：704篇
转载：42篇
评论：97条
(67)(76)(26)(29)(97)(125)(154)(160)(13)有哪位神会解超越方程么_数学吧_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0成为超级会员，使用一键签到本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：265,968贴子：
有哪位神会解超越方程么收藏
超越方程这门槛太难跨了
除了些特殊的外就只有Mr. Computer回了
不会 二分法什么的
登录百度帐号我的游戏推荐游戏
后查看最近玩过的游戏
使用签名档&&
为兴趣而生，贴吧更懂你。或您所在位置： &
&nbsp&&nbsp&nbsp&&nbsp
托梁拔柱节点试验研究与理论分析论文.pdf74页
本文档一共被下载：
次 ,您可免费全文在线阅读后下载本文档
文档加载中...广告还剩秒
需要金币：200 &&
--优秀硕士毕业论文，完美PDF格式，可在线免费浏览全文和下载,支持复制编辑,可为大学生本专业本院系本科专科大专和研究生学士硕士相关类学生提供毕业论文范文范例指导,也可为要代写发表职称论文的提供参考！！
你可能关注的文档：
··········
··········
节点研究一直是钢结构研究中的重点。根据节点处传递荷载的情况、所采用
的连接方法及其细部构造，按照节点的力学特性，可分为刚性节点、半刚性节点
和铰接节点。目前国内对于节点的研究主要集中在半刚性节点领域，对铰接节点
涉及较少，同时尚没有明确的关于工程实际中铰接节点的判定标准。本文以天津
钢管公司460管加工车间托梁拔柱工程为背景，围绕其中的新型改造节点一托梁
拔柱节点进行了以下几方面的试验研究和理论分析工作。
简述了节点设计原则，分析了节点研究现状，简介了本文选题的工程背景，
提出了本文的主要工作。
介绍了节点的连接方法以及国内外的节点连接分类标准，并根据本文所研究
的方向，着重介绍并分析了五种工程中常见的铰接节点的连接形式及优缺点，为
本文的研究工作奠定了理论基础。
结合工程背景提出了托梁拔柱节点改造的形式和构造，在ANSYS里建立改
造节点缩尺模型，重点考虑了柱头板对于提高节点承载力的有利影响，运用非线
性有限元法对节点进行理论分析，明确了节点的破坏形态，得到了节点的理论承
载力，为试验提供了参考依据。
为了对改造后节点的实际性能有更为深入的了解，设计并进行了改造节点缩
尺模型的拟静力试验。对节点的滞回曲线、骨架曲线、能量耗散系数、刚度退化
系数等指标进行了详细分析，并与有限元结果进行对比，说明同刚性节点相比，
改造后节点的耗能能力差，刚度退化严重，表明改造后的托梁拔柱节点具有明显
的铰接特性。
根据上述理论分析和试验研究，提出了简化的节点弹簧模型，得到了节点转
正在加载中，请稍后...matlab方程求教 - 程序语言 - 小木虫 - 学术 科研 第一站
&& 查看话题
matlab方程求教
函数部分：
function epsiul_m =epsiul( Lambda )
syms&&epsiul_m w
%%%%%%% parameters of gold
epsiul_L=5.9673;
w_D=*pi*1e12;
r_D=15.92*2*pi*1e12;
omigal_L=650.07*2*pi*1e12;
T_L=104.86*2*pi*1e12;
deltepsiul=1.09;
&&= solve('epsiul_m=epsiul_L-w_D^2/(w*(w+j*r_D))-deltepsiul*omigal_L^2/((w^2-omigal_L^2)+j*T_L*w)','w=2*pi*c/(Lambda*(sqrt(real(epsiul_m)/epsiul_0)))');
epsiul_0=(1/(36*pi))*1e-9;
程序部分：
Lambda=1.55e-6;
epsiul_m =epsiul( Lambda);
n_gold=sqrt(epsiul_m/epsiul_0);
为什么运行后最后出来的结果还是包含epsiul_m的表达式呢，怎么修改直接出来一个数值，而不包含epsiul_m和其它已知数的表达式，，比如最后结果连pi都未带入，本来结果应该就是一个数值呢，而不是一个表达式呢。我不知道我哪里出错了。
也没人理俺:o Matlab中solve函数主要是用来求解线性方程组的解析解或者精确解。主要用法有四种，返回参数至多只有一个，所以楼主的“ = solve('epsiul_m=epsiul_L-w_D^2/(w*(w+j*r_D))-deltepsiul*omigal_L^2/((w^2-omigal_L^2)+j*T_L*w)','w=2*pi*c/(Lambda*(sqrt(real(epsiul_m)/epsiul_0)))');”我不知道是怎么运行结果的，是不应该有两个返回参数的吧。 : Originally posted by 兵封年代 at
Matlab中solve函数主要是用来求解线性方程组的解析解或者精确解。主要用法有四种，返回参数至多只有一个，所以楼主的“ = solve('epsiul_m=epsiul_L-w_D^2/(w*(w+j*r_D))-deltepsiul*omigal_L^2/((w^2-omigal_L^2)+ ... syms epsiulm w
epsiul_L=5.9673;
w_D=*pi*1e12;
r_D=15.92*2*pi*1e12;
omigal_L=650.07*2*pi*1e12;
T_L=104.86*2*pi*1e12;
deltepsiul=1.09;
epsiul_0=(1/(36*pi))*1e-9;
Lambda=1.55e-6;
= solve('epsiulm-epsiul_L-w_D^2/(w*(w+i*r_D))-deltepsiul*omigal_L^2/((w^2-omigal_L^2)+i*T_L*w)=0','(w*Lambda/(2*pi*c))^2-real(epsiulm)/epsiul_0=0');
n_gold=sqrt(epsiulm/epsiul_0);
不用函数其实应该是这样子的。但是出来的结果总是还包含epsiulm。本来结果应该直接就是一个数值的 : Originally posted by susulucky at
syms epsiulm w
epsiul_L=5.9673;
w_D=*pi*1e12;
r_D=15.92*2*pi*1e12;
omigal_L=650.07*2*pi*1e12;
T_L=104.86*2*pi*1e12;
deltepsiul=1.09;
epsiul_0=(1/(36*pi))*1e-9;
Lambda=1.5 ... 我说的意思是solve函数没有这种用法，它的返回参数至多只有一个。
var cpro_id = 'u1216994';
欢迎监督和反馈：本帖内容由
提供，小木虫为个人免费站点，仅提供交流平台，不对该内容负责。欢迎协助我们监督管理，共同维护互联网健康，如果您对该内容有异议，请立即发邮件到
联系通知管理员，也可以通过QQ周知，我们的QQ号为：8835100
我们保证在1个工作日内给予处理和答复，谢谢您的监督。
小木虫，学术科研第一站，为中国学术科研研究提供免费动力
欢迎监督，发现不妥请立即
E-mail： & QQ：8835100解超越方程，新手请各位帮帮忙，十分感谢！_matlab吧_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0成为超级会员，使用一键签到本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：54,016贴子：
解超越方程，新手请各位帮帮忙，十分感谢！收藏
方程为x*sqrt(1-y)=2*atan(1/0.98*sqrt(y/(1-y)))第一问：怎么求解这个方程，求代码。第二问：如何做出xy的图，x的取值范围在0到20，求代码。明天就要交作业了，恳请大神帮个忙！
这个挺简单啊！
solve命令直接打就行——
y=solve('x*sqrt(1-y)=2*atan(1/0.98*sqrt(y/(1-y)))','x')x=1:20;for i=1:length(x)
z(i)=subs(y,x(i));endzz是复数，无法绘出平面图
y=solve('x*sqrt(1-y)=2*atan(1/0.98*sqrt(y/(1-y)))','x')x=1:.01:20;for i=1:length(x)z(i)=subs(y,x(i));endzplot3(x,real(z),imag(z))
登录百度帐号我的游戏推荐游戏
后查看最近玩过的游戏
使用签名档&&
为兴趣而生，贴吧更懂你。或