题目47：找出最小的具有四个不同质数因子的四个连续整数

欧拉计划

Distinct primes factors

The first two consecutive numbers to have two distinct prime factors are:

14 = 2 × 7
15 = 3 × 5

The first three consecutive numbers to have three distinct prime factors are:



Find the first four consecutive integers to have four distinct prime factors. What is the first of these numbers?

题目：

最小的两个具有两个不同质数因子的连续整数是：

14 = 2 × 7
15 = 3 × 5

最小的三个具有三个不同质数因子的连续整数是：



找出最小的四个具有四个不同质数因子的整数。它们之中的第一个是多少？

迷雾少年

怎么我出来了好多。。
16:2 2 2 2
24:2 2 2 3
36:2 2 3 3
40:2 2 2 5
54:2 3 3 3
56:2 2 2 7
60:2 2 3 5
81:3 3 3 3
84:2 2 3 7
88:2 2 2 11
90:2 3 3 5
100:2 2 5 5
104:2 2 2 13
126:2 3 3 7
..........

1. 
   
2. #include "myLib\\myLib.h"
   
3. #include <vector>
   
4. #pragma  comment(lib,"myLib\\myLib.lib")
   
5. 
   
6. //分解质因数并且那啥
   
7. bool div(int n)
   
8. {
   
9.         vector<int> result;
   
10.         int num = 0;
    
11.         int t = n;
    
12.         while (1!=n)
    
13.         {
    
14.                 for (int i =2;i<=n;i++)
    
15.                 {
    
16.                         if(n%i==0)
    
17.                         {
    
18.                                 //i就是其中一个因数
    
19.                                 //如果有一个因数不是质数 false
    
20.                                 if(!iszhishu(i))
    
21.                                         return false;
    
22.                                 result.push_back(i);
    
23.                                 num++;
    
24. 
    
25.                                 if(num>4)//超过4个质因数了 咋办
    
26.                                         return false;
    
27. 
    
28.                                 n/=i;
    
29.                                 i=1;
    
30. 
    
31.                         }
    
32.                 }
    
33.         }
    
34.         if(num!=4) return false;
    
35. 
    
36.         std::cout<<t<<":";
    
37.         for (int k=0;k<result.size();k++)
    
38.         {
    
39.                 std::cout<<result[k]<<" ";
    
40.         }
    
41.         std::cout<<endl;
    
42.         return true;
    
43. 
    
44. }
    
45. int main(void)
    
46. {
    
47.        
    
48.     for (int i = 1;i<10000;i++)
    
49.     {
    
50.                 div(i);
    
51.                   
    
52.     }
    
53.         return 0;
    
54. }

复制代码

迷雾少年

iszhishu 是个判断是否质数的函数 没带上

jerryxjr1220

迷雾少年 发表于 2016-8-30 21:56
怎么我出来了好多。。
16:2 2 2 2
24:2 2 2 3

看错题目，是4个连续数字，不是3个。
答案应该是134043

1. 
   
2. class henum:
   
3.     def __init__(self,num):
   
4.         henum.num = 0
   
5.         henum.yue = []
   
6. 
   
7. def yueshu(num):
   
8.     yuelist = []
   
9.     for i in range(2,int(num**0.5+1)):
   
10.         if isPrime(i):
    
11.             if num%i == 0:
    
12.                 yuelist.append(i)
    
13.                 num1 = int(num/i)
    
14.                 if isPrime(num1):
    
15.                     yuelist.append(num1)
    
16.                 else:
    
17.                     yuelist = yuelist + yueshu(num1)
    
18.     return yuelist
    
19. 
    
20. def isPrime(num):
    
21.     flag = 1
    
22.     if num == 2:
    
23.         return flag
    
24.     elif num % 2 == 0:
    
25.         flag = 0
    
26.         return flag
    
27.     else:
    
28.         for i in range(3,int(num**0.5+1),2):
    
29.             if num % i == 0:
    
30.                 flag = 0
    
31.                 break
    
32.         return flag
    
33. 
    
34. hnumlist = []
    
35. for j in range(100000,200000):
    
36.     if isPrime(j) == 0:
    
37.         h=henum(j)
    
38.         h.num = j
    
39.         h.yue = list(set(yueshu(j)))
    
40.         hnumlist.append(h)
    
41. 
    
42. for k in range(3,len(hnumlist)):
    
43.     if hnumlist[k-3].num == hnumlist[k-2].num -1 and hnumlist[k-2].num == hnumlist[k-1].num -1 and hnumlist[k-1].num == hnumlist[k].num -1:
    
44.         if len(hnumlist[k-3].yue) == 4 and len(hnumlist[k-2].yue) == 4 and len(hnumlist[k-1].yue) == 4 and len(hnumlist[k].yue) == 4:
    
45.                 print(hnumlist[k-3].num, hnumlist[k-3].yue)
    
46.                 break
    

复制代码

芒果加黄桃

1. # encoding:utf-8
   
2. # 找出最小的不能写作一个质数与一个平方数的二倍之和的奇合数
   
3. from time import time
   
4. from math import sqrt
   
5. def getPrimes(N):
   
6.     primes = [True] * N
   
7.     primes[0], primes[1] = False, False
   
8.     for i, prime in enumerate(primes):
   
9.         if prime:
   
10.             for j in range(i * i, N, i):
    
11.                 primes[j] = False
    
12.     return [k for k, p in enumerate(primes) if p]
    
13. def findPrimeFactors(n, l_pr, s_pr):
    
14.     if n in s_pr:
    
15.         return [n]
    
16.     sqr_n = int(sqrt(n)) + 1
    
17.     result = list()
    
18.     for pr in l_pr:
    
19.         if n % pr == 0:
    
20.             result.append(pr)
    
21.             result.extend(findPrimeFactors(int(n / pr), l_pr, s_pr))
    
22.             break
    
23.         if pr > sqr_n:
    
24.             break
    
25.     return result
    
26. def euler047(N=150000):
    
27.     l_primes = getPrimes(N)
    
28.     s_primes = set(l_primes)
    
29.     l_result = dict()
    
30.     for each in range(647, N):
    
31.         tmp = findPrimeFactors(each, l_primes, s_primes)
    
32.         if len(set(tmp)) == 4:
    
33.             l_result[each] = tmp
    
34.         else:
    
35.             l_result.clear()
    
36.         if len(l_result) == 4:
    
37.             print(l_result)
    
38.             break
    
39. if __name__ == '__main__':
    
40.     start = time()
    
41.     euler047()
    
42.     print('cost %.6f sec' % (time() - start))
    

复制代码

{134043: [3, 7, 13, 491], 134044: [2, 2, 23, 31, 47], 134045: [5, 17, 19, 83], 134046: [2, 3, 3, 11, 677]}
cost 0.955096 sec

愤怒的大头菇

1. import math
   
2. 
   
3. 
   
4. def Prime(x):
   
5.     if x > 1:
   
6.         if x == 2:
   
7.             return True
   
8.         if x % 2 == 0:
   
9.             return False
   
10.         for i in range(3, int(math.sqrt(x)+1), 2):
    
11.             if x % i == 0:
    
12.                 return False
    
13.         return True
    
14. 
    
15. 
    
16. def cycle(a):
    
17.     list1 = []
    
18.     while True:
    
19.         tmp = a
    
20.         i = 2
    
21.         while True:
    
22.             if Prime(i):
    
23.                 if tmp % i == 0:
    
24.                     if i not in list1:
    
25.                         list1.append(i)
    
26.                     tmp = tmp/i
    
27.                     i = 2
    
28.                     continue
    
29.                 elif Prime(tmp) or tmp == 1:
    
30.                     if tmp not in list1:
    
31.                         if Prime(tmp):
    
32.                             list1.append(tmp)
    
33.                     return list1
    
34.                 else:
    
35.                     i += 1
    
36.             else:
    
37.                 i += 1
    
38. 
    
39. 
    
40. i = 2
    
41. while True:
    
42.     if not Prime(i):
    
43.         if len(cycle(i)) == 4 and len(cycle(i+1)) == 4 and len(cycle(i+2)) == 4 and len(cycle(i+3)) == 4:
    
44.             print(i)
    
45.             break
    
46.         else:
    
47.             i += 1
    
48.     else:
    
49.         i += 1

复制代码

答案:134043

渡风

此代码使用matlab编程
Problem47所用时间为691.2743秒
Problem47的答案为134043

1. %% Problem47.m
   
2. %最后编辑时间：2017-06-13 13:16 版本1.0
   
3. %找出连续四个数，四个数的没有相同的质因子，输出它们中的第一个
   
4. 
   
5. % Problem47所用时间为691.2743秒
   
6. % Problem47的答案为134043
   
7. 
   
8. function Output = Problem47()
   
9. tic
   
10. 
    
11. Start = 644;
    
12. Meet = 0;
    
13. process = 1;
    
14. while (Meet == 0)
    
15.     Meet = 1;
    
16.     Start = Start + process;
    
17.     disp(Start)
    
18.     AllNum = Start : Start + 3;
    
19.     for ii = 1:4
    
20.         if Isprime(AllNum(ii)) == 1
    
21.             process = ii;
    
22.             Meet = 0;
    
23.             break
    
24.         else
    
25.             if RealFactorNum(AllNum(ii)) < 4
    
26.                process = ii;
    
27.                Meet = 0;
    
28.                break
    
29.             end
    
30.         end
    
31.     end      
    
32. end
    
33. 
    
34. Output = AllNum(1);
    
35. toc
    
36. 
    
37. disp('此代码使用matlab编程')
    
38. disp(['Problem47所用时间为',num2str(toc),'秒'])
    
39. disp(['Problem47的答案为',num2str(Output)])
    
40. end
    
41. %% Isprime.m
    
42. %判定一个数是否为质数
    
43. function Judge = Isprime(n)
    
44. if nargin == 0
    
45.     n = 64;
    
46. end
    
47. 
    
48. if n == 2
    
49.     Judge = 1;
    
50. else
    
51.     if mod(n,2) == 0
    
52.         Judge = 0;
    
53.     else
    
54.         Judge = 1;
    
55.         for ii = 3: 2 :floor(sqrt(n)) + 1
    
56.             if mod(n,ii) == 0
    
57.                 Judge = 0;
    
58.                 break
    
59.             end
    
60.         end
    
61.     end
    
62. end
    
63. end
    
64. 
    
65. %% 输入一个数得到其真因子的个数
    
66. % 若输入为质数，则返回0
    
67. function Output = RealFactorNum(n)
    
68. if nargin == 0
    
69. n = 64;
    
70. end
    
71. % if isprime(n) == 1
    
72. %     Output = 0;
    
73. % else
    
74. %    PossiFactor = 1:n-1;
    
75.     Judge = ones(1,n-1);
    
76.     Judge(1) = 0;
    
77.     for ii = 1:n-1
    
78.         if Judge(ii) == 1
    
79.             for jj = 2*ii : ii : n-1
    
80.                 Judge(jj) = 0;
    
81.             end
    
82.         end
    
83.     end
    
84.     PossiFactor = find(Judge == 1);
    
85.     PrimeFactor = [];
    
86.     while (isprime(n) == 0)
    
87.         for kk = 1:length(PossiFactor)
    
88.             if mod(n, PossiFactor(kk)) == 0
    
89.                 PrimeFactor = [PrimeFactor, PossiFactor(kk)];
    
90.                 n = n/PossiFactor(kk);
    
91.                 break
    
92.             end
    
93.         end
    
94.     end
    
95.     PrimeFactor = [PrimeFactor, n];
    
96.     Output = length(unique(PrimeFactor));
    
97. end
    
98.    
    
99.    

复制代码

najin

用的matlab：
结果是：
      134043

时间已过 4.581318 秒。
>>

zzzgod

1. #include <iostream>
   
2. #include <cmath>
   
3. 
   
4. using namespace std;
   
5. 
   
6. bool div(int num)
   
7. {
   
8.         int count=0;
   
9.         for(int i=2;i<=num;i++)
   
10.         {
    
11.                 if(num%i==0)
    
12.                 {
    
13.                         num/=i;
    
14.                         count++;
    
15.                         while(num%i==0)
    
16.                                 num/=i;
    
17.                 }
    
18.                 if(count==5)
    
19.                         return false;
    
20.         }
    
21.         if(count==4)
    
22.                 return true;
    
23.         return false;
    
24. }
    
25. 
    
26. int main()
    
27. {
    
28.         int count=0;
    
29.         for(int i=2;i<=100000000;i++)
    
30.         {
    
31.                 if(div(i))
    
32.                         count++;
    
33.                 else
    
34.                         count=0;
    
35.                 if(count==4)
    
36.                 {
    
37.                         cout<<i-3;
    
38.                         return 0;
    
39.                 }
    
40.         }
    
41. }

复制代码

k往事如烟k

134043

1. def isPrime(num):
   
2.     if num <= 1:
   
3.         return False
   
4.     else:
   
5.         for i in range(2, int(num**0.5+1)):
   
6.             if num % i == 0:
   
7.                 return False
   
8.         else:
   
9.             return True
   
10. 
    
11. def PrimeFactorList(num):
    
12.     if isPrime(num):
    
13.         return []
    
14.     else:
    
15.         PrimeFactorList = []
    
16.         for i in range(2, int(num**0.5+1)):
    
17.             while num % i == 0:
    
18.                 num //= i
    
19.                 PrimeFactorList.append(i)
    
20.                 if isPrime(num):
    
21.                     PrimeFactorList.append(num)
    
22.         return list(set(PrimeFactorList))
    
23. 
    
24. for num1 in range(2*3*5*7, 1000000):
    
25.     num2 = num1 + 1
    
26.     num3 = num1 + 2
    
27.     num4 = num1 + 3
    
28.     if len(PrimeFactorList(num1)) == 4:
    
29.         if len(PrimeFactorList(num2)) == 4:
    
30.             if len(PrimeFactorList(num3)) == 4:
    
31.                 if len(PrimeFactorList(num4)) == 4:
    
32.                     print(num1)
    
33.                     break

复制代码

王小召

借用楼上答案， 本题重点：对于大数据查找效率考量，速度从快到慢：set > dict > list 如果本题第16行用 if num in b_nums（也就是：in + list模式） 做替换，用时可能会多百倍，复杂度会从O(log N) 编程了 O(N)！!

结果是：{134043: [3, 7, 13, 491], 134044: [2, 2, 23, 31, 47], 134045: [5, 17, 19, 83], 134046: [2, 3, 3, 11, 677]}
用时：0.6084039 秒

1. 
   
2. import time
   
3. 
   
4. 
   
5. def prime(N):
   
6.     bool_num = [1] * N
   
7.     bool_num[0], bool_num[1] = 0, 0
   
8.     for i, j in enumerate(bool_num):
   
9.         if j:
   
10.             for x in range(i**2, N, i):
    
11.                 bool_num[x] = 0
    
12.     return [i for i, j in enumerate(bool_num) if j]
    
13. 
    
14. 
    
15. def composite_numbers(num, b_nums, s_nums, result = []):
    
16.     if num in s_nums:
    
17.         result.append(num)
    
18.         return [num]
    
19. 
    
20.     sqr_n = num**0.5 + 1
    
21.     for each_prime in b_nums:
    
22.         if not num % each_prime:
    
23.             result.append(each_prime)
    
24.             composite_numbers(int(num/each_prime), b_nums, s_nums, result)
    
25.             break
    
26.         elif each_prime > sqr_n:
    
27.             break
    
28.     return result
    
29. 
    
30. 
    
31. def cal_(N):
    
32.     b_nums = prime(N)
    
33.     s_nums = set(b_nums)
    
34.     d_result = {}
    
35.     for number in range(647, N):
    
36.         tmp = composite_numbers(number, b_nums, s_nums, [])
    
37.         if len(set(tmp)) == 4:
    
38.             d_result[number] = tmp
    
39.         else:
    
40.             d_result.clear()
    
41.         if len(d_result) == 4:
    
42.             return d_result
    
43. 
    
44. print("结果是：{}\n用时：{} 秒".format(cal_(150000), time.process_time()))
    

复制代码

永恒的蓝色梦想

1. from time import process_time as t
   
2. t1=t()
   
3. from itertools import count
   
4. def f(a):
   
5.   l=set()
   
6.   if a%2==0:
   
7.     while a%2==0:a//=2
   
8.     l.add(2)
   
9.   while a!=1:
   
10.     for i in range(3,int(a**0.5+1),2):
    
11.       if a%i==0:
    
12.         while a%i==0:a//=i
    
13.         l.add(i)
    
14.         break
    
15.     else:
    
16.       l.add(int(a))
    
17.       break
    
18.   return len(l)==4
    
19. c=count(1)
    
20. while 1:
    
21.   n=c.__next__()
    
22.   if f(n):
    
23.     l=n
    
24.     for i in range(3):
    
25.       if not f(c.__next__()):break
    
26.     else:break
    
27. print(f'运行结果：{l}\n运行时间：{t()-t1}s')

复制代码

运行结果：134043
运行时间：0.734375s

永恒的蓝色梦想

王小召 发表于 2019-6-24 14:28
借用楼上答案， 本题重点：对于大数据查找效率考量，速度从快到慢：set > dict > list 如果本题第16行用 if ...

老哥说实话我看不懂你的代码……

debuggerzh

134043

Process returned 0 (0x0)   execution time : 1.258 s
Press any key to continue.

1. #include<iostream>
   
2. #include<cstdio>
   
3. #include<set>
   
4. using namespace std;
   
5. 
   
6. const int M = 5e5;
   
7. int cnt = 0;
   
8. bool prime[M];
   
9. int factor[M];
   
10. 
    
11. void euler_sieve(){
    
12.   prime[1] = false;
    
13.   for (int i = 2;i <= M;i++) prime[i] = true;
    
14. 
    
15.   for (int i = 2;i <= M;i++){
    
16.     if (prime[i]) factor[++cnt] = i;
    
17.     for (int j = 1;j <= cnt && i*factor[j] < M;j++){
    
18.       prime[i*factor[j]] = false;
    
19.       if (i % factor[j] == 0) break;
    
20.     }
    
21.   }
    
22. }
    
23. 
    
24. void solve(int x,set<int> & s){
    
25.   for (int i = 1;x != 1;i++){
    
26.     while (x % factor[i] == 0){
    
27.       s.insert(factor[i]);
    
28.       x /= factor[i];
    
29.     }
    
30.   }
    
31. }
    
32. 
    
33. bool judge(set<int> s[]){
    
34.   for (int i = 0;i < 4;i++)
    
35.     if (s[i].size() != 4) return false;
    
36. 
    
37.   for (int i = 0;i < 3;i++){
    
38.     for (int j = i+1;j < 4;j++){
    
39.       if (s[i] == s[j]) return false;
    
40.     }
    
41.   }
    
42.   return true;
    
43. }
    
44. 
    
45. int main(){
    
46.   euler_sieve();
    
47. 
    
48.   for (int i = 2;;i++){
    
49.     set<int> s[4];
    
50.     if (prime[i] || prime[i+1] || prime[i+2] || prime[i+3]) continue;
    
51. 
    
52.     for (int j = 0;j < 4;j++) solve(i+j,s[j]);
    
53. 
    
54.     if (judge(s)) {cout << i << endl; break;}
    
55.   }
    
56.   return 0;
    
57. }
    

复制代码

4444567

1. #找到一个数的所有质因数，运用元组
   
2. #找满足要求的下一个数
   
3. #两个列表，一个装质因数，一个装整数，用长度判断
   
4. def prime(n):
   
5.     if n == 2 or n==3:
   
6.         return True
   
7.     if n%2==0:
   
8.         return False
   
9.     else:
   
10.         for i in range(3,int(n**0.5)+1,2):
    
11.             if n%i==0:
    
12.                 return False
    
13.         return True
    
14.    
    
15. def zys(n):
    
16.     l = []
    
17.     for i in range(2,int(n/30)+1):#因为找四个不同质因数，那前三个质因数最小是2、3、5
    
18.         if n%i == 0 and prime(i):
    
19.             l.append(i)
    
20.     return l
    
21. 
    
22. n = 647
    
23. while n < 150000:
    
24.     if len(zys(n)) ==4:
    
25.         if len(zys(n+1)) ==4:
    
26.             if len(zys(n+2)) ==4:
    
27.                 if len(zys(n+3)) ==4:
    
28.                     print(n,zys(n))
    
29.                     print(n+1,zys(n+1))
    
30.                     print(n+2,zys(n+2))
    
31.                     print(n+3,zys(n+3))
    
32.                     break
    
33.     n += 1
    
34. 
    
35.    
    

复制代码

是有点菜，但是幸好可以算出来

a1351468657

迷雾少年 发表于 2016-8-30 21:56
怎么我出来了好多。。
16:2 2 2 2
24:2 2 2 3

你再看看题目

a1351468657

1. #include <stdio.h>
   
2. #include <time.h>
   
3. 
   
4. int check_num(int);
   
5. int check_num(int num)//判断是否有四个不同的质数因子
   
6. {
   
7.         int i, j, k, count = 0;
   
8.         k = num;
   
9.         for (i = 2; i <= num / 2; i++)
   
10.         {
    
11.                 if (k % i == 0)
    
12.                 {
    
13.                         count++;
    
14.                 }
    
15.                 while (k % i == 0)
    
16.                 {
    
17.                         k /= i;
    
18.                 }               
    
19.         }
    
20.         if (count == 4)
    
21.         {
    
22.                 return 1;
    
23.         }
    
24.         else
    
25.         {
    
26.                 return 0;
    
27.         }
    
28. }
    
29. 
    
30. main()
    
31. {
    
32.         int i = 644, j;
    
33.         int begin = time(NULL), end;
    
34. 
    
35.         while (1)
    
36.         {
    
37.                 i++;
    
38.                 j = i;
    
39.                 if (check_num(j))
    
40.                 {
    
41.                         if (check_num(++j))
    
42.                         {
    
43.                                 if (check_num(++j))
    
44.                                 {
    
45.                                         if (check_num(++j))
    
46.                                         {
    
47.                                                 printf("%d\n", i);
    
48.                                                 break;
    
49.                                         }
    
50.                                 }
    
51.                         }
    
52.                 }
    
53.                
    
54.         }
    
55.         end = time(NULL);
    
56.         printf("time = %ds", end - begin);
    
57. }

复制代码

134043
time = 28s

guosl

1. /*
   
2. 欧拉问题47
   
3. 答案：134043
   
4. 耗时：0.0228895秒
   
5. */
   
6. #include <iostream>
   
7. #include <algorithm>
   
8. #include <cstring>
   
9. #include <cmath>
   
10. using namespace std;
    
11. 
    
12. const int INF = 0x7fffffff;
    
13. const int nSteps = 1000;
    
14. 
    
15. bool chkpcount(int n)//计算n中的素因数的个数是否是4
    
16. {
    
17.   if (n == 1 || n == 2)
    
18.     return false;
    
19.   int nCount = 0;//记录素因数的个数
    
20.   if (n % 2 == 0)
    
21.   {
    
22.     ++nCount;
    
23.     while (n % 2 == 0)
    
24.       n /= 2;
    
25.   }
    
26.   int k = 3;
    
27.   while (n > 1 && nCount <= 4)
    
28.   {
    
29.     int d = (int)sqrt((double)n);
    
30.     bool bFlag = true;
    
31.     for (; k <= d; k += 2)//找出当前n的最小素因数
    
32.     {
    
33.       if ((n % k) == 0)//找到一个比n小的素因数
    
34.       {
    
35.         bFlag = false;
    
36.         ++nCount;
    
37.         while (n % k == 0)//将n中的所有因数k除掉
    
38.           n /= k;
    
39.         break;
    
40.       }
    
41.     }
    
42.     if (bFlag)//没有找到比n小的素因数，故n本身是一个素数
    
43.     {
    
44.       ++nCount;
    
45.       break;
    
46.     }
    
47.   }
    
48.   return (nCount == 4);
    
49. }
    
50. 
    
51. int main(void)
    
52. {
    
53.   char p[nSteps + 5];
    
54.   int k = 647;
    
55.   int nResult = INF;
    
56.   p[nSteps + 4] = 0;
    
57.   while (true)
    
58.   {
    
59.     memset(p, '0', (nSteps + 4) * sizeof(char));
    
60.     for (int i = 0; i < nSteps + 4; ++i)//+4是为了解决前后重叠的问题
    
61.     {
    
62.       int j = i + k;
    
63.       if (chkpcount(j))
    
64.         p[i] = '1';//标记哪些数是有4个素数因数
    
65.     }
    
66.     char* pi = strstr(p, "1111");//查找是否有4个连续的存在
    
67.     if (pi != NULL)//存在
    
68.     {
    
69.       int ip = int(pi - p) + k;
    
70.       nResult = min(nResult, ip);//记录
    
71.       break;//退出循环
    
72.     }
    
73.     else
    
74.       k += nSteps;
    
75.   }
    
76.   cout << nResult << endl;
    
77.   return 0;
    
78. }

复制代码

叶落了

#利用一个数可以表示唯一的质数积形式，迭代后看范围类是否只有4个质数满足
#include<stdio.h>
#include<math.h>

int determine0(int i)
{
        int l;
        l=pow(i,0.5);
        for(int m=2;m<=l;m++)
        {
                if(i%m==0)
                {
                        return 0;
                }
        }
        return 1;
}

int determine1(int i)
{
        int num=0;
        //int k=pow(i,0.5);
        int k=i/30;
        for(int m=2;m<=k;m++)//for(int m=2;m<=k;m++)想错了，这里的范围弄错  了
        {
                if(determine0(m)&&i%m==0)
                {
                        num=num+1;
                        if(num>4)
                        {
                                return 0;
                        }
                }
        }
        if(num==4)
        {
                return 1;
        }
        return 0;
}
int main(void)
{
        int i=1;
        int x=100;
        while(i)
        {
                if(determine1(x)&&determine1(x+1)&&determine1(x+2)&&determine1(x+3))
                {
                        i=0;
                        printf("%d\n",x);
                }
                x=x+1;
        }
        return 0;
}