- 注册时间
- 2018-12-8
- 最后登录
- 1970-1-1
- 威望
- 星
- 金币
- 枚
- 贡献
- 分
- 经验
- 点
- 鲜花
- 朵
- 魅力
- 点
- 上传
- 次
- 下载
- 次
- 积分
- 3549
- 在线时间
- 小时
|
发表于 2019-12-2 21:07:07
|
显示全部楼层
本帖最后由 dlpg070 于 2019-12-2 21:10 编辑
根据王守恩27#数据,归纳出如下2个公式
1 通项公式
2 递推公式
给出代码和模拟27#结果
-------------
- Clear["Global`*"];
- time1 = AbsoluteTime[];
- t[n_, k_] := (n - k)*n^(k - 1);
- (*
- A058127=Table[t[n,k],{n,1,10},{k,0,n-1}]//Flatten; (*Jean-Fran?ois \
- Alcover,Dec 03 2013*)
- A058127a=Table[{"n=",n," k=",k," ",t[n,k],"\n"},{n,1,10},{k,0,n-1}]
- Print["A058127 ",A058127]
- *)
- sol = Table[t[n, n - 1], {n, 2, 12}];
- time2 = AbsoluteTime[];
- Print["通项耗时 ", time2 - time1, " 秒"]
- Print["通项 输出 ", sol]
- time1 = AbsoluteTime[];
- a[2] = 1;
- a[n_] := a[n] =
- Module[{sum = 0},
- sum = (n - 1)*a[n - 1] +
- Sum[Binomial[n - 1, i]*t[n - 1, n - i - 1], {i, 2, n - 1}];
- sum];
- sol1 = Table[a[n], {n, 2, 12}];
- time2 = AbsoluteTime[];
- Print["递推耗时 ", time2 - time1, " 秒"]
- Print["递推 输出 ", sol]
- Print["模拟王守恩27#: "]
- time1 = AbsoluteTime[];
- Do[(* n *)
- str = ToString[Binomial[n - 1, 1]] <> "x" <>
- ToString[t[n - 1, n - 1 - 1]];
- Do[(*i *)
- str = str <> "+" <> ToString[Binomial[n - 1, i]] <> "x" <>
- ToString[t[n - 1, n - i - 1]], {i, 2, n - 1}];
- str = str <> "=" <> ToString[a[n]];
- Print[str];
- Print[" " <> ToString[n] <> " 岛 " <> ToString[n - 1] <> " 桥有" <>
- ToString[a[n]] <> "种建桥方案"];
-
- Print[" "],
- {n, 2, 12}]
- time2 = AbsoluteTime[];
- Print["模拟耗时 ", time2 - time1, " 秒"]
- Print["--- end ---"]
复制代码
计算结果:
通项耗时 0.0010001 秒
通项 输出 {1,3,16,125,1296,16807,262144,4782969,100000000,2357947691,61917364224}
递推耗时 0.0030002 秒
递推 输出 {1,3,16,125,1296,16807,262144,4782969,100000000,2357947691,61917364224}
模拟王守恩27#:
1x1=1
2 岛 1 桥有1种建桥方案
2x1+1x1=3
3 岛 2 桥有3种建桥方案
3x3+3x2+1x1=16
4 岛 3 桥有16种建桥方案
4x16+6x8+4x3+1x1=125
5 岛 4 桥有125种建桥方案
5x125+10x50+10x15+5x4+1x1=1296
6 岛 5 桥有1296种建桥方案
6x1296+15x432+20x108+15x24+6x5+1x1=16807
7 岛 6 桥有16807种建桥方案
7x16807+21x4802+35x1029+35x196+21x35+7x6+1x1=262144
8 岛 7 桥有262144种建桥方案
8x262144+28x65536+56x12288+70x2048+56x320+28x48+8x7+1x1=4782969
9 岛 8 桥有4782969种建桥方案
9x4782969+36x1062882+84x177147+126x26244+126x3645+84x486+36x63+9x8+1x1=100000000
10 岛 9 桥有100000000种建桥方案
10x100000000+45x20000000+120x3000000+210x400000+252x50000+210x6000+120x700+45x80+10x9+1x1=2357947691
11 岛 10 桥有2357947691种建桥方案
11x2357947691+55x428717762+165x58461513+330x7086244+462x805255+462x87846+330x9317+165x968+55x99+11x10+1x1=61917364224
12 岛 11 桥有61917364224种建桥方案
模拟耗时 0.0030002 秒
--- end --- |
|