যদিও আমি এই সমস্যাটির একটি মন্টি কার্লো রেজোলিউশনের 6 দশমিক জায়গার যথাযথতা সহ whuber দ্বারা চিহ্নিত করাটির ব্যবহারিক অসম্ভবতা (বা অবাস্তবতা) সম্পর্কে আপত্তি করি না , তবে আমি মনে করব যে ছয় অঙ্কের নির্ভুলতার সাথে একটি রেজোলিউশন অর্জন করা যেতে পারে।
প্রথমত, নিম্নলিখিত Glen_b , কণা একটি নিশ্চল শাসন বিনিময়যোগ্য হয়, অত: পর তা যথেষ্ট (হিসাবে পর্যাপ্ততা ) বিভিন্ন কোষের দখল নিরীক্ষণ করতে, এই মারকভ প্রক্রিয়া পাশাপাশি গঠন। পরবর্তী সময়ে দখলগুলির বিতরণ পদক্ষেপ বর্তমান সময়ে টি দখলকারীদের দ্বারা নির্ধারিত সম্পন্ন হয় । ট্রান্সজিশন ম্যাট্রিক্স কে লেখা অবশ্যই অবাস্তব তবে এই রূপান্তরটির অনুকরণটি সোজা।t + 1টিকে
দ্বিতীয়ত, যেমন নজরে shabbychef , এক 450 বিজোড় (এমনকি) স্কোয়ারের দখল প্রক্রিয়া, যা যখন শুধুমাত্র এমনকি বার বিবেচনা করা বিজোড় স্কোয়ারের থাকে, অর্থাত অনুসরণ করতে পারেন ছক মার্কভ ম্যাট্রিক্স ।কে2
তৃতীয়ত, মূল সমস্যা শুধুমাত্র শূন্য ফ্রিকোয়েন্সি বিবেচনায় পি 0 পর 50 মার্কভ রূপান্তরের। প্রদত্ত শুরুর স্থান মার্কভ চেন নিশ্চল সম্ভাব্যতা বিতরণের জন্য একটি খুব উচ্চ মূল্য আছে ( এক্স ( T ) ) , এবং প্রদত্ত যে সব কোষ জুড়ে আপনাকে একটি গড়, ফোকাস পি 0 = 1পি^050( এক্স( টি ))আমরা বিবেচনা করতে পারি যেt=50সময়েশৃঙ্খলার(এক্স(টি))উপলব্ধি স্থির সম্ভাবনা বন্টন থেকে প্রাপ্তি। এই কম্পিউটিং খরচ একটি প্রধান হ্রাস এনেছে, হিসাবে আমরা এই নিশ্চল বন্টন থেকে সরাসরি অনুকরণ করতেπ, যা এমনকি কোণায় সম্ভাব্যতা সঙ্গে একটি MULTINOMIAL বন্টন সমানুপাতিক 2, 3, এবং 4, অন্যান্য কোষ কিনারায়, এবং ভেতরের কোষ যথাক্রমে
পি^0= 1450Σi = 1450আমি0( এক্স( 50 )আমি)
( এক্স( টি ))t = 50π
স্পষ্টতই, স্থিতিশীল বিতরণ
সমান খালি ঘরগুলির প্রত্যাশিত সংখ্যাকে 166.1069 এর সমান সরবরাহ করে ,
Σi = 1450( 1 - π)আমি)450
166.1069
pot=rep(c(rep(c(0,1),15),rep(c(1,0),15)),15)*c(2,
rep(3,28),2,rep(c(3,rep(4,28),3),28),2,rep(3,28),2)
pot=pot/sum(pot)
sum((1-pot)^450)-450
[1] 166.1069
যা মন্টি কার্লো অনুমানের খুব কাছাকাছি [10⁸ সিমুলেশনগুলির উপর ভিত্তি করে, যা আমার মেশিনে 14 ঘন্টা সময় নিয়েছিল]। তবে 6 দশমিকের জন্য পর্যাপ্ত নয়।166,11
হুবহু দ্বারা মন্তব্য করা হয়েছে , প্রশ্নের সঠিক উত্তর দেওয়ার জন্য অনুমানগুলি 2 দ্বারা গুণিত করা দরকার, সুতরাং চূড়ান্ত মান 332.2137,