C - 柱爷的下凡

Time Limit: 400/400MS (Java/Others)     Memory Limit: 65535/65535KB (Java/Others)

Submit  Status

下凡的柱爷想只创造3种硬币,要求这3种硬币可以组成[1,N][1,N]的所有整数,并且表示一个[1,N][1,N]中的整数所用硬币的平均数量最少.

你能帮柱爷解决他的小小问题吗?

Input

第一行一个正整数TT,表示测试组数.

接下来TT行,每行一个正整数NN.

数据保证:

  • 1≤T≤2001≤T≤200

  • 1≤N≤2001≤N≤200

Output

输出一共有TT行,每行三个整数AA BB CC表示柱爷希望的三种不同硬币的面值.如果有多组解,请保证A尽可能小,如果仍有多组解,请保证B尽量小,如果仍有多组解,请保证C尽量小. A<B<CA<B<C

Sample input and output

Sample Input Sample Output
1
1
1 2 3
1
7
1 2 5

Hint

当N=7N=7时

  • 11元,需要11个11元
  • 22元,需要11个22元
  • 33元,需要11个11元和11个22元
  • 44元,需要22个22元
  • 55元,需要11个55元
  • 66元,需要11个11元和11个55元
  • 77元,需要11个22元和11个55元
  • 平均需要约1.57
function testmain1
close all
T = 14560000;
T = 0.8*T;
%北京市总人口(Total Population)
%1456.00万人(2003年)
%已确诊病例累计  现有疑似病例  死亡累计  治愈出院累计
d=[339  402  18  33
482  610  25  43
588  666  28  46
693  782  35  55
774  863  39  64
877  954  42  73
988  1093  48  76
1114  1255  56  78
1199  1275  59   78
1347  1358  66  83
1440  1408  75  90
1553  1415  82  100
1636  1468  91  109
1741  1493  96  115
1803  1537  100  118
1897  1510  103  121
1960  1523  107  134
2049  1514  110  141
2136  1486  112  152
2177  1425  114  168
2227  1397  116  175
2265  1411  120  186
2304  1378  129  208
2347  1338  134  244
2370  1308  139  252
2388  1317  140  257
2405  1265  141  273
2420  1250  145  307
2434  1250  147  332
2437  1249  150  349
2444  1225  154  395
2444  1221  156  447
2456  1205  158  528
2465  1179  160  582
2490  1134  163  667
2499  1105  167  704
2504  1069  168  747
2512  1005  172  828
2514  941  175  866
2517  803  176  928
2520  760  177  1006
2521  747  181  1087
2522  739  181  1124
2522  734  181  1157
2522  724  181  1189
2522  718  181  1263
2522  716  181  1321
2522  713  183  1403
2523  668  183  1446
2522  550  184  1543
2522  451  184  1653
2522  351  186  1747
2523  257  186  1821
2523  155  187  1876
2522  71  187  1944
2522  4  189  1994
2522  3  189  2015
2521  3  190  2053
2521  5  190  2120
2521  4  191  2154
2521  3  191  2171
2521  3  191  2189
2521  2  191  2231
2521  2  191  2257
2521  2  191  2277];
S = T - d(:,1)-d(:,3)-d(:,4);I = d(:,1)+d(:,2); % 感染人数
S1 = d(:,2); % 现有疑似病例
R = d(:,3)+d(:,4);
t = (1:length(I))';
subplot(1,2,1);
h=plot(t,I,'o-',t,S1,'*');
legend(h,'I','S')xdata = [I(1:end-1) S(1:end-1)]
ydata = diff(I)./diff(t)
%ydata = ydata((1:30),:);coef = lsqcurvefit(@ffun,[1 1],xdata,ydata)
format long
k = coef(1)
h = coef(2)
%z=kIS-hI
% lsqcurvefit xdata(:,1) -> I, xdata(:,2)->S
[T,Y]=ode23(@dfun,[0:60],[S(1) I(1) R(1)]);
%Y(:,2) 为I的数值解
subplot(1,2,2);
%plot(xdata,ydata,'ro',xdata,ffun(coef,xdata),'k*');
plot(t,I,'ro',T',Y(:,2)','k*')function ydata=ffun(coef,xdata)
k = coef(1); h = coef(2);
ydata = k*xdata(:,1).*xdata(:,2)-h*xdata(:,1);function dy=dfun(t,y) %y (1)-->S, y(2)-->I
dy = zeros(2,1); S=y(1);  I = y(2); R=y(3);
k = 1.738382084343759e-05;
h = 0.004435622880419;
dy(1) = -k.*I.*S;
dy(2)=k.*I.*S-h.*I;
dy(3)=h.*I;

CDOJ1323柱爷的下凡相关推荐

  1. 【CDOJ 1323】柱爷的下凡

    很有趣的题,写代码来打表然后o(1)可过. #include<cstdio> #include<cstring> #include<iostream> using ...

  2. cdoj 1354 柱爷很忙

    题目传送门 柱爷很忙 Time Limit: 1000/1000MS (Java/Others)     Memory Limit: 65535/65535KB (Java/Others) Submi ...

  3. CDOJ 1347柱爷的矩阵(二维dp)

    柱爷的矩阵 Time Limit: 125/125MS (Java/Others)     Memory Limit: 65535/65535KB (Java/Others) Submit  Stat ...

  4. UESTC-2016dp专题 F 柱爷与三叉戟不得不说的故事

    柱爷与三叉戟不得不说的故事 Time Limit: 500/500MS (Java/Others)     Memory Limit: 65535/65535KB (Java/Others) 题目链接 ...

  5. CDOJ 1355 柱爷与三叉戟不得不说的故事(枚举子集状压dp)

    柱爷与三叉戟不得不说的故事 Time Limit: 500/500MS (Java/Others)     Memory Limit: 65535/65535KB (Java/Others) Subm ...

  6. CDOJ 1321柱爷的恋爱 (区间dp)

    柱爷的恋爱 Time Limit: 1000/1000MS (Java/Others)     Memory Limit: 65535/65535KB (Java/Others) Submit  St ...

  7. 2016 UESTC Training for Dynamic Programming F - 柱爷与三叉戟不得不说的故事 压位dp

    F - 柱爷与三叉戟不得不说的故事 Time Limit: 500/500MS (Java/Others)     Memory Limit: 65535/65535KB (Java/Others) ...

  8. cdoj 1322 柱爷把妹(吃惊高清重制版)

    题目传送门 柱爷把妹(吃惊高清重制版) Time Limit: 2000/2000MS (Java/Others)     Memory Limit: 125000/125000KB (Java/Ot ...

  9. 2017四川省赛总结

    可能这个省赛是很纠结的东西啊,好像提前一个月 其它省都举办了省赛,然后sc慢一点?刚好省赛又在蓝桥杯国赛后面一周,所以五月份就很忙的啊,最主要的是有一门考试刚好跟这两个时间冲突,刚开始跟蓝桥杯冲突,后 ...

最新文章

  1. [保护模式]PAE模式
  2. Java学习小程序(4)数列求和
  3. 俺的新书《Sencha Touch实战》终于出版了
  4. Drools DMN最新开源引擎性能改进
  5. uniapp动态修改样式_掌握Photoshop图层样式技术
  6. 苹果cmsV10高仿草民电影网在线影视网站模板 带手机版
  7. min java_LeetCode算法题-Min Stack(Java实现)
  8. 记一次Project插件开发
  9. SQL Server单表查询语句
  10. Docker 私有仓库搭建
  11. mysql网络异常_网络连接配置出现异常_网络连接配置无法修复_Mysql网络连接的性能配置项...
  12. java中重写equals方法
  13. Atitit 发帖机实现(2)---usrQBN2243 文本解析到对象协议规范
  14. SQL语句:查询多表更新数据
  15. NAS和私有云盘的区别?1篇文章说清楚
  16. 朴素贝叶斯之邮件分类
  17. DELL笔记本E5400刷BIOS激活win7
  18. 努力和家境哪个更重要
  19. 王者荣耀国际服——传说对决
  20. 【软考-软件设计师】(五).计算机网络

热门文章

  1. python 代理ip池_GitHub - xuan525/proxy_pool: Python爬虫代理IP池(proxy pool)
  2. 扶蕾wlt~学习计划——MATLATB
  3. 关闭windows server服务器的密码过期提示
  4. Spring5学习详细笔记
  5. 软件测试种类(一):包括功能性测试,可靠性测试,强度测试等
  6. 如何将原始SNP信息转化为0,1,2的矩阵形式
  7. ExcludeClipRect和无闪烁图像
  8. Ubuntu16.04无法卸载VSCode
  9. Delphi的RTTIVMT
  10. 互联网金融学习总结(4)——大数据风控的九种维度