サタンプライム

彼らは誰なの？
彼らはしているPrimes含む666
これらサタン-素数です：[46663,266677,666599,666683,616669]
これらはありません。[462667,665669,36363631,555]

プロット

6661より大きいすべての数字の後ろにはサタンプリムスがいます

チャレンジ

整数n>6661が与えられた場合、サタンプライムを背後（または等しい）で、それ自体に最も近いものを見つけます。

例

Integerのn=30000背後には3つのSatan-Primes（SP）があります[6661, 16661, 26669]。
あなたのコード26669は、その後ろに最も近いものを返す必要があります

テストケース

入力->出力

6662->6661    
10000->6661    
66697->66697 (a SP returns himself)  
328765->326663  
678987->676661
969696->966677

ルール

Yorコードはn、あなたの言語の範囲内で動作するはずです。

これはcode-golfなので、バイト単位の最短回答が勝ちです！

Mathematica、82バイト

Last@Select[Prime@Range@PrimePi@#,!FreeQ[Subsequences[IntegerDigits@#],{6,6,6}]&]&

ニーム、9バイト

>ͻ:D+6S𝕚÷

説明：

>         Increment input
 ͻ        Start infinite loop
  :        Previous prime
   D       Duplicate
    +6     Push 666
      S    Swap
       𝕚   See if 666 is a substring of the top of the stack
        ÷  If true, break

オンラインでお試しください！

ゼリー、10 9バイト

@Dennisのおかげで10％節約！

ÆRwÐf666Ṫ

オンラインでお試しください！

説明

ÆR          # All primes in range [2, input]
   Ðf      # Keep those which satisfy
  w        # truthy if y is in x
     666   #           ^ (this is y)
        Ṫ  # Tail (take the last element)

Pyth、15 14バイト

Daveの助けを借りて1バイトを保存しました。

969696私のマシンのメモリエラーとそれ以上のものですが、十分なメモリが与えられていれば問題ありません。

ef&/`T*3\6P_TS

ここで試してみるか、テストスイートをご覧ください。

どうやって？

ef＆/ `T * 3 \ 6P_TSQ-完全なプログラム、最後に暗黙的な入力（Q）

             SQ-範囲[1、Q]
 f-フィルター。
          P_T-プライムですか？
  ＆-そして
   / `T * 3 \ 6-666が含まれています。
e-最後の要素。
                -結果を暗黙的に出力します。

Pyth、14バイト

ef/`T*\63fP_TS

ここで試してみてください！

05AB1E、11バイト

ƒNpN666å*iN

オンラインでお試しください！

Bash + Core Utils、51 49バイト

seq $1|tac|factor|awk 'NF==2&&/666/&&!a--&&$0=$2'

コマンドライン引数を取ります。数値が大きいと非常に遅くなる可能性があります。

Mathematica、64 62 61 53バイト

#//.i_/;!PrimeQ@i||ToString@i~StringFreeQ~"666":>i-1&

@KellyLowderのおかげで-1バイト

-8バイト（すごい）@Notatreeのおかげ

説明

入力してください。これらの条件下でデクリメントします：

• 入力が素数ではない、または

• 入力の数字に3つの6が連続して含まれていません。

サタンの素数に達するまで、それを繰り返します。

Perl 5、47バイト

46バイトのコード+ 1 -p

{$f=0|sqrt;1while$_%$f--;/666/*!$f||$_--*redo}

オンラインでお試しください！

ジャプト、14バイト

õ fj w æ_sø666

試して

そこにとして見ていました、50％の時間ベースのボーナス：を完了テストケース969696半秒下で。

説明

integerの暗黙的な入力U。

õ

から整数の配列を生成1しUます。

fj

フィルター（f）プライム。

w

逆。

æ_

以下1をチェックする関数を介して渡されたときに、真理値（この場合）を返す最初の要素を返します。

sø666

文字列に変換された整数（s）には（ø）666。

より高速な代替案、15バイト

繰り返しますが、元々は時間ベースのボーナスがあったので、これは代替の、はるかに高速なソリューションで、これ以上ゴルフをすることはできません。

U-@j *U´sø666}a

試して

PowerShell、128バイト

param($n)function f($a){for($i=2;$a-gt1){if(!($a%$i)){$i;$a/=$i}else{$i++}}}for(){if($n-match666-and($n-eq(f $n))){$n;exit}$n--}

オンラインでお試しください！

PowerShellには素因数分解が組み込まれていないため、これはPrime Factors Buddiesに関する私の答えからコードを借用しています。

入力を受け取り$n、入力function fの要素を計算する新しいものを宣言します$a。入力$aが素数の場合、これは単にを返します$a。

プログラムの主要部分は無限for()ループです。ループ内では、$n -matchesがesに反する666かどうか、および$n素数（つまり、の$nすべての要素に一致するかどうか）をチェックし$nます。そうである場合$n、パイプラインとexitに配置し、暗黙的な出力を使用します。それ以外の場合は、デクリメント$n--してループを継続します。

パイソン2、77の 76バイト

編集：@ Mr.Xcoderのおかげで-1バイト

実行時間が遅い、で実行 O(n^2)

lambda x:max(q for q in range(x+1)if"666"in`q`*all(q%t for t in range(2,q)))

オンラインでお試しください！

別の76バイトソリューション

lambda x:max(q*("666"in`q`*all(q%t for t in range(2,q)))for q in range(x+1))

オンラインでお試しください！

SymPy 73バイト

lambda x:max(q for q in primerange(0,x+1)if"666"in`q`)
from sympy import*

オンラインでお試しください！

PowerShell、71 69 64バイト

param($s)for(;$s-notmatch666-or(2..($s/2)|?{!($s%$_)});$s--){}$s

オンラインでお試しください！

• 私のマシンでは328765に30秒かかりますが、Tio.runの60秒の制限がタイムアウトします。

• 678987には約1.5分かかります。

• 969696には約4.5分かかります。

MATL、16バイト

ZqP"@V'666'Xf?@.

MATL Onlineでお試しください

説明

         Implicitly grab input (n)
Zq       Compute the primes up to n (output is in increasing order)
P        Flip the array (so larger primes come first)
"        For each prime
  @V     Convert it to a string
  '666'  Push the string literal '666' to the stack
  Xf     Find the location of '666' in the prime
  ?      If it was present...
    @.   Push it to the stack and break
         Implicitly display the stack contents

C ++ 389バイト

#include <iostream>
#include <boost/multiprecision/cpp_int.hpp>
#include <boost/multiprecision/miller_rabin.hpp>
using namespace boost::random;typedef boost::multiprecision::cpp_int Z;int main(int,char**v){mt19937 m(clock());independent_bits_engine<mt11213b,256,Z>g(m);Z n{v[1]},p;while(p++<=n)if(miller_rabin_test(p,25,g)&&p.convert_to<std::string>().find( "666" )!=-1)std::cout<<p<<" ";}

これは完全なプログラムです！
コンパイルするにはBoostが必要です。（または、お気に入りのオンラインC ++シェルにコピーして貼り付けます。）
コマンドラインからnを引数として実行します。

ゴルフをしていない：

#include <iostream>
#include <boost/multiprecision/cpp_int.hpp>
#include <boost/multiprecision/miller_rabin.hpp>
using namespace boost::random;

typedef boost::multiprecision::cpp_int integer;

int main( int argc, char** argv )
{
  mt19937 mt( clock() );
  independent_bits_engine <mt11213b, 256, integer> rng( mt );

  integer input {argv[ 1 ]};
  integer possible;

  while (possible++ <= input)
    if (
      // is_prime( possible )
      miller_rabin_test( possible, 25, rng )
    && 
      // possible has "666" in it
      (possible.convert_to <std::string> ().find( "666" ) != std::string::npos))

    std::cout << possible << " ";
}

ショートカットは、乱数テストの観点から作成されました。元のコードは6661で可能な素数のテストを開始し、2ずつ増加しました。また、nposの代わりに（-1）があるため、コンパイラの警告が表示されます。

それでも、これは非常に迅速に実行されます。私の古いAMD Sempron 130で、1,000,000未満の214のサタンプライムをすべて見つけるのに約40秒しかかかりませんでした。

：^ D

Bash + bsd-gamesパッケージ、33

• @FedericoPoloniのおかげで2バイト節約されました。

primes 2 $[$1+1]|grep 666|tail -1

オンラインでお試しください。

Pythonの3、85の83 80バイト

HalvardはPython 2で行われているため、4バイト短くなっています。

lambda k:max(x for x in range(k+1)if"666"in str(x)*all(x%i for i in range(2,x)))

オンラインでお試しください！

少し時間がかかりますが、そのO(n^2)複雑さのために非常に遅くなります。

JavaScript（ES6）、55 54バイト

@ThePirateBayのおかげで-1バイト。

f=n=>/666/.test(d=n)&eval("while(n%--d);d<2")?n:f(n-1)

入力が大きいと非常に遅くなります。このコードのゴルフの答えから適応した素数テスト。

タイミング

• 入力に10000は10秒かかりました
• 入力に328765は3分かかりました
• 入力に678987は9分かかりました
• 入力に969696は16分かかりました

テスト

これらの一部は、ブラウザを数分間ハングさせます。

f=n=>/666/.test(d=n)&eval("while(n%--d);d<2")?n:f(n-1)

function test(n) {
  O.value="Working..."
  setTimeout(_=>{
    let t = Date.now()
    O.value=`f(${n}) = ${f(n)} in ${(Date.now()-t)/1000}s`
  }, 10)
}

Tests<br>
<button onclick="test(6662)">6662</button>
<button onclick="test(10000)">10000</button>
<button onclick="test(328765)">328765</button>
<button onclick="test(678987)">678987</button>
<button onclick="test(969696)">6662</button><br>
Result<br>
<input id=O type=text size=25 readonly>

より高速なバージョン、56バイト

各テストケースを1秒以内に完了します。

f=n=>/666/.test(n)&&eval("for(d=2;n%d++;);d>n")?n:f(n-1)

;[6662, 10000, 328765, 678987, 969696].forEach(n=>console.log(`f(${n}) -> ${f(n)}`))

ルビー、67、66、58、56のバイト

の+7バイトを含む-rprime

->z{z.downto(1).find{|x|/666/=~x.to_s&&x.prime?}}

かなり高速で、最大~2^52約1秒2^64で5分以内に値を計算します（2011 MBP、Ruby 2.3.1）。

スタックス、10 バイト

ü>:Ñb/VP6─

実行してデバッグする

説明（開梱）：

^w:pc$666$#! Full program, implicit input-parsing
^            Increment input
 w           do-while:
  :p           Previous prime
    c$         Copy and stringify
      666$     Push "666"
          #    Number of occurences
           !   Logical not
             Implicit output

PHP、148バイト

<?php $p=[2];$s=[];for($i=3;$i<=$argv[1];$i++){foreach($p as $q)if($i%$q===0)continue 2;$p[]=$i;if(strpos($i,'666')!==false)$s[]=$i;}echo end($s);?>

オンラインでお試しください！

Perl 6、35バイト

my&f={/666/&&.is-prime??$_!!f $_-1}

オンラインでお試しください！

簡単な再帰的ソリューション。

C＃（.NET Core）、117 115 112バイト

f=>{for(int i=f;i>1;i--){int p=1,j=2;while(j<i)if(i%j++<1)p=0;if(p>0&$"{i}".Contains("666"))return i;}return 0;}

オンラインでお試しください！

• 不要なブレースを削除することで2バイト節約されました。
• int宣言を組み合わせることで3バイト節約されます。

これはもっと短くできると確信しています。たぶん、再帰的にfunc f外側のforループを呼び出して削除することによって。

再帰的アプローチ、85バイト

i=>{int p=1,j=2;while(j<i)if(i%j++<1)p=0;return p>0&$"{i}".Contains("666")?i:f(--i);}

オンラインでお試しください！

Func<int,int> f = null最初のセットを設定する必要があるため、このアプローチがcode-golfの範囲内に収まるかどうかはわかりfませんが、これは再度呼び出されますが、バイトにはカウントされません。明確化をいただければ幸いです。