• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>

            ArcTan

            dfs
            隨筆 - 16, 文章 - 117, 評論 - 6, 引用 - 0
            數據加載中……

            poj 2689(篩法求素數的應用)

            http://poj.org/problem?id=2689

            篩法篩去區間內的合數。
            開哥說是水題一個!不過,我WA了一天。二分查找的錯!!!
            141MS AC
            #include<stdio.h>
            #include
            <string.h>
            #include
            <math.h>
            #define inf 500000
            unsigned prime[
            50005],b[5000005],a[500005],tot;
            unsigned L,U;
            int b_search(int l,int r,int x)
            {
                
            int mid;
                
            while (l<=r)
                {
                    mid
            =(l+r)/2;
                    
            if (prime[mid]==x)
                        
            return mid;
                    
            else
                        
            if (prime[mid]<x)
                            l
            =mid+1;
                        
            else
                            r
            =mid-1;
                }
                
            return l;
            }
            int Getprime()
            {
                
            int i,j;
                memset(b,
            0,sizeof(b));
                tot
            =0;
                i
            =2;
                
            while (i<inf)
                {
                    
            while (b[i])    i++;
                    prime[
            ++tot]=i;
                    j
            =i;
                    
            while (j<=inf)
                    {
                        b[j]
            =1;
                        j
            +=i;
                    }
                }
                tot
            --;
            }
            int main()
            {
                unsigned 
            int i,j,k,total;
                unsigned 
            int max,min,maxi,mini;
                Getprime();

                
            while (scanf("%d%d",&L,&U)==2)
                {
                    memset(b,
            0,sizeof(b));
                    k
            =b_search(1,tot,(int)sqrt(U*1.0)+1.0);

                    i
            =1;
                    
            while (i<=k+1)     //這里改為k+1就AC了!!!嗚嗚嗚啎
                    {
                        j
            =L%prime[i];
                        
            if (j>0)
                            j
            =prime[i]-j;
                        
            while (j<=U-L)
                        {
                            
            if (j+L!=prime[i])
                                b[j]
            =1;
                            j
            +=prime[i];
                        }
                        i
            ++;
                    }

                    total
            =0;
                    
            for (i=0;i<=U-L ;i++ )
                        
            if (!b[i]&&L+i>=2)
                            a[
            ++total]=L+i;

                    
            if (total<2)
                        printf(
            "There are no adjacent primes.\n");
                    
            else
                    {
                        max
            =min=a[2]-a[1];maxi=1;mini=1;
                        
            for (i=3;i<=total ;i++ )
                        {
                            
            if (a[i]-a[i-1]>max)
                            {
                                max
            =a[i]-a[i-1];
                                maxi
            =i-1;
                            }
                            
            if (a[i]-a[i-1]<min)
                            {
                                min
            =a[i]-a[i-1];
                                mini
            =i-1;
                            }
                        }
                        printf(
            "%d,%d are closest, %d,%d are most distant.\n",a[mini],a[mini+1],a[maxi],a[maxi+1]);
                    }
                }
                
            return 0;
            }

            posted on 2012-04-22 20:41 wangs 閱讀(404) 評論(0)  編輯 收藏 引用 所屬分類: ACM-數學

            久久九九亚洲精品| 狠狠色丁香久久婷婷综合图片| 久久不见久久见免费视频7| 久久久久久九九99精品| 色综合久久精品中文字幕首页| 欧美伊人久久大香线蕉综合69| 久久久久亚洲精品天堂| 久久久久99精品成人片| 亚洲日本va中文字幕久久| 伊人久久免费视频| 久久99精品久久久久久久不卡| 日产久久强奸免费的看| 国产精品一区二区久久精品| 亚洲精品无码久久久久AV麻豆| 久久99国产综合精品| 久久天天婷婷五月俺也去| 伊人久久大香线焦综合四虎| 日韩人妻无码一区二区三区久久| 久久综合色区| 国产精品久久久久乳精品爆| 久久99精品久久久久久动态图| 亚洲精品第一综合99久久| 色噜噜狠狠先锋影音久久| 国内精品人妻无码久久久影院 | A级毛片无码久久精品免费 | 久久久一本精品99久久精品88| 久久高清一级毛片| 国产成人久久久精品二区三区 | 丁香狠狠色婷婷久久综合| 久久精品中文字幕一区| 亚洲乱码日产精品a级毛片久久| 国产精品成人久久久久三级午夜电影| 国内精品久久久久影院一蜜桃 | 久久天天躁狠狠躁夜夜不卡| 97久久精品人人做人人爽| 欧美一区二区精品久久| 99久久精品免费看国产| 久久成人精品| 久久精品国产久精国产果冻传媒 | 久久久久人妻精品一区| 国产成人精品免费久久久久|