結果

問題 No.407 鴨等素数間隔列の数え上げ
コンテスト
ユーザー 小指が強い人
提出日時 2016-08-08 00:07:25
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
TLE  
実行時間 -
コード長 3,440 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,185 ms
コンパイル使用メモリ 167,044 KB
実行使用メモリ 10,020 KB
最終ジャッジ日時 2024-11-07 10:03:43
合計ジャッジ時間 6,592 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 5
other AC * 14 TLE * 2 -- * 15
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include "bits/stdc++.h"
#include <gmp.h>
using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef vector<int> veci;
typedef vector<ll> vecll;
typedef vector<string> vecs;
template<class T,class U> using Hash=unordered_map<T,U>;

#define REP(i, a, n) for(ll i = a; i < n; i++)
#define RREP(i, a, n) for(ll i = n-1; i >= a; i--)
#define rep(i, n) REP(i, 0, n)
#define rrep(i, n) RREP(i, 0, n)
#define MD 1000000007
#define _SPLIT " "

template<class T> T read(){T a;cin >> a;return a;}
template<class T> void read(T& a){cin >> a;}
template<class T,class ...Args> void read(T& a, Args&... args){cin >> a; read(args...);} 
template<class T> void rarr(T a, int n){for(int i = 0; i < n; i++) {cin >> a[i];}}
template<class T> void write(T a){cout << setprecision(12) << a << endl;}
template<class T,class ...Args> void write(T a, Args... args){cout << setprecision(12) << a << _SPLIT; write(args...);}
template<class T> void warr(vector<T> a, const char* c = " "){cout << a[0];for(int i = 1; i < (int)a.size(); i++)cout << c << a[i];cout << endl;;}
template<class T> void warr(T a, int n, const char* c = " "){cout << a[0];for(int i = 1; i < n; i++)cout << c << a[i];cout << endl;}
void split(string s, string delim, veci& result){result.clear();string::size_type pos = 0;while(pos != string::npos){string::size_type p = s.find(delim, pos);if(p == string::npos){result.push_back(atoi(s.substr(pos).data()));break;}else {result.push_back(atoi(s.substr(pos, p - pos).data()));}pos = p + delim.size();}}
void split(string s, string delim, vecs& result){result.clear();string::size_type pos = 0;while(pos != string::npos){string::size_type p = s.find(delim, pos);if(p == string::npos){result.push_back(s.substr(pos));break;}else {result.push_back(s.substr(pos, p - pos));}pos = p + delim.size();}}
ll gcd(ll a, ll b){while(true){ll k = a % b;if(k == 0)return b;a = b;b = k;}}
ll comb(ll n, ll m){ll p=1;m=min(m,n-m);for(ll i=1;i<=m;i++){p*=n-i+1;p/=i;}return p;}

ull power( const ull& a, ull k, const ull& n )
{
  if ( a == 0 || n == 0 ) return( 0 );
  if ( k == 0 ) return( 1 % n );

  ull currentMod = a % n;
  ull currentValue = ( ( k & 1 ) > 0 ) ? currentMod : 1;

  for ( k >>= 1 ; k > 0 ; k >>= 1 ) {
    currentMod = ( currentMod * currentMod ) % n;
    if ( ( k & 1 ) > 0 )
      currentValue = ( currentValue * currentMod ) % n;
  }

  return( currentValue );
}

bool MillerRabinTest( const ull p )
{
  if ( p == 1 ) return( false );
  if ( p == 2 ) return( true );
  if ( ( p & 1 ) == 0 ) return( false );

  ull q( p - 1 );
  ull k = 0;
  while ( ( q & 1 ) == 0 ) {
    q >>= 1;
    k += 1;
  }

  ull a( 2UL ); // a = 2
  ull gap( p );
  gap /= 100;
  if ( gap == 0 ) gap += 1;
  ull pm( p - 1 ); // p - 1

  while ( a < p ) {
    ull pow( power( a, q, p ) );
    a += gap;
    if ( pow == 1 ) continue;  // a^q _ 1 ( mod p )
    if ( pow == pm ) continue; // a^q _ -1 ( mod p )
    ull i( 1UL );
    while ( i < k ) {
      pow = power( pow, 2, p );
      if ( pow == pm ) break; // a^(2^i)q _ -1 ( mod p )
      i += 1;
    }
    if ( i == k ) return( false );
  }

  return( true );
}

int main(void)
{
    veci v;
    int n,l;
    cin>>n>>l;
    ll res=0;
    int i=2;
    while(true){
        if(!MillerRabinTest(i)){
            i++;
            continue;
        }
        int pm=i;
        int mn=(n-1)*pm;
        if(mn>l)break;
        res+=l-mn+1;
        i++;
    }
    write(res);
    return 0;
}
0