「算法竞赛入门经典」「第二章」

例题2-1 aabb(P20)

输出所有形似aabb的4位完全平方数,下面这种方法不用开根,如果使用开根的方法,比如利用if(sqrt(n)==floor(sqrt(n)))这种方法判断是否是完全平方数,因为浮点数的计算存在误差,所以要谨慎使用,一般牵涉到浮点数的floor运算,都要进行四舍五入,floor(x+0.5)

#include "stdafx.h"
#include <stdio.h>

int main(int argc, char* argv[])
{
    for(int i=0;;i++){
        int x = i*i;
        if(x<1000)continue;
        if(x>9999)break;
        int high = x / 100;
        int low = x %100;
        if(high/10==high%10 && low/10==low%10)
            printf("%d\n",x);
    }
    return 0;
}
例题2-2 aabb(P23)

对于任意大于1的自然数n,若n为奇数,则n=3n+1,否则n=n/2。经过若干次这样的变换,最终n一定会变为1。输入n,输出变换次数,n<=10^9。

我还在用VC6.0写程序,98年产生的VC6.0,long long类型是C99(1999年)中才新加的,所以并不支持,会报错。而且VC6.0应该是32位环境,所以用long也是不行的。

16位系统:long是4字节,int是2字节
32位系统:long是4字节,int是4字节
64位系统:long是8字节,int是4字节

__int64(两个下划线)则与long long相同,__int64对应I64d,long long对应lld

// lt22.cpp : Defines the entry point for the console application.
// 3n+1问题

#include "stdafx.h"
#include <stdio.h>

int main(int argc, char* argv[])
{
    __int64 n;
    int count=0;
    while(scanf("%I64d",&n)){
        count=0;
        while(1){
            if(n==1)break;
            if(n%2==1) n=3*n+1;
            else n = n/2;
            count++;
        }
        printf("%d\n",count);
    }
    return 0;
}

例题2-3 计算1-1/3+1/5-1/7+...,直到最后一项小于10^-6(P24)
// lt23.cpp : Defines the entry point for the console application.
//计算1-1/3+1/5-1/7+...直到最后一项小于10^-6

#include "stdafx.h"
#include <stdio.h>

int main(int argc, char* argv[])
{
    float x=1;
    float sum=0;
    int count=0;
    do{
        if(count%2==0){
            sum+=x;
        }
        else if(count%2==1){
            sum-=x;
        }
        count++;
        x = 1.0/(2*count+1); //如果此处写成1/(2*count+1)结果将永远是1
    }while(x>1e-6);
    printf("%f\n",sum);
    return 0;
}
例题2-4 阶乘之和(P25)

输入n,计算S=1!+2!+...+n!的末6位。n<=10^6。
陷阱就在阶乘很容易就会产生乘法溢出,所以,有一个规律。

要计算只包含加法,减法和乘法的整数表达式除以正整数n的余数,可以在每步计算之后对n取余,结果不变。

#include <stdio.h>

int main(int argc, char* argv[])
{
    int n;
    int sum=0;
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        int s=1;
        for(int j=1;j<=i;j++){
            s = (s*j%1000000);
        }
        sum = (sum+s)%1000000;
    }
    printf("%d\n",sum);
    return 0;
}
习题2-2 韩信点兵

果断暴力循环求解

#include <stdio.h>
int main(int argc, char* argv[])
{
    int a,b,c;
    int count=0;
    while(scanf("%d%d%d",&a,&b,&c)==3){
        count++;
        for(int i=10;i<=100;i++){
            if(i%3==a && i%5==b && i%7==c){
                printf("Case %d: %d",count,i);
                break;
            }
        }
        if(i==101)printf("Case %d: No answer",count);
    }
    return 0;
}
习题2-4 子序列之和

本题陷阱两处。首先要在代码中将被除数1写成浮点形式1.0,第二正整数的平方有可能会导致溢出,所以不能直接使用1.0/(m*m),要改成1.0/m/m,才没有问题。

#include <stdio.h>

int main(int argc, char* argv[])
{
    int n,m;
    int count=0;
    double sum=0;
    while(scanf("%d%d",&n,&m)&&(n||m)){
        sum=0;
        count++;
        for(int i=n;i<=m;i++){
            sum += 1.0/i/i;
        }
        printf("Case %d: %.5lf\n",count,sum);
    }
    return 0;
}
习题2-5 分数化小数

本题说要精确到小数点后c位,由于之前认为printf("%5.3f",x)中的‘5’和‘3’都只能是常数而不能是变量,所以曾尝试将a/b的结果乘以10的c次方再除以10的c次方,利用%g出去之后多余的零,但其实并不符合题意。后来了解到其中是可以有变量的,比如printf("%*.*lf",x,y,z)其中的第一个*对应后面的x,代表整个数输出宽度,第二个*对应后面的y,代表小数点后保留位数,而要输出的数是z。

#include <stdio.h>

int main(int argc, char* argv[])
{
    int a,b,c;
    int count=0;
    while(scanf("%d%d%d",&a,&b,&c)&&(a+b+c)){
        count++;
        double ans = (double)a/(double)b;
        printf("Case %d: %.*lf\n",count,c,ans);
    }
    return 0;
}
习题2-6 排列

暴力求解相当麻烦。所以设置数组,看最后数组和是否为9。很巧妙。

#include <stdio.h>

int main(int argc, char* argv[])
{
    int x,y,z;
    int a[10]={0};
    for(x=100;x<330;x++){
        y = x*2;
        z = x*3;
        a[x/100] = a[x/10%10] = a[x%10] = 1;
        a[y/100] = a[y/10%10] = a[y%10] = 1;
        a[z/100] = a[z/10%10] = a[z%10] = 1;
        int sum=0;
        for(int i=1;i<=9;i++){
            sum += a[i];
        }
        if(sum==9)
            printf("%d %d %d\n",x,y,z);
        for(i=0;i<=9;i++){
            a[i]=0;
        }
    }
    return 0;
}
思考题 题目2
#include <stdio.h>

int main(int argc, char* argv[])
{
    double i;
    for(i=0;i!=10&&i<11;i+=0.1){
        printf("%.1f\n",i);

    }
    return 0;
}

由上面代码输出结果看不出到底为什么会出现一直循环输出的错误。
但其实是因为循环中是浮点数运算,浮点数在计算机中采用二进制科学计数法来储存,绝大多数是不能精确表示的。
当两个浮点数都是由赋值得到的,如

a=3.3; 
b=3.3;

则a==b值为真。

而当a,b中至少有一个为经由运算得到,那么由于其精确度的原因,会出现这样的情况

a=1.7+2.7; 
b=4.4;

则a==b的值就不一定为真了。

若想进行粗略的比较,则可用

fabs(a-b)<1e-6 
floor(b+0.5)==floor(a+0.5) 

等限制误差的语句。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 158,847评论 4 362
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 67,208评论 1 292
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 108,587评论 0 243
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 43,942评论 0 205
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 52,332评论 3 287
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 40,587评论 1 218
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 31,853评论 2 312
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 30,568评论 0 198
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 34,273评论 1 242
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 30,542评论 2 246
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 32,033评论 1 260
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 28,373评论 2 253
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 33,031评论 3 236
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 26,073评论 0 8
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 26,830评论 0 195
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 35,628评论 2 274
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 35,537评论 2 269

推荐阅读更多精彩内容