gcd，lcm，ext_gcd，inv

Least Common Multiple http://acm.hdu.edu.cn/showproblem.php?pid=1019

 #include<cstdio>
 int gcd(int a,int b){
     return b?gcd(b,a%b):a;
 }
 int lcm(int a,int b){
     return a/gcd(a,b)*b;
 }
 int main(){
     int n,m,ans,x;
     while(~scanf("%d",&n)){
         while(n--){
             ans=;
             scanf("%d",&m);
             while(m--){
                 scanf("%d",&x);
                 ans=lcm(ans,x);
             }
             printf("%d\n",ans);
         }
     }
     return ;
 }

Turn the pokers http://acm.hdu.edu.cn/showproblem.php?pid=4869

 #include<cstdio>
 #include<algorithm>
 using namespace std;
 typedef __int64 LL;
 const int M=;
 const int mod=;
 int ext_gcd(int a,int b,int &x,int &y){//扩展gcd d=gcd(a,b)=a*x+b*y; return d,x,y;
     int t,ret;
     if(!b){
         x=;
         y=;
         return a;
     }
     ret=ext_gcd(b,a%b,x,y);
     t=x;
     x=y;
     y=t-a/b*y;
     return ret;
 }
 int inv(int a,int b,int c){//ext_gcd求逆元 (b/a)%c
     int x,y;
     ext_gcd(a,c,x,y);
     return (1LL*x*b%c+c)%c;
 }
 LL C[M];
 LL INV[M];
 int main() {
     for(int i=; i<M; i++) {
         INV[i]=inv(i,,mod);
     }
     int n,m,a;
     while(~scanf("%d%d",&n,&m)) {
         C[]=;
         for(int i=;i<=m;i++){
             C[i]=C[i-]*(m-i+)%mod*INV[i]%mod;
         }
         int L=,R=,nl,nr,tmp;
         for(int i=;i<n;i++){
             scanf("%d",&a);
             tmp=min(m-L,a);
             nr=L+tmp-(a-tmp);
             tmp=min(R,a);
             nl=R-tmp+(a-tmp);
             if(nl>nr) swap(nl,nr);
             if(L<=a&&a<=R){
                 if(L%==a%){
                     nl=;
                 }
                 else{
                     nl=min(nl,);
                 }
             }
             if((m-R)<=a&&a<=(m-L)){
                 if((m-L)%==a%){
                     nr=m;
                 }
                 else{
                     nr=max(nr,m-);
                 }
             }
             if(L>=a) nl=min(nl,L-a);
             if(m-R>=a) nr=max(nr,R+a);
             L=nl;
             R=nr;
         }
         int ans=;
         for(int i=L;i<=R;i+=){
             ans+=C[i];
             ans%=mod;
         }
         printf("%d\n",ans);
     }
     return ;
 }

 #include<cstdio>
 #include<algorithm>
 using namespace std;
 typedef __int64 LL;
 const int M=;
 const int mod=;
 LL C[M];
 LL INV[M];
 void inv_init(){//初始化%mod的乘法逆元
     INV[]=;
     for(int i=;i<M;i++){
         INV[i]=INV[mod%i]*(mod-mod/i)%mod;
     }
 }
 int main() {
     inv_init();
     int n,m,a;
     while(~scanf("%d%d",&n,&m)) {
         C[]=;
         for(int i=;i<=m;i++){
             C[i]=C[i-]*(m-i+)%mod*INV[i]%mod;
         }
         int L=,R=,nl,nr,tmp;
         for(int i=;i<n;i++){
             scanf("%d",&a);
             tmp=min(m-L,a);
             nr=L+tmp-(a-tmp);
             tmp=min(R,a);
             nl=R-tmp+(a-tmp);
             if(nl>nr) swap(nl,nr);
             if(L<=a&&a<=R){
                 if(L%==a%){
                     nl=;
                 }
                 else{
                     nl=min(nl,);
                 }
             }
             if((m-R)<=a&&a<=(m-L)){
                 if((m-L)%==a%){
                     nr=m;
                 }
                 else{
                     nr=max(nr,m-);
                 }
             }
             if(L>=a) nl=min(nl,L-a);
             if(m-R>=a) nr=max(nr,R+a);
             L=nl;
             R=nr;
         }
         int ans=;
         for(int i=L;i<=R;i+=){
             ans+=C[i];
             ans%=mod;
         }
         printf("%d\n",ans);
     }
     return ;
 }

end