আনুমানিক ব্রুনের কনস্ট্যান্ট


25

Brun ধ্রুবক মান যা বিপরীতকের এর সমষ্টি যমজ মৌলিক জোড়া ( 1/pএবং 1/(p+2)যেখানে pএবং p+2এগোয় উভয় প্রধানমন্ত্রী হয়)। এটা প্রায় 1.902160583104

একটি ইতিবাচক পূর্ণসংখ্যা দেওয়া হয়েছে N, যমজ দুটি প্রধান প্রাইম জুটির পারস্পরিক সংক্ষিপ্তসারকে সংক্ষিপ্ত করে ব্রুনের ধ্রুবক দেওয়া হয় যেখানে জুটির উভয় প্রাইম কম হয় Nএবং আনুমানিক আউটপুট দেয়।

বিধি

  • N আপনার ভাষার জন্য উপস্থাপনযোগ্য সীমার মধ্যে একটি ধনাত্মক পূর্ণসংখ্যা হবে।
  • ভাসমান-পয়েন্ট গাণিতিক ভুলগুলির কারণে কোনও সম্ভাব্য সমস্যা উপেক্ষা করে আপনার ভাষার ভাসমান-বিন্দু বাস্তবায়নের সীমাবদ্ধতার মধ্যে আউটপুটটি যথাযথ যথাযথ হওয়া উচিত। আপনার ভাষা যদি স্বেচ্ছাচারিতা-নির্ভুলতা গণিত করতে সক্ষম হয় তবে এটি অবশ্যই আইইইই 754 ডাবল-স্পষ্টতা পাটিগণিতের মতো কমপক্ষে যথাযথ হওয়া উচিত।
  • বিকল্পভাবে, একটি সঠিক ভগ্নাংশটি কোনও সামঞ্জস্যপূর্ণ, দ্ব্যর্থহীন বিন্যাসে আউটপুট হতে পারে।
  • যদি কোনও প্রধান একাধিক যমজ প্রধান জোড়া (উদাহরণস্বরূপ 5, উভয়ের অংশ (3, 5)এবং (5, 7)) প্রদর্শিত হয়, তবে এর পরস্পর প্রতিবার যোগফলকে অবদান রাখে।

পরীক্ষার মামলা

2 -> 0
6 -> 0.5333333333333333
10 -> 0.8761904761904762
13 -> 0.8761904761904762
100 -> 1.3309903657190867
620 -> 1.4999706034568274
100000 -> 1.67279958482774

সঠিক ভগ্নাংশ আউটপুট হতে পারে?
LegionMammal978

@ LegionMammal978 হ্যাঁ, আমি পরিষ্কার করব।
মেগো

পার্শ্ব দ্রষ্টব্য: ব্রুনের ধ্রুবকের জন্য মানটি 1.902160583104 ... কেবল অনুমান করা হয়; এমনকি প্রথম উল্লেখযোগ্য চিত্রও কঠোরভাবে গণনা করা হয়নি (এটি এটি 2 এর চেয়ে বড় বা কম কিনা তাও জানা যায়নি)।
গ্রেগ মার্টিন

@ গ্রেগমার্টিন এটি সত্য হওয়া সত্ত্বেও এটি বর্তমানে আমাদের সেরা অনুমানও।
মেগো

5 হ'ল একমাত্র প্রাইম যা দুটি প্রধান যুগল হিসাবে উপস্থিত হয়
খ্রিস্টান সিভারস 22'17

উত্তর:


25

পাইথন 3 , 78 77 75 70 68 62 বাইট

f=lambda n,k=3,m=1,j=0:k<n and-m%k*j*2/k+f(n,k+2,m*k**4,m%k/k)

@ Xnor কে 2 4 বাইট গলফ করার জন্য এবং আরও 4 টির পথ সুগম করার জন্য ধন্যবাদ!

এটি অনলাইন চেষ্টা করুন!

পটভূমি

উইলসনের উপপাদ্যটি স্মরণ করুন যে সমস্ত সংখ্যার জন্য কে> 1 ,

যেখানে a ≡ b (mod d) এর অর্থ হ'ল a - b সমানভাবে d দ্বারা বিভাজ্য , অর্থাত্ a এবং b এর সাথে অবশিষ্টাংশ থাকে যখন d দ্বারা ভাগ হয় ।

ইন ডাবল, অস্বাভাবিক-, উপ-এবং সুপার ফ্যাক্টরিয়ালগুলির জন্য উইলসন উপপাদ্য , লেখক ডবল ফ্যাক্টরিয়ালগুলির, যার উপর এই উত্তরটি তৈরী করে জন্য সরলীকরণ প্রমাণ করা। দ্বৈত গৌণিক একটি পূর্ণসংখ্যা এর ট ≥ 0 দ্বারা সংজ্ঞায়িত করা হয়

উপরোক্ত কাগজের 4 টি উপপাদ্যটি নীচে বর্ণিত।

উভয় পক্ষকে চতুর্থ শক্তিতে উন্নীত করে, আমরা এটি হ্রাস করি

সমস্ত বিজোড় primes জন্য পি । ২০০ Since সাল থেকে !! = 1 , সমতা পি = 2 এর জন্যও ধারণ করে ।

এখন, উইলসনের উপপাদ্যের সাথে একই কাজ করে তা প্রকাশ পায়

থেকে

এটা যে অনুসরণ করে

যখনই পি প্রধান হয়।

এখন, কে একটি বিজোড়, ধনাত্মক, সংমিশ্রণ পূর্ণসংখ্যার হতে দিন। সংজ্ঞা দ্বারা, সেখানে পূর্ণসংখ্যার অস্তিত্ব A, B> 1 যেমন যে ট = AB

যেহেতু কে বিজোড়, তাই এবং । সুতরাং, উভয় ক্রম 1, 3,…, কে - 2 এবং

যেখানে | বিভাজ্যতা নির্দেশ করে।

সমস্ত বিজোড় পূর্ণসংখ্যার জন্য সমষ্টি আপ <> 1

যেখানে P (ট) = 1 যদি মৌলিক এবং পি (ট) = 0 হলে যৌগিক হয়।

কিভাবে এটা কাজ করে

যখন ফাংশন f কে একক যুক্তি দিয়ে ডাকা হয়, k , m , এবং j 3 , 1 , এবং 0 হিসাবে আরম্ভ করা হয় ।

নোট করুন ((কে - 2) !!) 4 = 1 !! 4 = 1 = মি । আসলে, সমতা এম = ((কে - 2) !!) 4 সর্বদা ধরে রাখবে will একটি ভাসা এবং সর্বদা ((কে - 4)) 4 % (কে - 2) / (কে - 2) এর সমান হবে ।

যখন কে <এন , এর সঠিক যুক্তিটি andমূল্যায়ন করবে। যেহেতু ঞ = ((ট - 4) !!) 4 % (ট - 2) / (ট - 2) , প্রথম অনুচ্ছেদে প্রমাণিত হিসাবে, ঞ = 1 / (ট - 2) যদি ট - 2 প্রধানমন্ত্রী এবং ঞ = 0 না হলে। অনুরূপভাবে, যেহেতু মি% ট = ((ট - 2) !!) 4 সমান 1 যদি মৌলিক এবং 0 যদি না, -m% ট = ট - 1 যদি মৌলিক এবং -m% ট = 0 যদি না। অতএব, -m%k*j*2/kমূল্যায়ণ 2 (ট - 1) / (ট (ট - 2)) = ((ট - 2) + K) / (ট (ট - 2)) = 1 / ট + 1 টি / (ট - 2) যদি জোড় (কে - 2, কে)যমজ প্রাইম এবং 0 না থেকে থাকে if

উপরের গণনা করার পরে, আমরা ফলাফলটি পুনরাবৃত্তির কলটির রিটার্ন মানতে যুক্ত করি f(n,k+2,m*k**4,m%k/k) দ্বারা বৃদ্ধি পায় 2 তাই এটি শুধুমাত্র বিজোড় মান লাগে ‡ † , আমরা সংখ্যাবৃদ্ধি আছি দ্বারা 4 থেকে এমকে 4 = ((ট - 2) !!) 44 = (ট !!) 4 , এবং বর্তমান মান পাস মি% K / ট - যা সমান 1 / ট যদি 'পুরানো' একটি মৌলিক এবং 0 যদি না - প্যারামিটার হিসাবে ফাংশন কল করা হয়।

অবশেষে, একবার সমান বা তার চেয়ে অনেক বেশী এন , ফিরে আসবে মিথ্যা এবং পুনরাবৃত্তির স্টপ। চ (এন) এর রিটার্ন মান হ'ল সমস্ত 1 / কে + 1 / (কে - 2) এর সমষ্টি হবে ( যেমন - কে - 2, কে) হ'ল একটি যুগল প্রধান যুগল এবং কে <এন , পছন্দসই হিসাবে।


Back পটভূমি অনুচ্ছেদে ফলাফলগুলি কেবল বিজোড় পূর্ণসংখ্যার জন্য ধারণ করে। যেহেতু পূর্ণসংখ্যাগুলিও দু'জন প্রাইম হতে পারে না, তাই আমরা সেগুলি নিরাপদে এড়িয়ে যেতে পারি।


আমি মনে করি আপনার অভিব্যক্তি হিসাবে একই m%k*(j/k+j/(k-2))
xnor

হ্যাঁ, এটি কাজ করে। ধন্যবাদ!
ডেনিস


বিজোড় জন্য ((k-2)!!)^4 = p(k)মডুলো চমৎকার পর্যবেক্ষণ । আমি আপনার যুক্তি দিয়ে কাজ করি নি, তবে এখানে একটি উপস্থিত হয়েছি (এটি সংক্ষেপে একই হতে পারে)। সেটে কাজের মডুলো , সন্ধ্যাগুলি হ'ল প্রতিকূলতার নেতিবাচক। তাই ,। উইলসনের উপপাদ্য আমাদের তা বলে । যেহেতু , আমরা এবং তাই । ppp{1,2,..,p-1}prod(odds) = ± prod(evens)prod(all) = - p(k)prod(all) = prod(odds) * prod(evens) = prod(odds) * ± prod(evens)prod(odds)^2 = ±p(k)prod(odds)^4 = p(k)^2 = p(k)
xnor

নিস! আমি যোগফলটি একটি একক ভগ্নাংশ হিসাবে প্রকাশ করার চেষ্টা করেছি, তবে জেতে এর অংশীকরণের অংশটি আমার কাছে ঘটেনি। আবার ধন্যবাদ! আপনার প্রমাণটি কাগজ থেকে প্রাপ্ত একটি তুলনায় অনেক সহজ।
ডেনিস

7

জেলি , 15 14 বাইট

’ÆRµ_2fµ+2;µİS

এটি অনলাইন চেষ্টা করুন!

কিভাবে এটা কাজ করে

’ÆRµ_2fµ+2;µİS  Main link. Argument: n

’               Decrement; yield n-1.
 ÆR             Prime range; yield all primes in [1, ..., n-1].
   µ            New chain. Argument: r (prime range)
    _2          Subtract 2 from all primes.
      f         Filter; keep all p-2 that appear in r.
       µ        New chain. Argument: t (filtered range)
        +2      Add 2 to all primes in s.
          ;     Concatenate with s.
           µ    New chain. Argument: t (twin primes)
            İ   Take the inverses.
             S  Sum.

5

জেলি , 14 14 বাইট (@ ডেনিসের সামান্য সহায়তায়)

’ÆRṡ2_/2+$$ÐḟFİS

এটি অনলাইন চেষ্টা করুন!

আমার পূর্ববর্তী উত্তরটি উন্নত করার চেষ্টা করার সময়, আমি একটি সম্পূর্ণ ভিন্ন অ্যালগরিদম চিন্তা করেছিলাম এবং এটি কিছুটা সংক্ষিপ্ত আকারে আসে। আমি এটির জন্য একটি আলাদা পোস্ট ব্যবহার করছি, যেমন একটি উত্তরের জন্য এখানে একটি স্ট্যান্ডার্ড যা আলাদা কৌশল ব্যবহার করে।

ডেনিস প্রতিস্থাপন পরামর্শ _/2+$$Ðḟদিয়ে Iċ¥Ðf2; আমি একটি ডায়াডিক ফিল্টার সম্ভাবনা সম্পর্কে সম্পূর্ণরূপে ভুলে গিয়েছিলাম। এই হিসাবে, এই অ্যালগরিদম এখন ডেনিসের উত্তর ব্যবহৃত একটির সাথে সম্পর্কযুক্ত ties

ব্যাখ্যা

’ÆRṡ2Iċ¥Ðf2FİS
’                  Decrement.
 ÆR                Primes from 2 to the argument inclusive
                   (i.e. 2 to the original input exclusive).
   ṡ2              Take overlapping slices of size 2.
        Ðf         Keep only elements where the following is true:
       ¥           {the second parse of, which parses like this}
     Iċ   2          the differences (I) contain (ċ) 2
           F       Flatten.
            İ      Take 1/x {for every list element}.
             S     Sum.

2_/2+$$Ðḟহয়ে যেতে পারে Iċ¥Ðf2
ডেনিস

4

ব্র্যাচল্যাগ , 17 বাইট

{>I-₂:I{ṗ/₁}ᵐ}ᶠc+

এটি অনলাইন চেষ্টা করুন!

এটি চকচকে কোড পৃষ্ঠা সহ ব্র্যাচল্যাগের একেবারে নতুন সংস্করণ!

ব্যাখ্যা

{            }ᶠ        Find all valid outputs of the predicate in brackets
               c+      Output is the sum of that list after flattening it

 >I                    Input > I
   -₂:I                The list [I-2, I]
       {   }ᵐ          Map:
        ṗ/₁              Must be prime and the output is its inverse

3

এমএটিএল , 16 বাইট

liqZqtd2=)t2+h/s

এটি অনলাইন চেষ্টা করুন!

13উদাহরণ হিসাবে ইনপুট বিবেচনা করুন ।

l     % Push 1
      %   STACK: 1
i     % Input N
      %   STACK: 1, 13
q     % Subtract 1
      %   STACK: 1, 12
Zq    % Primes up to that
      %   STACK: 1, [2 3 5 7 11]
t     % Duplicate
      %   STACK: 1, [2 3 5 7 11], [2 3 5 7 11]
d     % Consecutive differences
      %   STACK: 1, [2 3 5 7 11], [1 2 2 4]
2=    % Compare with 2, element-wise
      %   STACK: 1, [2 3 5 7 11], [0 1 1 0]
)     % Use as logical index to select elements from array
      %   STACK: 1, [3 5]
t     % Duplicate
      %   STACK: 1, [3 5], [3 5]
2+    % Add 2, element-wise
      %   STACK: 1, [3 5], [5 7]
h     % Concatenate horizontally
      %   STACK: 1, [3 5 5 7]
/     % Divide, element-wise
      %   STACK: [0.3333 0.2 0.2 0.1429]
s     % Sum of array. Implicitly display
      %   STACK: 0.8762

2

গণিত, 48 47 বাইট

জংহওয়ান মিনকে 1 বাইট বাঁচানোর জন্য ধন্যবাদ!

If[PrimeQ/@(i&&(g=i-2)),1/i+1/g,0]~Sum~{i,#-1}&

নামহীন ফাংশন ইনপুট হিসাবে ইতিবাচক পূর্ণসংখ্যার গ্রহণ এবং একটি সঠিক ভগ্নাংশ ফিরে; উদাহরণস্বরূপ, If[PrimeQ/@(i&&(g=i-2)),1/i+1/g,0]~Sum~{i,#-1}&[10]ফেরৎ 92/105

If[PrimeQ/@(i&&(g=i-2)),1/i+1/g,0]উভয় iএবং i-2প্রধান কিনা তা পরীক্ষা করে , যদি তাদের 0না হয় এবং যদি না হয় তবে তাদের পারস্পরিক সংখ্যার যোগফল ফিরিয়ে দেয় । ~Sum~{i,#-1}&তারপরে ইনপুটটির iচেয়ে কম মানের সমস্ত মানের জন্য এই অবদানগুলির যোগফলটি প্রদান করে ।

পূর্ববর্তী জমা:

If[And@@PrimeQ@{i,g=i-2},1/i+1/g,0]~Sum~{i,#-1}&

এখন যে ঠিক স্পকি। আমি হাল ছেড়ে দিলাম। ⚐
LegionMammal978

আমি ভাবলাম "সঠিক ভগ্নাংশ" মানে ম্যাথমেটিকা ​​:)
গ্রেগ মার্টিন

-1 বাইট:If[PrimeQ/@(i&&(g=i-2)),1/i+1/g,0]~Sum~{i,#-1}&
জংহওয়ান মিন

N@কোডের সামনে দুটি বাইট যুক্ত করে কেউ একটি স্বেচ্ছাসেবী-নির্ভুলতা নম্বর পেতে পারে ।
জংহওয়ান মিন

কন্ডিশনে চমৎকার গল্ফিং! এটি সত্য যে Nএকটি বাস্তব সংখ্যায় দশমিকের প্রায় প্রত্যাবর্তন করে; তবে এটির জন্য 6 টিরও বেশি সিগ ডুমুর প্রদর্শন করতে অতিরিক্ত বাইটের প্রয়োজন হয় এবং কতগুলি সিগ ডুমুর প্রদর্শিত হয় তা এখনও ভগ্নাংশের চেয়ে কম সঠিক।
গ্রেগ মার্টিন

2

অক্টাভা, 45 বাইট

@(n)sum(all(isprime(m=[h=3:n-1;h-2]))*m'.^-1)

ব্যাখ্যা:

m=[h=3:n-1;h-2]             generate an concatenate two ranges 3:n-1 and 1:n-3
rec=m'.^-1                  transpose and reciprocal
idx=all(isprime(m))         create a logical [0 1 ..] array  if both ranges are prime set 1 else set 0
sum1 = idx * rec            matrix multiplication(extrat elements with logical index and sum along the first dimension)
sum(sum1)                   sum along the second dimension  

এটি অনলাইন চেষ্টা করুন!


2

জাভাস্ক্রিপ্ট (ES6), 67 66 বাইট

@ আরনাউল্ডকে ধন্যবাদ 1 বাইট সংরক্ষিত

f=n=>--n>1&&((p=x=>n%--x?p(x):x==1)(n)&&p(n-=2)&&1/n+++1/++n)+f(n)

falseপরীক্ষার ক্ষেত্রে আউটপুট 2, যা ডিফল্টরূপে অনুমোদিত

পরীক্ষার স্নিপেট


আমি মনে করি 1/n+++1/++nএকটি বাইট সংরক্ষণ করে
আর্নল্ড

@ আরনাউল্ড ধন্যবাদ কোনও কারণে আমি জানতাম না যে সর্বদা ত্রুটি ছুঁড়ে +++না ...
ETH প্রোডাকশনগুলি


1

জেলি , 19 বাইট

’ÆRḊµ_Æp=2Tịµ_2;µİS

এটি অনলাইন চেষ্টা করুন!

আমার অনুভূতি আছে যে এটি অসম্পূর্ণ, তবে আমি কীভাবে তা অবিলম্বে দেখতে পারছি না।

ব্যাখ্যা

’ÆRḊµ_Æp=2Tịµ_2;µİS
 ÆR                  Generate all primes from 2 to n inclusive
’                    Subtract 1
   Ḋ                 Remove first element
’ÆRḊ                 Generate all primes from 3 to n-1 exclusive

     _Æp             Subtract the previous prime (i.e. calculate the prime gap)
        =2           Compare to 2
          Tị         Take elements of the input where the comparison is true
     _Æp=2Tị         Filter a list of primes to the latter halves of prime pairs

             _2      Subtract 2
               ;     Append
             _2;     Append the list to the list with 2 subtracted from it
                 İ   Take reciprocals
                  S  Sum
                 İS  Take the sum of the reciprocals

এই µসমস্ত অংশকে পাইপলাইন-স্টাইলে একসাথে সংযুক্ত করুন, যার প্রতিটি তার আউটপুট হিসাবে আগে একের আউটপুট গ্রহণ করে।



1

পার্ল 6 , 59 51 বাইট

{sum 1 «/»grep((*-(2&0)).is-prime,^$_).flatmap:{$_-2,$_}}

{sum 1 «/»grep(*.all.is-prime,(-2..*Z ^$_)).flat}

-2..* Z ^$_তালিকার -2, -1, 0, 1, ...সাথে তালিকার অসীম তালিকাটি জিপ করে 0, 1, ... $_-1( $_ফাংশনটির পক্ষে যুক্তি হওয়ায়), তালিকাটি তৈরি করে(-2, 0), (-1, 1), (0, 2), ..., ($_-3, $_-1) । (স্পষ্টতই 3 এর চেয়ে কম সংখ্যক কোনওটিই প্রাইম জুটিতে থাকতে পারে না, তবে 3..* Z 5..^$_এটি কয়েক বাইট দীর্ঘ হয় এবং অতিরিক্ত সংখ্যার কোনওটিই প্রধান নয়))

দ্য grepনির্বাচন কেবলমাত্র সেই জোড়া যেখানে সমস্ত (যে, উভয়) সংখ্যার মৌলিক, আর flatচ্যাপ্টা তাদের সংখ্যার একটি প্লেইন তালিকায়।

«/»ডিভিশন হাইপারোপরেটর; ডান এবং 1বামে তালিকাটি সহ, এটি প্রধান জোড়াগুলির তালিকাটিকে তাদের পারস্পরিক ক্রিয়াকলাপে রূপান্তরিত করে, যার পরে সংক্ষিপ্তসার করা হয় sum


1

ক্লোজার, 147 বাইট

(fn[n](let[p #(if(> % 2)(<(.indexOf(for[a(range 2 %)](mod % a))0)0))](reduce +(for[a(range 2 n)](if(and(p a)(p(- a 2)))(+(/ 1 a)(/ 1(- a 2)))0)))))

এবং ক্লোজিউর সর্বশেষ হিসাবে সর্বশেষে মারা যায়।

Ungolfed:

; Returns the primality of a number.
(defn prime? [n]
  (if (> n 2)
    (< (.indexOf (for [a (range 2 n)] (mod n a)) 0) 0)))

; Calculates the actual Brun's Constant. ' (Stupid highlighter)
(defn brunsconst [n]
  ; Adds all of the entries together
  (reduce
    +
    ; For a in range(2, n):
    (for [a (range 2 n)]
      (let [b (- a 2)]
        ; If both a and a-2 are prime:
        (if (and (prime? a) (prime? b))
          ; Place (1/a + 1/a-2) on the array, else 0
          (+ (/ 1 a) (/ 1 b)) 0)))))


0

বাশ + জিএনইউ ইউটিলিটিস, 86 85 বাইট

for((k=4;k<$1;k++,j=k-2)){ [ `factor $k $j|wc -w` = 4 ]&&x=$x+1/$k+1/$j;};bc -l<<<0$x

এটি অনলাইন চেষ্টা করুন!

একটি বৃহত পাটিগণিত প্রকাশ করে এবং তারপরে এটি bc -lমূল্যায়নের জন্য ফিড দেয়।

সম্পাদনা করুন: নেস্টেড কমান্ড প্রতিস্থাপনের সাথে পুরানো সংস্করণ থেকে ভুল করে একটি $ (...) জুটিতে রেখে দেওয়া হয়েছে; বাইট সংরক্ষণ করতে ব্যাকটিক্সে পরিবর্তিত হয়েছে।


0

এপিএল NARS, 216 বাইট, 108 চর

  r←z n;h;i;k;v
  i←0⋄n-←1⋄h←1+⍳n-1⋄→B
A:k←i⊃h⋄h←k∪(0≠k∣h)/h
B:→A×⍳(⍴h)≥i+←1
  r←+/÷(v-2),v←(h=1⌽h+2)/h

এটি "ক্রাইভেলো দি ইরাতোস্টিন" ব্যবহার করবে 1. অনুরোধের প্রাইমগুলির 1..arg এ সাবলিস্টটি সন্ধান করার জন্য। টেস্ট:

  z¨2 6 10 13 100 620
0 0.5333333333 0.8761904762 0.8761904762 1.330990366 1.499970603 
  z 100000
1.672799585
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.