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=100010;
int a[M];
LL b[M],c[M];
int val[M];
vector<int> yueshu[M];
int main(){
    for(int i=1;i<=100000;i++){
        yueshu[i].clear();
        for(int j=1;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=0;i<n;i++){
            scanf("%d",&a[i]);
        }
        mt(val,0);
        for(int i=n-1;i>=0;i--){
            int now=a[i];
            if(!val[now]){
                b[i]=now;
            }
            else{
                b[i]=val[now];
            }
            int len=yueshu[now].size();
            for(int j=0;j<len;j++){
                int yue=yueshu[now][j];
                val[yue]=now;
            }
        }
        mt(val,0);
        for(int i=0;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=0;j<len;j++){
                int yue=yueshu[now][j];
                val[yue]=now;
            }
        }
        LL ans=0;
        for(int i=0;i<n;i++){
            ans+=b[i]*c[i];
        }
        printf("%I64d
",ans);
    }
    return 0;
}

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=6;
const int M=8;
class Matrix{
    typedef int typev;
public:
    typev val[M][M];
    void zero(){
        memset(val,0,sizeof(val));
    }
    void unit(){
        zero();
        for(int i=0;i<M;i++) val[i][i]=1;
    }
}D;
int LEN;
Matrix operator * (const Matrix &a,const Matrix &b){
    Matrix tmp;
    tmp.zero();
    for(int k=0;k<LEN;k++)
        for(int i=0;i<LEN;i++){
            if(a.val[i][k])
                for(int j=0;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&1) tmp=tmp*a;
        a=a*a;
        p>>=1;
    }
    return tmp;
}
struct mat {
    int n,m;
    int data[1024][1024];
}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 0;
        c.n=a.n,c.m=b.m;
        for (i=0; i<c.n; i++)
            for (j=0; j<c.m; j++)
                for (c.data[i][j]=k=0; k<a.m; k++)
                    c.data[i][j]+=a.data[i][k]*b.data[k][j];
        return 1;
    }
} 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=0;i<N;i++){
            for(int j=0;j<K;j++){
                scanf("%d",&A.data[i][j]);
            }
        }
        for(int i=0;i<K;i++){
            for(int j=0;j<N;j++){
                scanf("%d",&B.data[i][j]);
            }
        }
        gx.mul(C,B,A);
        LEN=C.n;
        for(int i=0;i<C.n;i++){
            for(int j=0;j<C.m;j++){
                D.val[i][j]=C.data[i][j];
            }
        }
        D=D^(N*N-1);
        for(int i=0;i<C.n;i++){
            for(int j=0;j<C.m;j++){
                C.data[i][j]=D.val[i][j];
            }
        }
        gx.mul(E,A,C);
        gx.mul(C,E,B);
        LL ans=0;
        for(int i=0;i<C.n;i++){
            for(int j=0;j<C.m;j++){
                ans+=C.data[i][j]%6;
            }
        }
        printf("%I64d
",ans);
    }
    return 0;
}

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

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

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

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=100010;
LL lazy[M],sum[M];
int main(){
    int n,m,x,y,d;
    while(~scanf("%d",&n),n){
        mt(lazy,0);
        scanf("%d",&m);
        while(m--){
            scanf("%d%d%d",&x,&y,&d);
            lazy[x]+=d;
            lazy[y+1]-=d;
        }
        LL now=0;
        for(int i=1;i<=n;i++){
            now+=lazy[i];
            sum[i]=now;
        }
        for(int i=n-1;i>=1;i--){
            sum[i]+=sum[i+1];
        }
        scanf("%d",&m);
        LL h;
        int ans=0;
        while(m--){
            scanf("%I64d%d",&h,&x);
            if(h>sum[x]) ans++;
        }
        printf("%d
",ans);
    }
    return 0;
}

end