2014 Multi-University Training Contest 9

官方解题报告：http://blog.sina.com.cn/s/blog_6bddecdc0102uzwm.html

Boring Sum http://acm.hdu.edu.cn/showproblem.php?pid=4961

这个二分也过了，on更快些。

 #include<cstdio>
 #include<cstring>
 #include<vector>
 #define mt(a,b) memset(a,b,sizeof(a))
 using namespace std;
 typedef __int64 LL;
 const int M=;
 int a[M];
 LL b[M],c[M];
 int val[M];
 vector<int> yueshu[M];
 int main(){
     for(int i=;i<=;i++){
         yueshu[i].clear();
         for(int j=;j*j<=i;j++){
             if(!(i%j)){
                 yueshu[i].push_back(j);
                 if(i/j!=j)
                     yueshu[i].push_back(i/j);
             }
         }
     }
     int n;
     while(~scanf("%d",&n),n){
         for(int i=;i<n;i++){
             scanf("%d",&a[i]);
         }
         mt(val,);
         for(int i=n-;i>=;i--){
             int now=a[i];
             if(!val[now]){
                 b[i]=now;
             }
             else{
                 b[i]=val[now];
             }
             int len=yueshu[now].size();
             for(int j=;j<len;j++){
                 int yue=yueshu[now][j];
                 val[yue]=now;
             }
         }
         mt(val,);
         for(int i=;i<n;i++){
             int now=a[i];
             if(!val[now]){
                 c[i]=now;
             }
             else{
                 c[i]=val[now];
             }
             int len=yueshu[now].size();
             for(int j=;j<len;j++){
                 int yue=yueshu[now][j];
                 val[yue]=now;
             }
         }
         LL ans=;
         for(int i=;i<n;i++){
             ans+=b[i]*c[i];
         }
         printf("%I64d\n",ans);
     }
     return ;
 }

Fast Matrix Calculation http://acm.hdu.edu.cn/showproblem.php?pid=4965

ababab矩阵等于a*（ba）^2*b,ba是kk的矩阵复杂度就k^3 log（n^2）

 #include<cstdio>
 #include<cstring>
 typedef __int64 LL;
 const int mod=;
 const int M=;
 class Matrix{
     typedef int typev;
 public:
     typev val[M][M];
     void zero(){
         memset(val,,sizeof(val));
     }
     void unit(){
         zero();
         for(int i=;i<M;i++) val[i][i]=;
     }
 }D;
 int LEN;
 Matrix operator * (const Matrix &a,const Matrix &b){
     Matrix tmp;
     tmp.zero();
     for(int k=;k<LEN;k++)
         for(int i=;i<LEN;i++){
             if(a.val[i][k])
                 for(int j=;j<LEN;j++){
                     tmp.val[i][j]+=a.val[i][k]*b.val[k][j];
                     tmp.val[i][j]%=mod;
                 }
         }
     return tmp;
 }
 Matrix operator ^ (Matrix &a,int p){
     Matrix tmp;
     tmp.unit();
     while(p){
         if(p&) tmp=tmp*a;
         a=a*a;
         p>>=;
     }
     return tmp;
 }
 struct mat {
     int n,m;
     int data[][];
 }A,B,C,E;
 class MatrixOp {
 public:
     int mul(mat& c,const mat& a,const mat& b) {  //c=a*b
         int i,j,k;
         if (a.m!=b.n) return ;
         c.n=a.n,c.m=b.m;
         for (i=; i<c.n; i++)
             for (j=; j<c.m; j++)
                 for (c.data[i][j]=k=; k<a.m; k++)
                     c.data[i][j]+=a.data[i][k]*b.data[k][j];
         return ;
     }
 } gx;
 int main(){
     int N,K;
     while(~scanf("%d%d",&N,&K),N|K){
         A.n=N;
         A.m=K;
         B.n=K;
         B.m=N;
         for(int i=;i<N;i++){
             for(int j=;j<K;j++){
                 scanf("%d",&A.data[i][j]);
             }
         }
         for(int i=;i<K;i++){
             for(int j=;j<N;j++){
                 scanf("%d",&B.data[i][j]);
             }
         }
         gx.mul(C,B,A);
         LEN=C.n;
         for(int i=;i<C.n;i++){
             for(int j=;j<C.m;j++){
                 D.val[i][j]=C.data[i][j];
             }
         }
         D=D^(N*N-);
         for(int i=;i<C.n;i++){
             for(int j=;j<C.m;j++){
                 C.data[i][j]=D.val[i][j];
             }
         }
         gx.mul(E,A,C);
         gx.mul(C,E,B);
         LL ans=;
         for(int i=;i<C.n;i++){
             for(int j=;j<C.m;j++){
                 ans+=C.data[i][j]%;
             }
         }
         printf("%I64d\n",ans);
     }
     return ;
 }

Improving the GPA  http://acm.hdu.edu.cn/showproblem.php?pid=4968

4个for循环，突出一个暴力0ms，硬是要写dfs，当练习了，15ms。

 #include<cstdio>
 double big,sma;
 int num[],judge;
 int tobig[]={,,,,};
 int tosma[]={,,,,};
 double val[]={4.0,3.5,3.0,2.5,2.0};
 void dfs(int id,int sum){
     if(id==){
         num[]=sum;
         int bigsum=,smasum=;
         for(int i=;i<;i++){
             bigsum+=tobig[i]*num[i];
             smasum+=tosma[i]*num[i];
         }
         if(smasum<=judge&&judge<=bigsum){
             double now=;
             for(int i=;i<;i++){
                 now+=val[i]*num[i];
             }
             big=big>now?big:now;
             sma=sma<now?sma:now;
         }
         return ;
     }
     for(int i=;i<=sum;i++){
         num[id]=i;
         dfs(id+,sum-i);
     }
 }
 int main(){
     int t,n,avg;
     while(~scanf("%d",&t)){
         while(t--){
             scanf("%d%d",&avg,&n);
             sma=1e9;
             big=;
             judge=avg*n;
             dfs(,n);
             printf("%.4f %.4f\n",sma/n,big/n);
         }
     }
     return ;
 }

Killing Monsters http://acm.hdu.edu.cn/showproblem.php?pid=4970

这个题目是来卡nlogn的。左++右--处理离线区间覆盖问题。on。

 #include<cstdio>
 #include<cstring>
 #define mt(a,b) memset(a,b,sizeof(a))
 typedef __int64 LL;
 const int M=;
 LL lazy[M],sum[M];
 int main(){
     int n,m,x,y,d;
     while(~scanf("%d",&n),n){
         mt(lazy,);
         scanf("%d",&m);
         while(m--){
             scanf("%d%d%d",&x,&y,&d);
             lazy[x]+=d;
             lazy[y+]-=d;
         }
         LL now=;
         for(int i=;i<=n;i++){
             now+=lazy[i];
             sum[i]=now;
         }
         for(int i=n-;i>=;i--){
             sum[i]+=sum[i+];
         }
         scanf("%d",&m);
         LL h;
         int ans=;
         while(m--){
             scanf("%I64d%d",&h,&x);
             if(h>sum[x]) ans++;
         }
         printf("%d\n",ans);
     }
     return ;
 }

end