纯SAS宏程序打印杨辉三角，完全数，相亲数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

打印杨辉三角（Pascal's triangle），

1. %MACRO TRIANGLE(LINENo);
   
2.     %do i=1 %to &lineno;
   
3.         %let str=;
   
4.         %do j=1 %to &i;
   
5.             %if &j=1 or &j=&i %then  
   
6.                 %let x&i.y&j=1;
   
7.             %else
   
8.                 %do;
   
9.                 %let upline=%eval(&i-1);
   
10.                 %let up1=%eval(&j-1);
    
11.                 %let x&i.y&j=%eval(&&x&upline.y&up1+&&x&upline.y&j);
    
12.                 %end;
    
13.             %let str=&str%qsysfunc(repeat(%str( ),%eval(5-%length(&&x&i.y&j))))&&x&i.y&j;
    
14.         %end;
    
15.         %let str=%qsysfunc(repeat(%str( ),%eval((&lineNo-&j+1)*3)))%str(&str);
    
16.         %put %str(&str);
    
17.     %end;
    
18. %MEND;
    
19. %TRIANGLE(15);

复制代码

纯宏程序打印完全数、完美数(Perfect number)

1. %macro divisor(Num,print=0);
   
2.         %local allDivisor ;
   
3.         %do Divisor=1 %to %eval(&Num/2+1);
   
4.         %if  %sysfunc(mod(&Num,&divisor))=0 %then %let allDivisor=&alldivisor+&divisor;
   
5.         %end;
   
6.         %if &print %then %put  &allDivisor;
   
7.         %eval(0&allDivisor)
   
8. %mend;
   
9. 
   
10. 
    
11. 
    
12. %macro complete(N);
    
13.         %local one;
    
14.         %do  one=1 %to &n;
    
15.         %if &one=%divisor(&one) %then %put one=&one %divisor(&one,print=1);
    
16.         %end;
    
17. 
    
18. %mend;
    
19. 
    
20. %complete(1000);

复制代码

纯宏程序打印相亲数（Amicable）

1. %macro divisor(Num,print=0);
   
2.         %local allDivisor ;
   
3.         %do Divisor=1 %to %eval(&Num/2+1);
   
4.         %if  %sysfunc(mod(&Num,&divisor))=0 %then %let allDivisor=&alldivisor+&divisor;
   
5.         %end;
   
6.         %if &print %then %put  &allDivisor;
   
7.         %eval(0&allDivisor)
   
8. %mend;
   
9. 
   
10. %macro love(N);
    
11.         %local one anotherone;
    
12.         %do  one=1 %to &n;
    
13.         %let anotherOne=%divisor(&one);
    
14.         %if   %divisor(&anotherOne)=&one and &one ne &anotherOne  %then %put one=&one %divisor(&one,print=1) anotherone=&anotherone %divisor(&Anotherone,print=1);
    
15.         %end;
    
16. %mend;
    
17. data _null_;
    
18.         call symput('ahuigefromtime',put(time(),time8.));
    
19. run;
    
20. 
    
21. %love(1500);
    
22. 
    
23. data _null_;
    
24.         diff=time()-input("&ahuigefromtime",time8.);
    
25.         put '################ time used :' diff time8.;
    
26. run;

复制代码

》》》》》》》》》》》》》》》》》》》》》》》》
ahuige

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：SAS宏程序 杨辉三角 SAS宏 完全数 宏程序 杨辉三角 程序 杨辉三角 完全数 相亲数 SAS宏程序

感谢分享