nyoj-素数环

直接枚举n!个效率不高,回溯法更优
#include <iostream> #include <cstdio> #include <cstring> using namespace std; const int Max = 10000; bool isp[Max]; int vis[Max]; int a[Max]= {1}; int n; void is_prime(int n) { memset(isp,true,sizeof(isp)); isp[0] = isp[1] = false; for(int i = 2; i <= 2*n; i++) if(isp[i] == true) //打点到2n的位置 { for(int j = i*2; j <= 2*n; j += i) isp[j] = false; } } void dfs(int cur) { if(cur == n && isp[a[0] + a[n-1]]) { for(int i = 0; i < n; i++) printf("%d",a[i]); printf("\n"); return; } for(int i = 2; i <= n; i++) if(!vis[i] && isp[i+a[cur-1]]) { a[cur] = i; vis[i] = 1; dfs(cur+1); vis[i] = 0; } } int main() { scanf("%d",&n); is_prime(n); dfs(1); return 0; }

 优化素数环

#include <iostream>
#include <cstring>

using namespace std;

const int Max = 40;
int visit[Max];
int a[Max];
int prime[Max];
int sign = 0;
int cnt = 0;

void prim()
{                            //数据是0~20,素数和最大18+19,不超过40,素数打表
     for(int i = 0; i <= Max;i++)
       prime[i] = 1; 
     prime[0] = prime[1] = 0;
     for(int  i = 2; i * i < Max; i++)
     if(prime[i] == 1)
       for(int  j = i + i; j <= Max; j = j+i)
          prime[j] = 0;
}  
void dfs(int x,int n,int step)
{    
     if(step == n&&prime[a[step-1]+a[0]])
     {
             sign = 1;
       for(int i = 0; i < n; i++)
         cout << a[i] << ' ';
         cout << endl;
       return;
       }
        
       for(int i = 2; i <= n;i++)
         if(/*prime[i]&&&*/prime[x+i]&&!visit[i]&&(x+i)%2==1)//剪枝,和前一个数的和不是素数的剪去
         {
            a[step] = i;
            visit[i]= 1;
            dfs(i,n,step+1);
            visit[i] = 0;
         }
            
}
int main()
{
  int n;
  while(cin >> n)
  {
            sign = 0;
            if(n == 0)
             return 0;
            cout << "Case " << ++cnt << ":" << endl;
                 if(n==1) //坑
                  {cout << 1 << endl;
                  continue;
                  }
             if(n%2==1)
               {       //奇数个数其中必有两个奇数相加为大于2的偶数,该偶数必定不是素数,无法组成素数环
                 cout << "No Answer" << endl;
                 continue;
} memset(visit,
0,sizeof(visit)); a[0] = 1; prim(); //打表 dfs(1,n,1); if(sign == 0) cout << "No Answer" << endl; } return 0; }

 

转载于:https://www.cnblogs.com/ekinzhang/p/4356684.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值