😊出会いの泉😊

最大は9回で、9回になるのは
0149,0589,0941,0985,1094,1490,4109,4901,
5098,5890,8509,8905,9014,9058,9410,9850
の16個です。
ただし、abcd,bcda,cdab,dabc,dcba,cbad,badc,adcbが同じ回数になりますので、
本質的には0149と0589の2個ですね。

(追記)
abcd+efgh=9999のときabcdの回数とefghの回数は同じなので、本質的には0149の1通りだけでした。(∵0149+9850=9999)

(追々記)
5桁で試したら、「1回以下で終わるもの」と「無限に終わらないもの」しかないようでした。
7桁も同じです。奇数桁では自明な解を除き0にならないのかも知れません。
6桁は最大4回(例:014523)、8桁は最大22回(例:00012448)、10桁は最大4回(例:0143014523)でした。

(さらに追記)
よく考えたら奇数桁では2回以上の解はないですね。
例えば5桁でもし2回以上の解があったとすると
最後が00000→その前はaaaaa (a=1～9)
その前がbcdefとすると
c=b±a, d=c±a, e=d±a, f=e±a, b=f±aなので
b=b±a±a±a±a±a
±a±a±a±a±a=0
∴±1±1±1±1±1=0
これはあり得ませんので奇数桁では2回以上の解はなく、
よって無限回を除外すると1回(全桁同じ数字から開始)が最大となります。

出所は、記憶が定かではありませんが、25年くらい前の、数学セミナーだと思います。
設問は、「４つの数（桁数関係なく）からはじめて、差をとると0になり、その操作の回数を10回以上にしてください。」
　0を除く一桁の数で、10回以上可能だったような。二桁かも。記憶は誤りでしたか？
　選択の数の桁数を増やせば、いくらでも回数を増やすことができる？

小学生に、出題する予定でした。

「四角形の数」数学の部屋
サイトがありました。

上限が9より大きくてよければ、10回以上は可能です。
「その操作の回数を10回以上」とのことなので最初の状態はカウントしません。
(0,2,6,13)
→(2,4,7,13)
→(2,3,6,11)
→(1,3,5,9)
→(2,2,4,8)
→(0,2,4,6)
→(2,2,2,6)
→(0,0,4,4)
→(0,4,0,4)
→(4,4,4,4)
→(0,0,0,0)
最小数と最大数の差が12以下のとき10回未満となります。
2桁の最大は13回(例:0,7,20,44)
3桁の最大は19回(例:0,81,230,504)
4桁の最大は25回(例:0,927,2632,5768)
でした。

数字を、大きくしても、なかなか、10回以上は、見つかりにくいです。
A〈B〈C〈D　かつ　A+B+C＝Dが、良い条件のようです。
他に、見つける、条件、ありませんか？

紹介して頂いた貴重な論文を拝見させて頂きました。
目が回る様な論理の展開で一つの式で表現するためには
大変な考察が必要なことが実感できました。

一般にO(0,0),P(n,m)の2点を結ぶ直線の下方(直線上を含む)の領域
だけを通過する格子路でOからPまでの最短路の総数G(n,m)を求める
プログラムをらすかるさんのアイデアをお借りして以前作成していた
のを思い出しました。
以下がそのプログラム(PARI/GPでのコード)と結果になります。
なお\記号は複数行に渡る記述のための繋ぎのためのものです。

gp > G(n,m)={M=matrix(n+1,m+1,i,j,if(j==1,1,i==1 && j>1,0))};\
for(x=2,n+1,for(y=2,m+1,if(m*(x-1)<n*(y-1),next,\
M[x,y]=M[x-1,y]+M[x,y-1])));M[n+1,m+1]

gp > for(n=2,9,print1(n"=>");for(m=1,30,print1(G(n,m)","));print)
2=>1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,11,11,12,12,13,13,14,14,15,15,16,
3=>1,2,5,5,7,12,12,15,22,22,26,35,35,40,51,51,57,70,70,77,92,92,100,117,117,
126,145,145,155,176,
4=>1,3,5,14,14,23,30,55,55,76,91,140,140,178,204,285,285,345,385,506,506,593,
650,819,819,938,1015,1240,1240,1396,
5=>1,3,7,14,42,42,66,99,143,273,273,364,476,612,969,969,1197,1463,1771,2530,
2530,2990,3510,4095,5481,5481,6293,7192,8184,10472,
6=>1,4,12,23,42,132,132,227,377,525,728,1428,1428,2010,2803,3504,4389,7084,
7084,9097,11654,13793,16380,23751,23751,28931,35246,40356,46376,62832,
7=>1,4,12,30,66,132,429,429,715,1144,1768,2652,3876,7752,7752,10659,14421,
19228,25300,32890,53820,53820,67860,84825,105183,129456,158224,231880,231880,
278256,
8=>1,5,15,55,99,227,429,1430,1430,2529,3978,7229,9690,14985,21318,43263,43263,
61600,82225,121637,148005,199238,254475,420732,420732,543806,672452,900239,
1043460,1307742,
9=>1,5,22,55,143,377,715,1430,4862,4862,8398,15090,22610,35530,58040,81719,
120175,246675,246675,345345,500449,650325,876525,1220135,1542684,2017356,
3362260,3362260,4289780,5630306,

6=>の場合がat氏の出力と一致すると思います。

一般の長方形の格子路でカタラン路のような数を求める時には一辺が
多くの約数を含むようなものについては一個の式で表すのにはどうしても複雑な場合分けでの式が重なってします。
私も一辺が6のものについてatさんとは異なる式となりましたが何とか式にしてみました。

G6(m)={k=m\6;L=6*k+1;S=(6+m)!/(m!*6!);}\
if(m%6==0,S/L,\
m%6==1,S/(L+6),\
m%6==2,(S-k*(k+1)*(3*k+2)*(6*k+7)*(9*k+7)/40)/(L+6),\
m%6==3,(S-k*(k+1)*(6*k+7)*(28*k^2+61*k+31)/30)/(L+6),\
m%6==4,(S-(k+1)*(3*k+4)*(6*k+7)*(57*k^2+133*k+70)/40)/(L+6),\
m%6==5,S/(L+10))

で100までを出力してみると
gp > for(m=1,100,print1(G6(m)",");if(m%10==0,print))
1,4,12,23,42,132,132,227,377,525,
728,1428,1428,2010,2803,3504,4389,7084,7084,9097,
11654,13793,16380,23751,23751,28931,35246,40356,46376,62832,
62832,73950,87143,97584,109668,141778,141778,162883,187453,206591,
228459,285384,285384,322046,364124,396510,433160,527085,527085,586638,
654240,705789,763686,910252,910252,1002037,1105317,1183487,1270752,1489488,
1489488,1625096,1776599,1890570,2017169,2331924,2331924,2525439,2740354,2901207,
3079140,3518515,3518515,3786757,4083170,4304066,4547556,5145336,5145336,5508104,
5907251,6203610,6529292,7324878,7324878,7805193,8331713,8721393,9148503,10187344,
10187344,10811692,11493880,11997356,12547920,13881945,13881945,14680520,15550580,16191123,

場合分けを一つの式にしてみましたが、あまり綺麗になりませんでした。
G6(m)={k=m\6;}((6+m)!/(m!*6!)-\
((m+1)%6\3)*((m%6+1)*k+(m%6\4)*(13*k+24))*\
((54*k^2+78*k+28)+((m%3+2)%4)*(k^2+11*k+8)+(2-m%6%4%3)*(11*k+9))*\
(k+1)*(6*k+7)/240)/(m+163036%(m%6+9))

管理人様
【更新履歴】にコメント有難うございます。

ここで、全く役に立たない定理を考えてみました。

『頂角が２０度の全ての二等辺三角形において、
「底辺と同じ長さの線分」を、頂点を始点として置くと、
終点と、その対角を結んだ線分との、なす角は１５０度となる』


知っていても役に立たないが、テストの時だけ速攻で使えるかな？
ヒントに「正三角形を作る」が無いとテスト時間内に解くには
時間が足りなくなると思う。
私なら「捨て問」です。

角の大きさ（４８）の問題、随分昔にここの掲示板に同じ問題が出されて、解いた覚えがあるんですよね。
なのでここのサイトのどこかに残っているのではないかと思って探し回ったのですが、
見つけられませんでした。
もし過去の掲示板のログが残っていたら、どこかにあるかも知れません。

らすかるさん
有難うございます。

あれれ？？
ホーム画面の【更新履歴】を見ると、
初回の投稿が、令和８年　４／２４になっているので・・・・・


私の大きな勘違いでしたら申し訳ありません　ｍ（＿）ｍ

すみません
４／２５の間違えでした。

らすかるさん
申し訳ありませんでした。

更新履歴の管理人様のコメントを見て、
掲示板の全ての投稿が履歴に
載っている訳ではなかったのですね。

このサイトを、よく理解していなくてスミマセンｍ（＿）ｍ

ちゃちゃっとプログラムを作って（最下段左端＜最下段右端という条件を付けて）調べたところ
（例示されたものも含めて）以下の4通りになりました。

___3
__4,7
_5,9,2
6,1,10,8

___3
__5,2
_4,9,7
6,10,1,8

___4
__2,6
_5,7,1
8,3,10,9

___4
__5,1
_2,7,6
8,10,3,9

ならば
___0
__2,8
_7,5,3
1,6,9,4
の様に
下の隣り合う2数の和の下桁の数が上の段にあり、積み上げ終わると
0～9の数が一通り揃う。
ということになる配列は他にあるか？

全部で以下の4通りだと思います。

___0
__2,8
_7,5,3
1,6,9,4

___0
__4,6
_9,5,1
2,7,8,3

___0
__2,8
_7,5,3
6,1,4,9

___0
__4,6
_9,5,1
7,2,3,8

では、最初の差分の方式で
1段(1のみ): 1通り
2段(1～3): 2通り
3段(1～6): 4通り
4段(1～10): 4通り
となりますが、5段(1～15)では何通りでしょう？

一通りのみでは。
_____5
____4,9
___7,11,2
__8,1,12,10
6,14,15,3,13

では6段には存在するか？
存在しないならその証明は？

5段の1通りは正解です。
証明はわかりませんが、6段・7段・8段では解はありませんでした。
「6段以上では解はない」という可能性もありますが、
さすがに6・7・8だけでは何とも言えないですね。
ちなみに8段の全探索には半日かかりました。

6段での証明が気になったので色々調べてみたら
Shichermanというこのパズルを提出した人らしい人物が
mod 2
では|a-b|≡a+b (mod 2)
から,
6個の異なる整数
a,b,c,d,e,fから
a+b,b+c,c+d,d+e,e+f
a+2*b+c,b+2*c+d,c+2*d+e,d+2*e+f
a+3*b+3*c+d,b+3*c+3*d+e,c+3*d+3*e+f
a+4*b+6*c+4*d+e,b+4*c+6*d+4*e+f
a+5*b+10*c+10*d+5*e+f
と和を作っていき、それまでのすべて現れる総和が
6*a+20*b+34*c+34*d+20*e+6*f
なのでその数字は偶数であることになる。
一方1~21(6段では全部の数は1+2+･･･+6=21)
の数の総和は21*22/2=231
で奇数である。
mod2では奇数、偶数は一致するはずなのでこれは矛盾し
如何なる6個の数でも構成は不可能となる。

で示していた。
そこで7段ではと思い
a,b,c,d,e,f,g
で生まれてくる2数の和による（mod 2での考察)構成での総和をみると
7*a+27*b+55*c+69*d+55*e+27*f+7*g≡1 (mod 2)
一方
28*29/2=406≡0 (mod 2)
より7段でも矛盾
8段では
a,b,c,d,e,f,g,h
からは総数
8*a+35*b+83*c+125*d+125*e+83*f+35*g+8*h≡0 (mod 2)
36*37/2=666≡0 (mod 2)
従って8段はこの手では証明ができないことになる。

さらに気になったのでAIを利用して尋ねると
n>5では一切存在できないことがHerbert Taylorにより証明が与えられており
かなりテクニカルなもので単純なパリティ計算だけではすまず、より精巧な
不変量や構造解析を行うとある。

なおこれが2018段においては不可能であることの証明を問う問題が
数学オリンピックに出題さているという。

やはり6段以上では解はなかったのですね。
2018段の問題をいきなり出されても解けないだろうなぁ・・・

今さらですが、7段の証明がちょっと違うような気がします。
（私の勘違いでしたらご容赦下さい）
> 7*a+27*b+55*c+69*d+55*e+27*f+7*g≡1 (mod 2)
これはa,b,c,d,e,f,gのうち奇数が奇数個なら成り立ちますが、
奇数が偶数個の場合は
7*a+27*b+55*c+69*d+55*e+27*f+7*g≡0 (mod 2)
となるのではないでしょうか。

a,b,c,･･･,gの奇遇によって変化してしまいますね。
単に奇数の係数が奇数個であったので1と判断してしまっておりました。
らすかるさんの指摘どうりですね。

Shichermanというと、ジッヒャーマンダイスのShichermanでしょうか。

左からでも右からでも同じ並びの数は「回文数」といいます。
5回を超えるものは山ほどありますが、以下は例です。

6回で終わるもの
1069, 1079, 1159, 1169, 1249, 1259, 1339, 1349, 1429, 1439,
1519, 1529, 1609, 1619, 1699, 1709, 1789, 1799, 1879, 1889,
(以下略)

7回で終わるもの
1394, 1484, 1574, 1664, 1754, 1844, 1898, 1934, 1988, 1992,
1994, 1999, 2393, 2483, 2573, 2663, 2753, 2843, 2897, 2933,
(以下略)

8回で終わるもの
1993, 1995, 2992, 2994, 3991, 3993, 4990, 4992, 5991, 6990,
8059, 8149, 8239, 8329, 8419, 8509, 8599, 8689, 8779, 8869,
(以下略)

9回で終わるもの
1397, 1487, 1577, 1667, 1757, 1847, 1937, 2396, 2486, 2576,
2666, 2756, 2846, 2936, 2999, 3395, 3485, 3575, 3665, 3755,
(以下略)

10回で終わるもの
9059, 9149, 9239, 9329, 9419, 9509, 9599, 9689, 9779, 9869, 9959

11回で終わるもの
7069, 7159, 7249, 7339, 7429, 7519, 7609, 7699, 7789, 7879,
7969, 8068, 8158, 8248, 8338, 8428, 8518, 8608, 8698, 8788,
(以下略)

12回で終わるもの
2069, 2159, 2249, 2339, 2429, 2519, 2609, 2699, 2789, 2879,
2969, 3068, 3158, 3248, 3338, 3428, 3518, 3608, 3698, 3788,
(以下略)

13回で終わるもの
1797, 1887, 1894, 1977, 1984, 2796, 2886, 2893, 2976, 2983,
3795, 3885, 3892, 3975, 3982, 4794, 4884, 4891, 4974, 4981,
(以下略)

14回で終わるもの
1991, 2990

15回で終わるもの
1998, 2997, 3996, 4995, 5994, 6079, 6169, 6259, 6349, 6439,
6529, 6619, 6709, 6799, 6889, 6979, 6993, 7078, 7168, 7258,
(以下略)

16回で終わるもの
1496, 1586, 1676, 1766, 1856, 1897, 1946, 1987, 2495, 2585,
2675, 2765, 2855, 2896, 2945, 2986, 3494, 3584, 3674, 3764,
(以下略)

17回で終わるもの
1792, 1882, 1972, 2791, 2881, 2971, 3790, 3880, 3970

18回で終わるもの
1798, 1888, 1978, 2797, 2887, 2977, 3796, 3886, 3976, 4795,
4885, 4975, 5794, 5884, 5974, 6793, 6883, 6973, 7792, 7882,
(以下略)

20回で終わるもの
6999, 7998, 8039, 8129, 8219, 8309, 8399, 8489, 8579, 8669,
8759, 8849, 8939, 8997, 9038, 9128, 9218, 9308, 9398, 9488,
(以下略)

21回で終わるもの
1297, 1387, 1477, 1567, 1657, 1747, 1837, 1927, 2296, 2386,
2476, 2566, 2656, 2746, 2836, 2926, 3295, 3385, 3475, 3565,
(以下略)

10000回でも終わらないもの
1495, 1497, 1585, 1587, 1675, 1677, 1765, 1767, 1855, 1857,
1945, 1947, 1997, 2494, 2496, 2584, 2586, 2674, 2676, 2764,
(以下略)
4桁で10000回でも終わらないものは、数回で
52514,83127,96558,97768,109989
のいずれかになる数です。

ピラミッド数（仮称）
121＝11×11
12321＝111×111
1234321＝1111×1111
123454321＝11111×11111
掛け算でつくれますが、
上の数を、操作１，２，３で作ることを考えてみました。
任意の数で、できるかは不明です。

10進法なら
111111111×111111111＝12345678987654321
までですが、例えば11進法にすると、
1111111111×1111111111＝123456789A987654321
となりますね。

操作1，2，3で
121＝56＋65
24300＋342＝24642（偶数なので）
12321＝24642÷2
2464000＋4642＝2468642
1234321＝2468642÷2
5の時は、繰り上がるので、難しいが、可能。

196問題
196＋691＝887
887＋788＝1675
1675＋5761＝7436
7436÷2＝3718
3718＋8173＝11891
11891＋19811＝31702
31702÷2＝15851

123454321=(246850000+58642)/2
12345654321=(24685650000+5658642)/2
1234567654321=(2468567650000+567658642)/2
123456787654321=(246856787650000+56787658642)/2
12345678987654321=(24685678987650000+5678987658642)/2

二桁の10からは、不思議なことが起こりました。

12345678910987654321を作るために
パターンに従がって、２倍した
24685678910987650000と逆転数5678901987658642を足して２で割れば、出来なくて
逆転数の代わりに、5678910987658642だとうまくいきます。

（246851250000＋51258642）÷2＝123451254321　？！

三ケ桁の数（100～999）について調べました。
A＝｛数字が、回文数｝
B＝｛逆転数との和が、回文数｝
C＝｛逆転数との和が偶数で、２で割ると、回文数｝
A（90）＋B（210）＋C（284）＝584　重複は避けました。

三桁の数の場合、逆転数との和が、偶数ならば、2で割った数は、回文数になるようです。どうでしょうか？4桁以上は、どうでしょうか？
操作３を使えば、４桁も、回文数になる回数が、大分減少します。

四ケタの数　ABCDで、A＋D,B＋Cが、共に偶数であれば、
2で割って、回文数
五ケタの数　ABCDEで、A＋E,B＋Dが、共に偶数であれば、
2で割って、回文数

試して頂いて有難うございます。
前提条件を、端折ってしまいました。
元の数ABCDの逆転数DCBAを、足して2で割った結果が、回文数に、なります。
５ケタも同様です。

誤差はほぼe/(2*5^(3^84))なので
小数点以下約8368428989068425943817590916445001887164桁(≒8.37×10^39桁)正しい
ということになるかと思います。

「小町算で無理数近似」の記事の話ですかね？

(1+0.2^(9^(7*6)))^5(^(3^84))

これは

(1+0.2^(9^(7*6)))^(5^(3^84))

ということですか？

はい。

ご教示をありがとうございました。

(1+1/N)^N
の形になっているのだとようやく気が付きました。
頭ワルすぎ‼️

三角形ABCの、外接円の半径をR、内接円の半径をｒとするとき、
　　Ｒ≧2ｒ　が成り立つ。等号は、正三角形のとき

三角形ＡＢＣの内部の点をＰ、フェルマー点をＦとするとき、
　　ＡＰ＋ＢＰ＋ＣＰ≧ＡＦ＋ＢＦ＋ＣＦ　
　　　　が成り立つ。　

△ABCの外接円のAにおける接線とBCの延長との交点（交わるとき）
をDとすると
BD：CD＝AB＾２：AC＾２　が成り立つ。

△ABCの外接円のB,Cにおける接線の交点をPとし、
直線APとBCの交点をDとすると
BD：CD＝AB＾２：AC＾２　が成り立つ。

△BCDに対して、∠BCD＝∠BADとなる、平行四辺形をとる。
△BCDの外接円の頂点Cにおける接線が直線AB、ADと交わる点をP、Qとすると
PC：CQ＝AP＾２；AQ＾２,
BP：DQ＝AP＾３：AQ＾３
が成り立つ。3乗が珍しい（私見）

成分が異なるだけでよいなら変なものも作り放題なので、
勝手に整数かつ非ゼロという条件を付けてみます。


目標の固有値が対角成分に並ぶ三角行列Aと正則行列Pを用いると、PAP^(-1)は目標の固有値を持つ。

Pの各成分が整数で行列式が+1または-1ならばP^(-1)の各成分も整数になる。

あとは適当に試して成分がすべて異なるようにすればよい。


(1)
A=[[-3,1],[0,7]], P=[[1,0],[-1,1]]
としてみると、
PAP^(-1)=[[-2,1],[9,6]]

(2)
A=[[1,3,2],[0,2,1],[0,0,3]], P=[[1,0,0],[3,1,0],[0,2,1]]
としてみると、
PAP^(-1)=[[4,-1,2],[12,-3,7],[18,-6,5]]

(3)
A=[[-4,3,2],[0,8,1],[0,0,9]], P=[[1,0,0],[1,1,0],[0,2,1]]
としてみると、
PAP^(-1)=[[-3,-1,2],[-9,5,3],[6,-6,11]]

こんな便利な方法があるのですね。
目的の固有方程式を満たすようにある意味力ずくで探していました。

そこで(3)の結果の行列S=[[-3,-1,2],[-9,5,3],[6,-6,11]]を使わせてもらって
f1(x)=1/312*x^3+3/104*x^2-29/156*x
と
f2(x)=1/312*(313*x^3--4047*x^2+1190*x+89856)
の2つの関数において、それぞれ
f1(S),f2(S)を計算させると結果は共に3次の正方行列M1,M2に集約されますが(f2の定数項では3次の単位行列を補う。)
それぞれM1,M2の固有値は何でしょうか？

とりあえず普通に計算しました。


Sの固有値λの固有ベクトルをvとする。
nを自然数とするとき、
S^nv
=S^(n-1)Sv
=S^(n-1)(λv)
=λ*S^(n-1)v
=λ*S^(n-2)Sv
=λ*S^(n-2)(λv)
=λ^2*S^(n-2)v
…
=λ^(n-1)*Sv
=λ^(n-1)*(λv)
=λ^n*v
なので、vはS^nの固有ベクトルでその固有値はλ^nである。


[1]
M1=1/312*S^3+3/104*S^2-29/156*S
に右からvを掛けると
M1v
=(1/312*S^3+3/104*S^2-29/156*S)v
=1/312*S^3v+3/104*S^2v-29/156*Sv
=1/312*λ^3*v+3/104*λ^2*v-29/156*λ*v
=(1/312*λ^3+3/104*λ^2-29/156*λ)*v
なので、vはM1の固有ベクトルであり、その固有値は
1/312*λ^3+3/104*λ^2-29/156*λ
である。
Sの固有値-4,8,9を代入するとM1の固有値は1,2,3となる。


[2]
単位行列をIとする。
M2=1/312*(313*S^3-4047*S^2+1190*S+89856*I)
に右からvを掛けると
M2v
=1/312*(313*S^3-4047*S^2+1190*S+89856*I)v
=1/312*(313*S^3v-4047*S^2v+1190*Sv+89856*v)
=1/312*(313*λ^3*v-4047*λ^2*v+1190*λ*v+89856*v)
=1/312*(313*λ^3-4047*λ^2+1190*λ+89856)*v
なので、vはM2の固有ベクトルであり、その固有値は
1/312*(313*λ^3-4047*λ^2+1190*λ+89856)
である。
Sの固有値-4,8,9を代入するとM2の固有値は1,2,5となる。

計算ありがとうございます。
[2]ではSの固有値9ではM2の固有値は3となりませんか？

たまたまフロベニウスの定理というものに出会い、本当にこんなことが起こるのか？
と思って色々固有値をもつ行列Sを使って実験をしている中で
f(x)の関数で作り上げるf(S)の行列Mの固有値をこちらが指定できるものに動かすことができる
f(x)はどんな関数として設定しておけばいいのかを探すのに
ラグランジュの補間法からのf1(x)
ファンデルモンドの行列式からのf2(x)
で験していたのがこの計算でした。
どうしてこんなことが成り立つのかは正しくりらひいさんが示されることで納得できました。

実験してみて一般に(x1,y1),(x2,y2),(x3,y3)を通る3次関数は
原点を通るものと、原点を通らないものと2通り存在できることが起こるんですね。

計算ミスしていたようです。
確かに、[2]ではSの固有値9でM2の固有値が3となりました。
失礼しました。


(-4,1), (8,2), (9,3) を通る3次関数は
f(x)=ax^3+bx^2+cx+d
とおいて、
f(-4)=1, f(8)=2, f(9)=3 より
-64a+16b-4c+d=1, 512a+64b+8c+d=2, 729a+81b+9c+d=3
を連立して解いた答え
b=-13a+11/156, c=4a-31/156, d=288a-12/13
（WolframAlphaに解いてもらった）
を使用して
f(x)=ax^3+(-13a+11/156)x^2+(4a-31/156)x+(288a-12/13)
とすればいいです。

a=1/312 を代入すれば f1(x) が得られ、
a=313/312 を代入すれば f2(x) が得られます。

地道に計算しようと思って式を立てて計算を進めたら、結構あっさり解けました。
A=[1,1;1,9], B=[1,4;5,2] (Pari/GPの形式) ですね。

(追記)
私が解いた何の工夫もない方法を整理すると
A=[a,b;b,c], B=[d,e;f,g]とおくと、条件から
(1)ad+bf=6 (2)ae+bg=6 (3)bd+cf=46 (4)be+cg=22
(5)ad+be=5 (6)bd+ce=37 (7)af+bg=7 (8)bf+cg=23
(7)-(2)からa(f-e)=1, (8)-(4)からb(f-e)=1, (3)-(6)からc(f-e)=9
なのでa=b=1,c=9,f=e+1。(2)(4)(5)に代入して
(2)e+g=6 (4)e+9g=22 (5)d+e=5
(4)-(2)からg=2, (2)からe=4, (5)からd=1, f=e+1からf=5
∴(a,b,c,d,e,f,g)=(1,1,9,1,4,5,2)

|A||B| = -144
なので、A,Bは正則である。

B^(-1) = [[a, b], [c, d]]
とおく。

一つ目の式に右からB^(-1)を掛けると
A = [[6a+6c, 6b+6d], [46a+22c, 46b+22d]]
となり、
二つ目の式に左からB^(-1)を掛けると
A = [[5a+7b, 37a+23b], [5c+7d, 37c+23d]]
となる。

6a+6c = 5a+7b, 6b+6d = 37a+23b, 46a+22c = 5c+7d, 46b+22d = 37c+23d
を連立して解くと、p,qを任意の実数として
a = 7p-6q, b = p, c = q, d = 46p-37q
となるので、
A = [[42p-30q, 282p-222q], [322p-254q, 1058p-814q]]
と書ける。

Aが対称行列のとき、
282p-222q = 322p-254q
を解くと、rを任意の実数として
p = 4r, q = 5r
と書けるので、
A = [[18r, 18r], [18r, 162r]]
および
B^(-1) = [[-2r, 4r], [5r, -r]]
となる。よって、
B = [[1/(18r), 4/(18r)], [5/(18r), 2/(18r)]]
となる。

Aの各成分が自然数になるのはmを自然数として r = m/18 のときであり、
Bの各成分が自然数になるのはnを自然数として r = 1/(18n) のときである。
m/18 = 1/(18n) より mn = 1 なので、m,nがともに自然数になるのは m = n = 1 のときだけである。
このとき r = 1/18 であり、
A = [[1, 1], [1, 9]], B = [[1, 4], [5, 2]]
となる。



><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><



|A||B| = -144
なので、A,Bは正則である。

A^(-1) = [[a, b], [c, d]]
とおく。

一つ目の式に左からA^(-1)を掛けると
B = [[6a+46b, 6a+22b], [6c+46d, 6c+22d]]
となり、
二つ目の式に右からA^(-1)を掛けると
B = [[5a+37c, 5b+37d], [7a+23c, 7b+23d]]
となる。

6a+46b = 5a+37c, 6a+22b = 5b+37d, 6c+46d = 7a+23c, 6c+22d = 7b+23d
を連立して解くと、p,qを任意の実数として
a = -46p+37q, b = p, c = q, d = -7p+6q
となるので、
B = [[-230p+222q, -254p+222q], [-322p+282q, -154p+138q]]
と書ける。

Aが対称行列のとき、A^(-1)も対称行列になるので、rを任意の実数として p = q = r と書けて
B = [[-8r, -32r], [-40r, -16r]]
および
A^(-1) = [[-9r, r], [r, -r]]
となる。よって、
A = [[-1/(8r), -1/(8r)], [-1/(8r), -9/(8r)]]
となる。

Aの各成分が自然数になるのはmを自然数として r = -1/(8m) のときであり、
Bの各成分が自然数になるのはnを自然数として r = -n/8 のときである。
-1/(8m) = -n/8 より mn = 1 なので、m,nがともに自然数になるのは m = n = 1 のときだけである。
このとき r = -1/8 であり、
A = [[1, 1], [1, 9]], B = [[1, 4], [5, 2]]
となる。