POJ_2034_anti-prime_sequences


#include <iostream>
#include <cstring>
#include <cmath>
#include <cstdio>
#include <vector>
#include <algorithm>
#include <windows.h>
using namespace std;

#define Max 3300
#define Size 3200

bool isPrime[Max];

int n,m,d;


void primeRadiation(){
    
    memset( isPrime, true, sizeof( isPrime ) );
    
    isPrime[1] = false;
    
    for( int i = 4; i <= Size; i += 2 )
        isPrime[i] = false;
        
    for( int i = 3; i <= sqrt( Size ); ++i ){
        if( isPrime[i] == true ){
            
            const int step = 2 * i;
            
            for( int j = i * i; j <= Size; j += step ) 
                isPrime[j] = false;
        }
    }
}


bool isAntiPrime( const vector< int >& vec, const int& dd ){
    
    const int end = vec.size()-1;
    
    for( int i = 0; i <= ( end - dd + 1 ); ++i ){
        
        int val = 0;
        
        for( int j = i; j <= ( i + dd - 1 ); ++j )
            val += vec[j];
    
        if( isPrime[val] == true )
            return false;
            
    }
    
    return true;
    
}


int main(){
    
    primeRadiation();
    
    while( cin >> n >> m >> d ){
        
        DWORD begin,end;
        vector< int > vec;
        
        begin = GetTickCount();
    
        for( int i = n; i <= m; ++i )
            vec.push_back( i );
            
        bool b = false;
        
        do{
            
            int dd = d;
            
            while( true ){
                
                b = isAntiPrime( vec, dd );
                
                if( b == false ) 
                    break;
                    
                dd--;
                
                if( dd == 1 )
                    break;
                    
            }
            
            if( b == true ){
                
                for( int i = 0; i < vec.size(); ++i ) 
                    cout << vec[i] << " ";
                cout << endl;
                
                break;
                
            }
            
        }while( next_permutation( vec.begin(), vec.end() ) );
        
        end = GetTickCount();
        
        cout << ( end - begin ) << endl;
        
        if( b == false )
            cout << "No anti-prime sequence exists." << endl;
    }
    
    return 0;
    
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值