uva10976_整数拆分（优化至3ms）

最新推荐文章于 2020-03-12 21:05:58 发布

tt2767

最新推荐文章于 2020-03-12 21:05:58 发布

阅读量1k

点赞数

分类专栏： ACM_做过的一些题文章标签： uva ACM 数论

本文链接：https://blog.csdn.net/tt2767/article/details/47731225

版权

ACM_做过的一些题专栏收录该内容

35 篇文章 0 订阅

订阅专栏

///
作者：tt2767
声明：本文遵循以下协议自由转载-非商用-非衍生-保持署名|Creative Commons BY-NC-ND 3.0
查看本文更新与讨论请点击：http://blog.csdn.net/tt2767
链接被删请百度： CSDN tt2767
///

PS：刘汝佳代码仓库中这题的代码逗了，我交上去居然RE了。。。。（逃

题解：
原书的意思是找出最大值于最小值枚举就好了，然而仔细观察的话就会发现，第一组分母较小的数的分母为k+1,而其他组数据这个数都不是素数，最后数据分母必为2*k，那么我们只要在非素数中查找符合条件的值即可。

PS：我又写了一份暴力的，果断T了。。。。。。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
//
#include<iostream>
#include<algorithm>
#include<string>
#include <iterator>
#include<sstream>
#include<functional>
#include<numeric>
///
#include<vector>
#include<map>
#include <stack>
#include<queue>
#include<set>
#include <bitset>
#include <list>
using namespace std;
#define lch(x)  ((x) << 1)
#define rch(x) ((x)<<1|1)
#define dad(x) ((x)>>1)
#define lowbit(x) ((x)&(-x))
typedef  long long int LL;
const int INF = ~0U>>1;
const double eps = 1e-6;
const long double PI = acos(0.0) * 2.0;
const int N = 10 + 100000;
int prime[N] = {0};
bool not_prime[N] = {1,1};
int make_prime(int n);
int b[N],a[N];
int main()
{
    //ios::sync_with_stdio(false);
#ifdef ONLINE_JUDGE
#else
    freopen("in.txt", "r", stdin);
    freopen("out.txt", "w", stdout);
#endif
    int k;
    make_prime(N);
    while(scanf("%d",&k)==1)
    {
        if(k==1)
        {
            puts("1");
            printf("1/1 = 1/2 + 1/2\n");
            continue;
        }

        int MAX = k*(k+1);
        int n = 0 ;
        for(int i  = k+2; i<2*k;i++)
            if(not_prime[i])   b[n++]=i;

        int cnt = 2,j;
        for( j = 0 ; j<n ;j++)
        {
            int x = b[j],y;
            y = x*k/(x-k);
            if( 1.0*x*y/(x*1.0+y*1.0)==1.0*k )
            {
                a[j] = y;
                cnt++;
            }
            else
                b[j]=a[j]=-1;
        }
        printf("%d\n",cnt);
        printf("1/%d = 1/%d + 1/%d\n",k,MAX,k+1);
        for(int i = 0; i<j ; i++)
            if(a[i]!=-1&&b[i]!=-1)
                printf("1/%d = 1/%d + 1/%d\n",k,a[i],b[i]);

        printf("1/%d = 1/%d + 1/%d\n",k,2*k,k*2);
    }
    return 0;
}

int make_prime(int n) //线性筛素数
{
    int num = 0;
    for(int i = 2 ; i < n ; i++)
    {
        if(!(not_prime[i]))
            prime[num++] = i;

        for(int j = 0 ; j < num && i*prime[j] < n ; j++)
        {
            not_prime[ prime[j]*i ] = 1;
            if(!(i%prime[j]))    break; //只筛去最小素数因子，以免重复筛去合数
        }
    }
    return num;
}