約 1,580,651 件
https://w.atwiki.jp/monosepia/pages/6447.html
いろいろブログ&サイト【航空機・船舶・自動車・鉄道など乗り物】 ● X-30(航空機)〔Wikipedia〕 ● 宇宙船 - GATAG | Free Images and Photos 2.0 ■ 非科学JAXA!成層圏内. 建屋宇宙混同放射線線量基準! 「~ようこそ!GALLERY*DISCOVERYへ~ (2013.6.16)」より / (※ 詳細はブログ記事で) ■ 宇宙航空機 「史跡探訪」より / 米国は秘密の新世代宇宙航空機を1990年代初頭から開発してきた。ロシアの情報筋は、米国がある航空基地で超秘密機の試験を行っていると暴露した。ロシアの代理人がこの新しい物体を自分の目で見て、写真を撮ろうと試みたが全部失敗した。米国人はその秘密兵器に信じられないくらい厳重に警備し、その試験は夜だけ行っている。しかし、その代理人は米国が、夜明け前の女神にちなんでオーロラと呼ぶこの新航空機に関するいくつかの情報を得た。 この謎の航空機は40キロ以上の高高度を飛ぶ。最新の追跡機でさえこんな高度を飛べない。このオーロラ機はまた科学者が言う6Mという途方もないスピードを出す。1Mは秒速340メートルに相当する。換言すると、オーロラは秒速2キロメートルで飛ぶがこれが制限速度ではない。この速度は極めて強力なエンジンと、胴体や翼に施された特別な設計によって達成された。 +続き この航空機はステルス技術を基にして造られた。その技術は米国防衛企業が以前に選定したものであり、地上や機上レーダが探知するのを困難にする。 本物のオーロラ機の写真はUFOの写真より珍しい。ロシアはこの謎の米航空機の写真を一枚だけ持っている。その写真はロシアの情報官が、オーロラ機が空中給油を受けている瞬間を撮影したものである。オーロラに給油する「空飛ぶ燃料タンク」(KC-135)の前方を、スパイ機がロシアを東から西に横切ったがレーダに探知されなかった。 米国はロケットに似た新しい航空機を開発した。この航空機はその目を将来に向けた米国の防衛科学者達が開発した。最も強力な追跡機は35~40キロメートルの高度を飛ぶ。ロシアはこの点では米国と太刀打ちできる。 ロシアはSU-2711,MiG29およびMiG31を持っている。しかしこれらの航空機は地球に近い環境内でのみ飛行できる。そこと宇宙空間の間にまだ開きがある。それが地球に近い環境の兵器開発が防衛にきわめて関連する理由である。それを先に手に入れるものが勝者になる。 ここで我々は全く新しい航空機(偵察機、戦闘機、旅客機等)を話題にしている。それらは大気圏から宇宙空間に入り込み、また地球に戻ることができる。米国はすでにこの種の航空機に長年取り組んできており、幾つかの「スターウオーズ」計画を進めてきた。米国はそんな航空機を長年設計試験してきた。彼らはすでにX-33、X-34、X-37及びX-38を持っている。 X-33:この超音速航空機は最新のチタン合金を用いて設計された。酸素・水素ジェットエンジンを装備している。この航空機は毎秒4,420メートルの速度に達した。グラフアイトとセラミック(熱的安定)技術がこのエンジンに使用されている。これら技術の幾つかはロシアから買うか、盗んだものである。最高飛行高度は120キロメートル。 X-34:より近代的な超音速機。この機体は独特の電子デバイスを備え、特殊な熱遮蔽コーテイングがされている。その特殊な特徴はこの機体が船からではなく、搭載母機から発進することにある。X-34は秒速5,440メートルに達し、最高高度は150キロメートルである。 X級の超音速機は無人と有人機や、垂直と水平離陸機もある。公式報道によれば、かような航空機は貨物を宇宙に運搬するのに必要な費用を大幅に減らせるという。現在では貨物1キロを運搬するのに2万2千ドルかかる。これらの新しい航空機は軍事衛星を軌道に投入するというような防衛目的に確かに役立つ。超音速宇宙機は超強力な望遠鏡やその他の装備品を搭載できる。科学者や工学者はこの機体にミサイル、レーザ兵器を搭載しようとしている。 別の言い方をすると、米国は初めて宇宙軍を創設しようとしている。それは近地球軌道を完全に支配しようとする。どんな対ミサイル防衛施設もそれらが全くの射程外にあるので、そのような航空機を打ち落とすことはできない。それに反して、米国の宇宙軍は地球上のいかなる目標もたたくことができる。 ソ連の工学者たちは米国に挑戦しようとし、ブランと呼ぶユニークな航空機を設計した。不幸にもこの航空機は成功せず米国との競争に敗れた。ソ連と米国間の軍事対立が終わったとき、ブランはモスクワの文化公園に展示物された。この最初のソ連の宇宙航空機は過去の権力と、こうあってほしい希望との退屈な思い出になった。 ロシア宇宙アカデミーの会員のセルゲイ・コルガノフ陸軍中将は米国が超音速宇宙航空機を開発しているのはロシアにもうひとつの脅威を与えていると語る。「私はロシアがこれら米国の宇宙機と戦えるか疑問に思っている。対ミサイル防衛施設は高度100キロメートルまでの目標を打ち落とすことができる。しかし、これらの新しい航空機はより高空を飛ぶことができる。ロシアは、私ははっきり言えないがそれを防衛する非常によいアイデアを持っている。米国はそれを夢にも見たことがないだろう。しかし、それらはアイデアの域を出るものでない。わが国の科学者たちは米国の宇宙機に対処できるどんなプロジェクトにも携わっていない」と語った。 ■ 夢のような宇宙航空機の試作機開発へ!スカイロンの開発は順調 「ZOOLBOX(2012.4.29)」より (※ 詳細はブログ記事で) ーーーーーーーーーー ■★ 音速の5倍の速さで宇宙に飛び立つ『宇宙航空機スカイロン』 「ロケットニュース24(2010.9.18)」より / 『衝撃映像! ブラジルで捕獲された宇宙人のリアル動画』というニュースが世間を賑わせているが、こちらはもっと現実的なお話。イギリスの宇宙開発企業が、今後10年以内に普通の空港の滑走路から宇宙に飛び立つ『宇宙航空機』を開発すると発表している。 同社によれば、宇宙航空機は水素と酸素を燃料にし、音速の5倍の速さで宇宙へと飛び立つという。しかも、従来のロケットに比べて約15分の1の費用で宇宙旅行を実現できるというのだ。 (※ 以下略、詳細はサイト記事で) .
https://w.atwiki.jp/zayin/pages/1028.html
プレイヤースキル 財力 厨房度 ★ ★★★★ ★★★★★★★★★★★★★★★ ランカークラス Class SS キルクラス Class SS デット数 0 所属部隊名 072st 発言の痛さ 消防 勝ち馬属性 カセより他国援軍のほうがうまい 戦闘スタイル 火皿 タグ 総評 ゲド戦記 本人への要望 どんどんカセを弱体化させてください 本人から すでに弱すぎる件 犬の工作という説も流れたが、乙wikiが水兵や金の派閥等と比べると適当なので、違った模様。 http //zoome.jp/seibou/diary/2/ 072stが所持する唯一のリア♀。 皿キャラはバンク専門らしい 何こいつどうみても小学生 -- 名無しさん (2009-01-10 02 13 14) おい!主戦にこれないkzさっさとキャラデリしろ -- 名無しさん (2009-01-10 02 58 50) 主戦に行ったところで押し引き激しいし追撃少ないからうまみがねぇよ(´ー`) -- 名無しさん (2009-01-10 07 54 10) 本鯖産不燃ゴミはどうせすぐいなくなる あくびでもしつつ眺めてりゃいい -- 名無しさん (2009-01-10 10 49 52) これで釣ってるつもりでいるから困る -- 名無しさん (2009-01-11 19 52 22) それでも釣れてるから本人ご満悦の様子 -- 名無しさん (2009-01-11 22 50 45) 財力4ってこいつ武道エンチャなのか -- 名無しさん (2009-01-13 04 40 24) 犬ってなんですか(´ー`) -- 名無しさん (2009-01-13 18 04 04) (U^ω^)わんわんお! -- 名無しさん (2009-01-18 16 19 35) わんわんお?Fにもそんな名前の部隊がいたような(´ー`)・・・ -- 名無しさん (2009-01-18 19 58 39) 賑わせている? お前が誰かに相手にされてるの見たのは初日だけだわ -- 名無しさん (2009-01-28 18 46 10) ああ、あれ初日で飽きました(´ー`) -- 名無しさん (2009-01-30 17 49 02) お前らちゅっちゅ -- 名無しさん (2009-04-30 02 28 08) お前らちゅっちゅしてあげる -- 名無しさん (2009-05-01 10 27 49) (^q^)<カセ民ちゃんちゅっちゅwwwwwwwww -- 名無しさん (2009-05-02 01 08 44) ☆3はないけどな けんごよりへた -- 名無しさん (2009-05-29 23 31 38) おいおいそこまで言うならキャラデリ杯挑んでこいよ^^ -- 名無しさん (2009-06-01 20 06 11) なんでこっちから挑むんだ・・・?あくまでこっちは自信があるならこいよということだが。 -- 名無しさん (2009-06-06 00 36 20) ↑のキャラデリが怖いってのが見え見えなんですけど? あと日本語ok? -- 名無しさん (2009-06-06 12 44 34) oi、編集してんじゃねーよ -- 名無しさん (2009-06-06 12 46 12) 結局口だけかよざこ -- 名無しさん (2009-06-06 23 05 54) キャラデリしたいのに怖いもないんだが・・・。FPS楽しいです -- 名無しさん (2009-08-16 07 37 01) ろくさすにゃんかわいいお^^ -- 名無しさん (2009-08-16 07 37 29) 口だけの人ってよくいるよね^^ こんな人が2chいるんですよね・・・ -- 名無しさん (2009-08-16 10 21 10) ブンドエリマンの方が断然おもしろい -- 名無しさん (2009-08-16 14 35 14) むしろお前はつまらない -- 名無しさん (2009-08-16 14 35 35) ↑byブンドエリマン -- 名無しさん (2009-08-16 14 36 09) むしろこんなとこ見て書き込みしてる奴は大概口だけだと思うのは僕だけですか? -- 名無しさん (2009-08-17 01 34 11) 直結したらいしいね -- 名無しさん (2010-11-09 01 34 13) ぱぴこ牛乳と直結成功。やるだけやって捨てるのが目に見えてるw -- 名無しさん (2010-11-09 13 27 50) @xRoxasAxelx フォロー待ってる☆ -- 名無しさん (2010-11-15 00 23 02) 名前 コメント
https://w.atwiki.jp/gensouiri/pages/1312.html
ひとことのところだけど、~(人)は嫌いみたいなの書くのはやめた方がいいと思う。作品自体は今後に期待。 -- (名無しさん) 2008-08-14 23 38 02 確かにそうですね。変更しておきました。本編は頑張って行くつもりなのでよろしくお願いします -- (佐藤高雄) 2008-08-19 09 22 08 新作(2話)できましたよ~ 25日の昼すぎにアップします -- (monmon(仮)改め sorai) 2008-08-25 06 05 20 25日はとっくに過ぎましたけどまだUpしないんですか? -- (名無しさん) 2008-08-28 02 07 53 それは多分、soraiさんの幻想入りのことだと思います。 自分は、遅いのでまだ全然できてないっス・・・・・・ -- (佐藤高雄) 2008-08-30 22 40 28 後何日くらいなのかな?出来るの楽しみにしてます!! -- (名無しさん) 2008-10-10 01 29 30 うむ。まだ何日後になるかわかりませんが必ず完成させますよ~。 こういうコメントはめちゃうれしい。作業スピードをチョい上げます。でも、学校の課題ががが -- (佐藤高雄) 2008-10-12 22 01 29 自分の好きな幻想入りの方も新作出したし、自分もそろそろ1話を出すぞっと。これなら今月中に出せそうだ -- (佐藤高雄) 2008-10-21 23 41 06 いよいよ完成ですね! -- (名無しさん) 2008-10-26 21 18 46 いよいよ完成ですね!楽しみにまってます! -- (名無しさん) 2008-10-26 21 19 44 ありがとうございます~。 残念な作品ですが見てやってくださいw -- (佐藤高雄) 2008-10-29 01 08 21 もう一週間たちましたが…許可は出なかったのですか? -- (名無しさん) 2008-11-10 11 59 31 できました!申し訳ないです! -- (佐藤高雄) 2008-11-11 02 05 47 おお!早速見て参ります! -- (名無しさん) 2008-11-11 03 08 28 出来る限りの努力はしました。ただ、他の方の幻想入りを見た後に自分の幻想入りを見ると泣きたくなるOTZ -- (佐藤高雄) 2008-11-12 01 37 41 初めまして。エスコン東方作者のraptorです。 今回はクロス予定の報告をしたいと思い書き込みました。 内容としては ・こちらの航空祭にてヘリ・車両類の展示 以上です。 "幻想入りした兵器"の一部として配置予定です。 よろしければ返信を。 -- (raptor) 2008-11-12 02 22 47 了解しました。ては、貴方の方のページに記入しておきます。 わくわくしてきたぜ! -- (佐藤高雄) 2008-11-12 23 05 35 ま・…まさかだれか考えると思って検索したら・・・ 本当にあった!!!(偶然検索した) -- (名無し) 2009-02-02 01 15 06 紫のキャラが崩壊している… これはファンから苦情来る前に直すべきだろ。 下手すれば荒らされかねない。 -- (名無しさん) 2009-06-09 19 38 09 エターナルおめでとうございます -- (トランペッター) 2009-07-29 18 04 27 続きまだですか? -- (名無しさん) 2009-08-29 22 13 52 てっきり小説の方かと・・・ -- (名無しさん) 2009-09-26 20 13 51 ご自由にどうぞ。てか、是非してください。 ヘリが飛んでるとか、背景の奥のほう で装甲車が走ってるでもいいんで ※ クロスオーバーしてくれた方は出来る限りこちらの動画でもクロスさせ ます。 何言ってるのこの作者・・・誰がわざわざ世界観ぶち壊す様な事すると思ってるんだか。 宣伝してもらおうと必死すぎ。 サークルの方へ、申し訳ありませんが返信が無かったので使わせていただきまし た。 削除の指 -- (名無しさん) 2015-03-02 11 53 42 ミスって送信してしまった。 サークルから返事貰ってないのに使うとか常識に囚われてはいけないのですね。 非表示も当然の結果かも知れないなぁ。 -- (名無しさん) 2015-03-02 11 55 48 序盤からブチぎれ賢者 紫の冷酷さを出そうとしたんだろうけど物凄い小物臭しかしなくてJAXA 人が下手に出りゃ←作者の国語力の低さよ -- (名無しさん) 2017-07-13 19 26 51 そもそもこの手の「妖怪相手に自衛隊ドンパチ」やりたいだけの話は理由付けがお粗末なんだな 紫の暇つぶし~ってアホかいな -- (名無しさん) 2017-07-13 19 35 31 見てないけどエライ荒れてるな 動画自体はもう消えて無いみたいだね ↑のコメを見る限りでは実弾兵器を持った自衛隊が妖怪と戦うわけね 調理次第じゃ面白い話になりそうなんだがなー 暇つぶしは過程として酷すぎるし、幻想郷を我が子の様に愛しているという紫のキャラクターも完全に崩壊している 俺が考えたのはありがちだけど・・・ 博麗大結界周辺の街で行方不明者が際立ってくる ↓ 超常現象に秀でた討伐隊が組まれ幻想郷を辛うじて発見するが、一人(霊能力者系で結界が認知できるなら職業はなんでもいいや)を残し死亡 ↓ 幻想郷発見!警察の特殊部隊が派遣されるが、これも壊滅 ↓ 自衛隊の派遣 俺はこっちの方が自然だと思う ただここまでの話をわかりやすく自然に導入部に入れて尚且つ視聴者を 飽きさせないようにするのは難しいと思うから、そこら辺が腕の見せ所だと思う -- (名無しさん) 2017-07-27 17 53 26 要は完全に戦争状態になるまでの流れが難しいのだと思う 付け加えるなら幻想郷の妖怪が「やりすぎた」って流れを作るまでが中々ストーリに組み込むのが難しそう -- (名無しさん) 2017-07-27 17 58 19 ↑↑それでだいたいあってるよ。 作者の支部絵にはルーミアが隊員の腕を食い千切るだのチルノが氷で戦車破壊だの、 まぁ単なる驚異みたいな感じの構想で描かれてたね 君のストーリー、彼の暇つぶし戦争よりかはまだ自然かもだけど、 そもそもなんで外の人間が勘付くレベルの大量神隠しなんてしちゃうの?という違和感があるね。 gateのパクリみたいな感じだし、東方の意味が、ね。 まぁそもそも論で隠れ里的舞台でわざわざ武力組織との全面戦争をさせる理由付けが難しいし、 割とよくある題材だけど高確率で戦国自衛隊の劣化が出来上がる宿命なのかも。 迷い込みました、でも規模的に理由付けが難しい。 一人だと詰むと思うのか、無限弾薬な能力持たせちゃう歴代作者達 -- (名無しさん) 2017-09-04 23 44 49
https://w.atwiki.jp/ce00582/pages/1100.html
class tax{ public static void main(String args[]){ double th[]= new double[101]; double c[]= new double[101]; double y[]= new double[101]; int opc[]= new int[101]; int opy[]= new int[101]; int opq[]= new int[101]; double u[][][]=new double[101][3][3]; double ww[][][]=new double[100][3][3]; double v[][][][]=new double[100][3][3][21]; int gotoc[][][][]=new int[100][3][3][21]; int gotoy[][][][]=new int[100][3][3][21]; int gotoq[][][][]=new int[100][3][3][21]; double endv[][]=new double[3][3]; int endc[][]=new int[3][3]; int endy[][]=new int[3][3]; int endq[][]=new int[3][3]; int s; double tl; double tr; double ls; double w; double cs; double ys; double b1; double b2; double tr1; double tr2; double tr3; double wel; double maxwel; double maxtl; double maxtr; int t; int n; int n1; int n2; int pp; int q; int j; double c1; double l1; double h; double u1; int ns1; int ns2; int qx; int nx1; int nx2; double vs; double v1; double maxv; double dc; double dy; int e; for (s=1;s 101;s++){ th[s]=0.02*s; } maxtl=0; maxwel=-999; maxtr=0; for (n=20;n 45;n++){ tl=0.01*n; tr1=0.01; tr2=0.02; b1=bud(tl,tr1,th); t=0; while (t 100) { b2=bud(tl,tr2,th); tr3=tr2-b2*(tr2-tr1)/(b2-b1); tr1=tr2; tr2=tr3; b1=b2; if (b2*b2 0.00001)t=1000; t=t+1; } tr=tr2; wel=seekwel(tl,tr,th); if (wel maxwel)maxtl=tl; if (wel maxwel)maxtr=tr; if (wel maxwel)maxwel=wel; } tl=maxtl; tr=maxtr; for (s=1;s 101;s++){ w=(1-tl)*th[s]; ls=(w-tr)/(2*w); if (ls 0)ls=0; c[s]=w*ls+tr; y[s]=th[s]*ls; } h=0.001; t=0; while(t 1000){ for (s=1;s 101;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s]; u1=ux(c1,l1); u[s][n1+1][n2+1]=u1; } } } for (s=1;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s+1]; u1=ux(c1,l1); ww[s][n1+1][n2+1]=u1; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ v[1][n1+1][n2+1][q+10]=-999; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ q=n2-n1; v[1][n1+1][n2+1][q+10]=u[1][n1+1][n2+1]; } } for (s=2;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ u1=u[s][n1+1][n2+1]; qx=q-n1+n2; pp=0; if (qx 10){ pp=100; } if (qx -10){ pp=100; } if (pp 50){ qx=0; } vs=-999; ns1=0; ns2=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[s-1][nx1+1][nx2+1][qx+10]; if (ww[s-1][nx1+1][nx2+1] u1){ v1=-999; } if (v1 vs){ ns1=nx1; } if (v1 vs){ ns2=nx2; } if (v1 vs){ vs=v1; } } } if (pp 50){ vs=-999; } gotoc[s][n1+1][n2+1][q+10]=ns1; gotoy[s][n1+1][n2+1][q+10]=ns2; gotoq[s][n1+1][n2+1][q+10]=qx; v[s][n1+1][n2+1][q+10]=vs; } } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ u1=u[100][n1+1][n2+1]; qx=n2-n1; vs=-999; ns1=0; ns2=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[99][nx1+1][nx2+1][qx+10]; if (ww[99][nx1+1][nx2+1] u1){ v1=-999; } if (v1 vs){ ns1=nx1; } if (v1 vs){ ns2=nx2; } if (v1 vs){ vs=v1; } } } endc[n1+1][n2+1]=ns1; endy[n1+1][n2+1]=ns2; endq[n1+1][n2+1]=qx; endv[n1+1][n2+1]=vs; } } maxv=-999; ns1=0; ns2=0; for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ if (endv[n1+1][n2+1] maxv){ ns1=n1; } if (endv[n1+1][n2+1] maxv){ ns2=n2; } if (endv[n1+1][n2+1] maxv){ maxv=endv[n1+1][n2+1]; } } } opc[100]=ns1; opy[100]=ns2; opc[99]=endc[opc[100]+1][opy[100]+1]; opy[99]=endy[opc[100]+1][opy[100]+1]; opq[99]=endq[opc[100]+1][opy[100]+1]; for (j=1;j 99;j++){ s=99-j; opc[s]=gotoc[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opy[s]=gotoy[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opq[s]=gotoq[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; } e=0; for (s=1;s 101;s++){ e=e+opc[s]*opc[s]+opy[s]*opy[s]; } for (s=1;s 101;s++){ c[s]=c[s]+opc[s]*h; y[s]=y[s]+opy[s]*h; } System.out.println(e); if (e 2){ h=h/2; } if (h 0.00001){ t=10000; } t=t+1; } for (s=1;s 100;s++){ dy=y[s+1]-y[s]; dc=c[s+1]-c[s]; if (dy 0)System.out.println(1-dc/dy); } } static double ux(double c1,double l1){ double up; double c3; double l3; int pp; pp=0; c3=c1; l3=l1; if (c3 0)pp=100; if (l3 0)pp=100; if (l3 1)pp=100; if (pp 50)c3=0.5; if (pp 50)l3=0.5; up=Math.log(c3)+Math.log(1-l3); if (pp 50)up=-999; return up; } static double bud(double tl,double tr,double th[]){ double bx; int s; double w1; double y1; double c1; double l1; bx=0; for (s=1;s 101;s++){ w1=(1-tl)*th[s]; l1=(w1-tr)/(2*w1); if (l1 0)l1=0; y1=th[s]*l1; c1=w1*l1+tr; bx=bx+y1-c1; } return bx; } static double seekwel(double tl,double tr,double th[]){ double sw; int s; double w1; double l1; double c1; sw=0; for (s=1;s 101;s++){ w1=(1-tl)*th[s]; l1=(w1-tr)/(2*w1); if (l1 0)l1=0; c1=w1*l1+tr; sw=sw+ux(c1,l1); } return sw; } }
https://w.atwiki.jp/ce00582/pages/1153.html
public class tax{ public static void main(String[] args) { double th1[]=new double[11]; double th2[]=new double[11]; double y1[][]=new double[11][11]; double y2[][]=new double[11][11]; double c[][]=new double[11][11]; double z[]=new double[11]; double u[][][][]=new double[11][3][3][3]; double w[][][][]=new double[10][3][3][3]; double rw[][][][]=new double[11][3][3][3]; double v[][][][][]=new double[11][3][3][3][11]; int goto1[][][][][]=new int[11][3][3][3][11]; int goto2[][][][][]=new int[11][3][3][3][11]; int gotoc[][][][][]=new int[11][3][3][3][11]; int gotoq[][][][][]=new int[11][3][3][3][11]; double endv[][][]=new double[3][3][3]; int endc[][][]=new int[3][3][3]; int end1[][][]=new int[3][3][3]; int end2[][][]=new int[3][3][3]; int endq[][][]=new int[3][3][3]; int opc[]=new int[11]; int op1[]=new int[11]; int op2[]=new int[11]; int opq[]=new int[11]; int s1,s2,j; double tk,tl,tr; double w1,w2,l1,l2; int i,t1,t2,t3,t4; double b1,b2,c1,we; double tr1,tr2,tr3; double maxwel,maxtk,maxtl,maxtr,maxu,maxv; int n1,n2,n3,q,qx; double u1,u2,v1; double h; double ww; int sig1,sig2,m1,m2,m3,mx1,mx2,mx3,nx1,nx2,nx3; int pp; int check; int sumep,ep; for (s1=1;s1 11;s1++){ th1[s1]=0.2*s1; } for (s2=1;s2 11;s2++){ th2[s2]=0.1*s2; } maxwel=-999; maxtk=0; maxtl=0; maxtr=0; for (n1=10;n1 40;n1++){ for (n2=5;n2 30;n2++){ tk=0.01*n1; tl=0.01*n2; tr=trs(tk,tl,th1,th2); we=wel(tk,tl,tr,th1,th2);; if (we maxwel)maxtk=tk; if (we maxwel)maxtl=tl; if (we maxwel)maxtr=tr; if (we maxwel)maxwel=we; } } tk=maxtk; tl=maxtl; tr=maxtr; for (s1=1;s1 11;s1++){ for (s2=1;s2 11;s2++){ w1=(1-tk)*th1[s1]; w2=(1-tl)*th2[s2]; l1=1-(w1+w2+tr)/(3*w1); l2=1-(w1+w2+tr)/(3*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=0; if (i==2)l2=(w2-tr)/(2*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=(w1-tr)/(2*w1); if (i==2)l2=0; y1[s1][s2]=th1[s1]*l1; y2[s1][s2]=th2[s2]*l2; c[s1][s2]=w1*l1+w2*l2+tr; } } t4=0; h=0.001; while(t4 100){ sumep=0; for (s2=1;s2 =10;s2++){ maxv=-999; t2=0; while(t2 100){ for (s1=1;s1 =10;s1++){ maxu=-999; for (sig1=-1;sig1 =1;sig1++){ for (sig2=-1;sig2 =1;sig2++){ m1=s1+sig1; m2=s2+sig2; if (m1 10)m1=10; if (m1 1)m1=1; if (m2 10)m2=10; if (m2 1)m2=1; l1=(y1[m1][m2])/th1[s1]; l2=(y2[m1][m2])/th2[s2]; c1=c[m1][m2]; u1=ux(c1,l1,l2); if (s2==m2)u1=-999; if (u1 maxu)maxu=u1; } } z[s1]=maxu; } for (s1=1;s1 =10;s1++){ for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ l1=(y1[s1][s2]+n1*h)/th1[s1]; l2=(y2[s1][s2]+n2*h)/th2[s2]; c1=c[s1][s2]+n3*h; u1=ux(c1,l1,l2); if (z[s1] u1)u1=-999; u[s1][n1+1][n2+1][n3+1]=u1; } } } } for (s1=1;s1 =10;s1++){ for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ check=0; for (sig1=-1;sig1 =1;sig1++){ for (sig2=-1;sig2 =1;sig2++){ m1=s1+sig1; m2=s2+sig2; if (m1 10)m1=10; if (m1 1)m1=1; if (m2 10)m2=10; if (m2 1)m2=1; l1=y1[m1][m2]/th1[m1]; l2=y2[m1][m2]/th2[m2]; c1=c[m1][m2]; u1=ux(c1,l1,l2); l1=(y1[s1][s2]+n1*h)/th1[m1]; l2=(y2[s1][s2]+n2*h)/th2[m2]; c1=c[s1][s2]+n3*h; u2=ux(c1,l1,l2); if (s2==m2)u2=-999; if (u2 u1)check=100; } } if(check 50)u[s1][n1+1][n2+1][n3+1]=-999; } } } } for (s1=1;s1 =9;s1++){ for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ l1=(y1[s1][s2]+n1*h)/th1[s1+1]; l2=(y2[s1][s2]+n2*h)/th2[s2]; c1=c[s1][s2]+n3*h; u1=ux(c1,l1,l2); w[s1][n1+1][n2+1][n3+1]=u1; } } } } for (s1=2;s1 =10;s1++){ for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ l1=(y1[s1][s2]+n1*h)/th1[s1-1]; l2=(y2[s1][s2]+n2*h)/th2[s2]; c1=c[s1][s2]+n3*h; u1=ux(c1,l1,l2); rw[s1][n1+1][n2+1][n3+1]=u1; } } } } for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ for (q=-5;q =5;q++){ v[1][n1+1][n2+1][n3+1][q+5]=-999; } } } } for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ q=n1+n2-n3; v[1][n1+1][n2+1][n3+1][q+5]=u[1][n1+1][n2+1][n3+1]; } } } for (s1=2;s1 =9;s1++){ for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ for (q=-5;q =5;q++){ u1=u[s1][n1+1][n2+1][n3+1]; qx=q-n1-n2+n3; pp=0; if (qx 5)pp=100; if (qx -5)pp=100; if (pp 10)qx=0; maxv=-999; mx1=0; mx2=0; mx3=0; for (m1=-1;m1 =1;m1++){ for (m2=-1;m2 =1;m2++){ for (m3=-1;m3 =1;m3++){ v1=u1+v[s1-1][m1+1][m2+1][m3+1][qx+5]; ww=w[s1-1][m1+1][m2+1][m3+1]; if (ww u1){ v1=-999; } if (rw[s1][n1+1][n2+1][n3+1] u[s1-1][m1+1][m2+1][m3+1])v1=-999; if (v1 maxv)mx1=m1; if (v1 maxv)mx2=m2; if (v1 maxv)mx3=m3; if (v1 maxv)maxv=v1; } } } if (pp 10)maxv=-999; v[s1][n1+1][n2+1][n3+1][q+5]=maxv; gotoc[s1][n1+1][n2+1][n3+1][q+5]=mx3; goto1[s1][n1+1][n2+1][n3+1][q+5]=mx1; goto2[s1][n1+1][n2+1][n3+1][q+5]=mx2; gotoq[s1][n1+1][n2+1][n3+1][q+5]=qx; } } } } } for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ u1=u[10][n1+1][n2+1][n3+1]; qx=n3-n1-n2; maxv=-999; mx1=0; mx2=0; mx3=0; for (m1=-1;m1 =1;m1++){ for (m2=-1;m2 =1;m2++){ for (m3=-1;m3 =1;m3++){ v1=u1+v[9][m1+1][m2+1][m3+1][qx+5]; ww=w[9][m1+1][m2+1][m3+1]; if (ww u1)v1=-999; if (rw[10][n1+1][n2+1][n3+1] u[9][m1+1][m2+1][m3+1])v1=-999; if (v1 maxv)mx1=m1; if (v1 maxv)mx2=m2; if (v1 maxv)mx3=m3; if (v1 maxv)maxv=v1; } } } endv[n1+1][n2+1][n3+1]=maxv; endc[n1+1][n2+1][n3+1]=mx3; end1[n1+1][n2+1][n3+1]=mx1; end2[n1+1][n2+1][n3+1]=mx2; endq[n1+1][n2+1][n3+1]=qx; } } } maxv=-999; nx1=0; nx2=0; nx3=0; for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ if (endv[n1+1][n2+1][n3+1] maxv)nx1=n1; if (endv[n1+1][n2+1][n3+1] maxv)nx2=n2; if (endv[n1+1][n2+1][n3+1] maxv)nx3=n3; if (endv[n1+1][n2+1][n3+1] maxv)maxv=endv[n1+1][n2+1][n3+1]; } } } opc[10]=nx3; op1[10]=nx1; op2[10]=nx2; opc[9]=endc[nx1+1][nx2+1][nx3+1]; op1[9]=end1[nx1+1][nx2+1][nx3+1]; op2[9]=end2[nx1+1][nx2+1][nx3+1]; opq[9]=endq[nx1+1][nx2+1][nx3+1]; for (j=1;j =8;j++){ s1=9-j; n1=op1[s1+1]; n2=op2[s1+1]; n3=opc[s1+1]; q=opq[s1+1]; opc[s1]=gotoc[s1+1][n1+1][n2+1][n3+1][q+5]; op1[s1]=goto1[s1+1][n1+1][n2+1][n3+1][q+5]; op2[s1]=goto2[s1+1][n1+1][n2+1][n3+1][q+5]; opq[s1]=gotoq[s1+1][n1+1][n2+1][n3+1][q+5]; } ep=0; for (s1=1;s1 =10;s1++){ ep=ep+opc[s1]*opc[s1]+op1[s1]*op1[s1]+op2[s1]*op2[s1]; } sumep=sumep+ep; for (s1=1;s1 =10;s1++){ c[s1][s2]=c[s1][s2]+opc[s1]*h; y1[s1][s2]=y1[s1][s2]+op1[s1]*h; y2[s1][s2]=y2[s1][s2]+op2[s1]*h; } if (ep 2)t2=1000; t2=t2+1; } System.out.println(maxv); } for (s1=1;s1 =10;s1++){ maxv=-999; t3=0; while(t3 100){ for (s2=1;s2 =10;s2++){ maxu=-999; for (sig1=-1;sig1 =1;sig1++){ for (sig2=-1;sig2 =1;sig2++){ m1=s1+sig1; m2=s2+sig2; if (m1 10)m1=10; if (m1 1)m1=1; if (m2 10)m2=10; if (m2 1)m2=1; l1=(y1[m1][m2])/th1[s1]; l2=(y2[m1][m2])/th2[s2]; c1=c[m1][m2]; u1=ux(c1,l1,l2); if (s1==m1)u1=-999; if (u1 maxu)maxu=u1; } } z[s2]=maxu; } for (s2=1;s2 =10;s2++){ for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ l1=(y1[s1][s2]+n1*h)/th1[s1]; l2=(y2[s1][s2]+n2*h)/th2[s2]; c1=c[s1][s2]+n3*h; u1=ux(c1,l1,l2); if (z[s2] u1)u1=-999; u[s2][n1+1][n2+1][n3+1]=u1; } } } } for (s2=1;s2 =10;s2++){ for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ check=0; for (sig1=-1;sig1 =1;sig1++){ for (sig2=-1;sig2 =1;sig2++){ m1=s1+sig1; m2=s2+sig2; if (m1 10)m1=10; if (m1 1)m1=1; if (m2 10)m2=10; if (m2 1)m2=1; l1=y1[m1][m2]/th1[m1]; l2=y2[m1][m2]/th2[m2]; c1=c[m1][m2]; u1=ux(c1,l1,l2); l1=(y1[s1][s2]+n1*h)/th1[m1]; l2=(y2[s1][s2]+n2*h)/th2[m2]; c1=c[s1][s2]+n3*h; u2=ux(c1,l1,l2); if (s1==m1)u2=-999; if (u2 u1)check=100; } } if(check 50)u[s2][n1+1][n2+1][n3+1]=-999; } } } } for (s2=1;s2 =9;s2++){ for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ l1=(y1[s1][s2]+n1*h)/th1[s1]; l2=(y2[s1][s2]+n2*h)/th2[s2+1]; c1=c[s1][s2]+n3*h; u1=ux(c1,l1,l2); w[s2][n1+1][n2+1][n3+1]=u1; } } } } for (s2=2;s2 =10;s2++){ for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ l1=(y1[s1][s2]+n1*h)/th1[s1]; l2=(y2[s1][s2]+n2*h)/th2[s2-1]; c1=c[s1][s2]+n3*h; u1=ux(c1,l1,l2); rw[s2][n1+1][n2+1][n3+1]=u1; } } } } for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ for (q=-5;q =5;q++){ v[1][n1+1][n2+1][n3+1][q+5]=-999; } } } } for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ q=n1+n2-n3; v[1][n1+1][n2+1][n3+1][q+5]=u[1][n1+1][n2+1][n3+1]; } } } for (s2=2;s2 =9;s2++){ for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ for (q=-5;q =5;q++){ u1=u[s2][n1+1][n2+1][n3+1]; qx=q-n1-n2+n3; pp=0; if (qx 5)pp=100; if (qx -5)pp=100; if (pp 10)qx=0; maxv=-999; mx1=0; mx2=0; mx3=0; for (m1=-1;m1 =1;m1++){ for (m2=-1;m2 =1;m2++){ for (m3=-1;m3 =1;m3++){ v1=u1+v[s2-1][m1+1][m2+1][m3+1][qx+5]; ww=w[s2-1][m1+1][m2+1][m3+1]; if (ww u1)v1=-999; if (rw[s2][n1+1][n2+1][n3+1] u[s2-1][m1+1][m2+1][m3+1])v1=-999; if (v1 maxv)mx1=m1; if (v1 maxv)mx2=m2; if (v1 maxv)mx3=m3; if (v1 maxv)maxv=v1; } } } if (pp 10)maxv=-999; v[s2][n1+1][n2+1][n3+1][q+5]=maxv; gotoc[s2][n1+1][n2+1][n3+1][q+5]=mx3; goto1[s2][n1+1][n2+1][n3+1][q+5]=mx1; goto2[s2][n1+1][n2+1][n3+1][q+5]=mx2; gotoq[s2][n1+1][n2+1][n3+1][q+5]=qx; } } } } } for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ u1=u[10][n1+1][n2+1][n3+1]; qx=n3-n1-n2; maxv=-999; mx1=0; mx2=0; mx3=0; for (m1=-1;m1 =1;m1++){ for (m2=-1;m2 =1;m2++){ for (m3=-1;m3 =1;m3++){ v1=u1+v[9][m1+1][m2+1][m3+1][qx+5]; ww=w[9][m1+1][m2+1][m3+1]; if (ww u1)v1=-999; if (rw[10][n1+1][n2+1][n3+1] u[9][m1+1][m2+1][m3+1])v1=-999; if (v1 maxv)mx1=m1; if (v1 maxv)mx2=m2; if (v1 maxv)mx3=m3; if (v1 maxv)maxv=v1; } } } endv[n1+1][n2+1][n3+1]=maxv; endc[n1+1][n2+1][n3+1]=mx3; end1[n1+1][n2+1][n3+1]=mx1; end2[n1+1][n2+1][n3+1]=mx2; endq[n1+1][n2+1][n3+1]=qx; } } } maxv=-999; nx1=0; nx2=0; nx3=0; for (n1=-1;n1 =1;n1++){ for (n2=-1;n2 =1;n2++){ for (n3=-1;n3 =1;n3++){ if (endv[n1+1][n2+1][n3+1] maxv)nx1=n1; if (endv[n1+1][n2+1][n3+1] maxv)nx2=n2; if (endv[n1+1][n2+1][n3+1] maxv)nx3=n3; if (endv[n1+1][n2+1][n3+1] maxv)maxv=endv[n1+1][n2+1][n3+1]; } } } opc[10]=nx3; op1[10]=nx1; op2[10]=nx2; opc[9]=endc[nx1+1][nx2+1][nx3+1]; op1[9]=end1[nx1+1][nx2+1][nx3+1]; op2[9]=end2[nx1+1][nx2+1][nx3+1]; opq[9]=endq[nx1+1][nx2+1][nx3+1]; for (j=1;j =8;j++){ s2=9-j; n1=op1[s2+1]; n2=op2[s2+1]; n3=opc[s2+1]; q=opq[s2+1]; opc[s2]=gotoc[s2+1][n1+1][n2+1][n3+1][q+5]; op1[s2]=goto1[s2+1][n1+1][n2+1][n3+1][q+5]; op2[s2]=goto2[s2+1][n1+1][n2+1][n3+1][q+5]; opq[s2]=gotoq[s2+1][n1+1][n2+1][n3+1][q+5]; } ep=0; for (s2=1;s2 =10;s2++){ ep=ep+opc[s2]*opc[s2]+op1[s2]*op1[s2]+op2[s2]*op2[s2]; } sumep=sumep+ep; for (s2=1;s2 =10;s2++){ c[s1][s2]=c[s1][s2]+opc[s2]*h; y1[s1][s2]=y1[s1][s2]+op1[s2]*h; y2[s1][s2]=y2[s1][s2]+op2[s2]*h; } if (ep 2)t3=1000; t3=t3+1; } System.out.println(maxv); } System.out.println(sumep); if (sumep 5)h=h/2; if (h 0.0001)t4=1000; t4=t4+1; } } static double bud(double tk,double tl,double tr,double th1[],double th2[]){ double b1,w1,w2,l1,l2,c1; int s1,s2,i; b1=0; for (s1=1;s1 11;s1++){ for (s2=1;s2 11;s2++){ w1=(1-tk)*th1[s1]; w2=(1-tl)*th2[s2]; l1=1-(w1+w2+tr)/(3*w1); l2=1-(w1+w2+tr)/(3*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=0; if (i==2)l2=(w2-tr)/(2*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=(w1-tr)/(2*w1); if (i==2)l2=0; c1=w1*l1+w2*l2+tr; b1=b1+th1[s1]*l1+th2[s2]*l2-c1; } } return b1; } static double wel(double tk,double tl,double tr,double th1[],double th2[]){ double we,w1,w2,l1,l2,c1; int s1,s2,i; we=0; for (s1=1;s1 11;s1++){ for (s2=1;s2 11;s2++){ w1=(1-tk)*th1[s1]; w2=(1-tl)*th2[s2]; l1=1-(w1+w2+tr)/(3*w1); l2=1-(w1+w2+tr)/(3*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=0; if (i==2)l2=(w2-tr)/(2*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=(w1-tr)/(2*w1); if (i==2)l2=0; c1=w1*l1+w2*l2+tr; we=we+Math.log(c1)+Math.log(1-l1)+Math.log(1-l2); } } return we; } static double trs(double tk,double tl,double[] th1,double[] th2){ double tr1,tr2,tr3,b1,b2; int t1; tr1=0.1; tr2=0.05; b1=bud(tk,tl,tr1,th1,th2); t1=0; while(t1 100){ b2=bud(tk,tl,tr2,th1,th2); tr3=tr2-b2*(tr2-tr1)/(b2-b1); b1=b2; tr1=tr2; tr2=tr3; if (b2*b2 0.001)t1=1000; t1=t1+1; } return tr2; } static double ux(double c1,double l1,double l2){ int pp; double u1,cp1,lp1,lp2; pp=0; cp1=c1; lp1=l1; lp2=l2; if (lp1 0)pp=100; if (lp1 0.99)pp=100; if (lp2 0)pp=100; if (lp2 0.99)pp=100; if (pp 50)lp1=0.5; if (pp 50)lp2=0.5; u1=Math.log(cp1)+Math.log(1-lp1)+Math.log(1-lp2); if (pp 50)u1=-999; return u1; } }
https://w.atwiki.jp/ce00582/pages/1094.html
class tax{ public static void main(String args[]){ double th[]= new double[101]; int s; double tl; double tr; double ls; double w; double cs; double ys; double b1; double b2; double tr1; double tr2; double tr3; double wel; double maxwel; double maxtl; int t; int n; for (s=1;s 101;s++){ th[s]=0.02*s; } maxtl=0; maxwel=-999; for (n=20;n 45;n++){ tl=0.01*n; tr1=0.01; tr2=0.02; tr=tr1; b1=bud(tl,tr1); t=0; while (t 100) { b2=bud(tl,tr2); tr3=tr2-b2*(tr2-tr1)/(b2-b1); tr1=tr2; tr2=tr3; b1=b2; if (b2*b2 0.00001){ t=1000; } t=t+1; } tr=tr2; wel=seekwel(tl,tr); if (wel maxwel){ maxtl=tl; } if (wel maxwel){ maxwel=wel; } } System.out.println(maxtl); } static double ux(double c1,double l1){ double u1; u1=Math.log(c1)+Math.log(1-l1); return u1; } static double bud(double tl,double tr){ double bx; int s; double w; double y1; double c1; double l1; double th[]=new double[101]; for (s=1;s 101;s++){ th[s]=0.02*s; } bx=0; for (s=1;s 101;s++){ w=(1-tl)*th[s]; l1=(w-tr)/(2*w); if (l1 0)l1=0; y1=th[s]*l1; c1=w*l1+tr; bx=bx+y1-c1; } return bx; } static double seekwel(double tl,double tr){ double sw; int s; double w; double l1; double c1; double th[]=new double[101]; for (s=1;s 101;s++){ th[s]=0.02*s; } sw=0; for (s=1;s 101;s++){ w=(1-tl)*th[s]; l1=(w-tr)/(2*w); if (l1 0)l1=0; c1=w*l1+tr; sw=sw+ux(c1,l1); } return sw; } }
https://w.atwiki.jp/relaxation/pages/11.html
コメントプラグイン @wikiのwikiモードでは #comment() と入力することでコメントフォームを簡単に作成することができます。 詳しくはこちらをご覧ください。 =>http //atwiki.jp/guide/17_60_ja.html たとえば、#comment() と入力すると以下のように表示されます。 名前 コメント
https://w.atwiki.jp/ce00582/pages/1098.html
class tax{ public static void main(String args[]){ double th[]= new double[101]; int s; double tl; double tr; double ls; double w1; double cs; double ys; double b1; double b2; double tr1; double tr2; double tr3; double wel; double maxwel; double maxtl; double maxtr; double maxv; int t; int n; double h; double c[]= new double[101]; double y[]= new double[101]; double u[][][]= new double[101][3][3]; double w[][][]= new double[100][3][3]; int n1; int n2; int gotoc[][][][]=new int[100][3][3][21]; int gotoy[][][][]=new int[100][3][3][21]; int gotoq[][][][]=new int[100][3][3][21]; double endv[][]=new double[3][3]; int endc[][]=new int[3][3]; int endy[][]=new int[3][3]; int endq[][]=new int[3][3]; double v[][][][]=new double[100][3][3][21]; int opc[]= new int[101]; int opy[]= new int[101]; int opq[]= new int[101]; int ns1; int ns2; int qx; int nx1; int nx2; double vs; double v1; double u1; double dc; double dy; int pp; int q; int e; int j; for (s=1;s 101;s++){ th[s]=0.02*s; } maxtl=0; maxwel=-999; maxtr=0; for (n=20;n 45;n++){ tl=0.01*n; tr1=0.01; tr2=0.02; tr=tr1; b1=bud(tl,tr1,th); t=0; while (t 100) { b2=bud(tl,tr2,th); tr3=tr2-b2*(tr2-tr1)/(b2-b1); tr1=tr2; tr2=tr3; b1=b2; if (b2*b2 0.00001)t=1000; t=t+1; } tr=tr2; wel=seekwel(tl,tr,th); if (wel maxwel)maxtl=tl; if (wel maxwel)maxtr=tr; if (wel maxwel)maxwel=wel; } tl=maxtl; tr=maxtr; for (s=1;s 101;s++){ w1=(1-tl)*th[s]; ls=(w1-tr)/(2*w1); if (ls 0)ls=0; c[s]=w1*ls+tr; y[s]=th[s]*ls; } t=0; h=0.001; while(t 1000){ for (s=1;s 101;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ cs=c[s]+n1*h; ls=(y[s]+n2*h)/th[s]; u[s][n1+1][n2+1]=ux(cs,ls); } } } for (s=1;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ cs=c[s]+n1*h; ls=(y[s]+n2*h)/th[s+1]; w[s][n1+1][n2+1]=ux(cs,ls); } } } for (s=2;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ u1=u[s][n1+1][n2+1]; qx=q-n1+n2; pp=0; if (qx 10)pp=100; if (qx -10)pp=100; if (pp 50)qx=0; vs=-999; ns1=0; ns2=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[s-1][nx1+1][nx2+1][qx+10]; if (w[s-1][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)ns1=nx1; if (v1 vs)ns2=nx2; if (v1 vs)vs=v1; } } if (pp 50)vs=-999; gotoc[s][n1+1][n2+1][q+10]=ns1; gotoy[s][n1+1][n2+1][q+10]=ns2; gotoq[s][n1+1][n2+1][q+10]=qx; v[s][n1+1][n2+1][q+10]=vs; } } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ u1=u[100][n1+1][n2+1]; qx=n2-n1; vs=-999; ns1=0; ns2=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[99][nx1+1][nx2+1][qx+10]; if (w[99][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)ns1=nx1; if (v1 vs)ns2=nx2; if (v1 vs)vs=v1; } } endc[n1+1][n2+1]=ns1; endy[n1+1][n2+1]=ns2; endq[n1+1][n2+1]=qx; endv[n1+1][n2+1]=vs; } } maxv=-999; ns1=0; ns2=0; for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ if (endv[n1+1][n2+1] maxv)ns1=n1; if (endv[n1+1][n2+1] maxv)ns2=n2; if (endv[n1+1][n2+1] maxv)maxv=endv[n1+1][n2+1]; } } opc[100]=ns1; opy[100]=ns2; opc[99]=endc[opc[100]+1][opy[100]+1]; opy[99]=endy[opc[100]+1][opy[100]+1]; opq[99]=endq[opc[100]+1][opy[100]+1]; for (j=1;j 99;j++){ s=99-j; opc[s]=gotoc[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opy[s]=gotoy[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opq[s]=gotoq[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; } e=0; for (s=1;s 101;s++){ e=e+opc[s]*opc[s]+opy[s]*opy[s]; } for (s=1;s 101;s++){ c[s]=c[s]+opc[s]*h; y[s]=y[s]+opy[s]*h; } if (e 1)h=h/2; if (h 0.0001)t=10000; System.out.print(t); System.out.print(","); System.out.println(maxv); t=t+1; } for (s=1;s 95;s++){ dc=c[s+1]-c[s]; dy=y[s+1]-y[s]; if (dy 0)System.out.println(1-dc/dy); } } static double ux(double c1,double l1){ double u1; double c3; double l3; int pp; pp=0; c3=c1; l3=l1; if (c3 0)pp=100; if (l3 0)pp=100; if (l3 1)pp=100; if (pp 50)c3=0.5; if (pp 50)l3=0.5; u1=Math.log(c3)+Math.log(1-l3); if (pp 50)u1=-999; return u1; } static double bud(double tl,double tr,double th[]){ double bx; int s; double w; double y1; double c1; double l1; bx=0; for (s=1;s 101;s++){ w=(1-tl)*th[s]; l1=(w-tr)/(2*w); if (l1 0)l1=0; y1=th[s]*l1; c1=w*l1+tr; bx=bx+y1-c1; } return bx; } static double seekwel(double tl,double tr,double th[]){ double sw; int s; double w; double l1; double c1; sw=0; for (s=1;s 101;s++){ w=(1-tl)*th[s]; l1=(w-tr)/(2*w); if (l1 0)l1=0; c1=w*l1+tr; sw=sw+ux(c1,l1); } return sw; } }
https://w.atwiki.jp/wasaxaza/pages/5.html
まとめサイト作成支援ツールについて @wikiにはまとめサイト作成を支援するツールがあります。 また、 #matome_list と入力することで、注目の掲示板が一覧表示されます。 利用例)#matome_listと入力すると下記のように表示されます #matome_list
https://w.atwiki.jp/wasaxaza/pages/12.html
人気商品一覧 @wikiのwikiモードでは #price_list(カテゴリ名) と入力することで、あるカテゴリの売れ筋商品のリストを表示することができます。 カテゴリには以下のキーワードがご利用できます。 キーワード 表示される内容 ps3 PlayStation3 ps2 PlayStation3 psp PSP wii Wii xbox XBOX nds Nintendo DS desctop-pc デスクトップパソコン note-pc ノートパソコン mp3player デジタルオーディオプレイヤー kaden 家電 aircon エアコン camera カメラ game-toy ゲーム・おもちゃ全般 all 指定無し 空白の場合はランダムな商品が表示されます。 ※このプラグインは価格比較サイト@PRICEのデータを利用しています。 たとえば、 #price_list(game-toy) と入力すると以下のように表示されます。 ゲーム・おもちゃ全般の売れ筋商品 #price_list ノートパソコンの売れ筋商品 #price_list 人気商品リスト #price_list