2014 ACM/ICPC Asia Regional Shanghai Online

Tree http://acm.hdu.edu.cn/showproblem.php?pid=5044

树链剖分，区间更新的时候要用on的左++右--的标记方法，要手动扩栈，用c++交，综合以上的条件可过。

#include<cstdio>
#include<cstring>
#include<algorithm>
#pragma comment(linker, "/STACK:36777216")
#define mt(a,b) memset(a,b,sizeof(a))
using namespace std;
typedef __int64 LL;
const int M=100010;
struct G{
    struct E{
        int v,next;
    }e[M<<1];
    int le,head[M];
    void init(){
        le=0;
        mt(head,-1);
    }
    void add(int u,int v){
        e[le].v=v;
        e[le].next=head[u];
        head[u]=le++;
    }
}g;
int n,val[M],fa[M],dep[M],num[M],son[M],top[M],sid[M];
void dfs(int u){
    num[u]=1;
    son[u]=0;
    for(int i=g.head[u];~i;i=g.e[i].next){
        int v=g.e[i].v;
        if(v!=fa[u]){
            fa[v]=u;
            dep[v]=dep[u]+1;
            dfs(v);
            if(num[son[u]]<num[v]) son[u]=v;
            num[u]+=num[v];
        }
    }
}
void get(int u,int Top){
    sid[u]=++n;
    top[u]=Top;
    if(son[u]) get(son[u],top[u]);
    for(int i=g.head[u];~i;i=g.e[i].next){
        int v=g.e[i].v;
        if(v!=fa[u]&&v!=son[u]){
            get(v,v);
        }
    }
}
LL lazy[M],ans[M];
void query(){
    LL now=0;
    for(int i=1;i<=n;i++){
        now+=lazy[i];
        ans[i]=now;
    }
}
void update(int x,int y,int z){
    lazy[x]+=z;
    lazy[y+1]-=z;
}
void worknode(int x,int y,int z){
    int tx=top[x],ty=top[y];
    while(tx!=ty){
        if(dep[tx]<dep[ty]){
            swap(tx,ty);
            swap(x,y);
        }
        update(sid[tx],sid[x],z);
        x=fa[tx];
        tx=top[x];
    }
    if(dep[x]>dep[y]) swap(x,y);
    update(sid[x],sid[y],z);
}
void workedge(int x,int y,int z){
    int tx=top[x],ty=top[y];
    while(tx!=ty){
        if(dep[tx]<dep[ty]){
            swap(tx,ty);
            swap(x,y);
        }
        update(sid[tx],sid[x],z);
        x=fa[tx];
        tx=top[x];
    }
    if(x==y) return ;
    if(dep[x]>dep[y]) swap(x,y);
    update(sid[son[x]],sid[y],z);
}
struct IN{
    int type,u,v,w;
}in[M];
struct EDGE{
    int u,v;
}edge[M];
int main(){
    int t,N,m,u,v;
    while(~scanf("%d",&t)){
        int cas=1;
        while(t--){
            scanf("%d%d",&N,&m);
            g.init();
            for(int i=0;i<N-1;i++){
                scanf("%d%d",&u,&v);
                g.add(u,v);
                g.add(v,u);
                edge[i].u=u;
                edge[i].v=v;
            }
            n=fa[1]=dep[1]=num[0]=0;
            dfs(1);
            get(1,1);
            char op[8];
            for(int i=0;i<m;i++){
                scanf("%s%d%d%d",op,&in[i].u,&in[i].v,&in[i].w);
                if(op[3]=='1') in[i].type=1;
                else           in[i].type=2;
            }
            mt(lazy,0);
            for(int i=0;i<m;i++){
                if(in[i].type&1) worknode(in[i].u,in[i].v,in[i].w);
            }
            printf("Case #%d:
",cas++);
            query();
            for(int i=1;i<=N;i++){
                if(i>1) printf(" ");
                printf("%I64d",ans[sid[i]]);
            }
            puts("");
            mt(lazy,0);
            for(int i=0;i<m;i++){
                if(in[i].type==2) workedge(in[i].u,in[i].v,in[i].w);
            }
            query();
            for(int i=0;i<N-1;i++){
                int u=edge[i].u;
                int v=edge[i].v;
                if(dep[u]<dep[v]){
                    swap(u,v);
                }
                if(i) printf(" ");
                printf("%I64d",ans[sid[u]]);
            }
            puts("");
        }
    }
    return 0;
}

Contest http://acm.hdu.edu.cn/showproblem.php?pid=5045

根据题目要求，任意时刻都不能出现两个人做题数量的差大于1，也就是差要<=1，也就是m个问题要分成长度为n的等长的几块，每块里面一人做一个。

问题转化为n个人n题目，如何选一种排列来获得最大价值，若暴力枚举，复杂度是10！，所以用状态压缩，降为2^10，从000推到111.

转移时看哪个人没做过题，那么就可以让他做，状态就是他这位变为1，值就是当前状态的值加上这个人做当前的题目可获得的价值，当前的题目正好等于已做过题的人的个数+1.

#include<cstdio>
#include<cstring>
#include<algorithm>
#define mt(a,b) memset(a,b,sizeof(a))
using namespace std;
double a[16][1024],dp[1<<10];
int one(int x){
    int res=0;
    while(x){
        if(x&1) res++;
        x>>=1;
    }
    return res;
}
int main(){
    int t,n,m;
    while(~scanf("%d",&t)){
        int cas=1;
        while(t--){
            scanf("%d%d",&n,&m);
            for(int i=0;i<n;i++){
                for(int j=0;j<m;j++){
                    scanf("%lf",&a[i][j]);
                }
            }
            int s=0,num,all=1<<n;
            bool flag=true;
            double ans=0;
            while(flag){
                flag=false;
                if(s+n<m){
                    flag=true;
                    num=n;
                }
                else{
                    num=m-s;
                }
                mt(dp,0);
                for(int i=0;i<all;i++){
                    int yi=one(i);
                    if(yi>=num) continue;
                    for(int j=0;j<n;j++){
                        if((i>>j)&1) continue;
                        int next=i|(1<<j);
                        dp[next]=max(dp[next],dp[i]+a[j][s+yi]);
                    }
                }
                double big=0;
                if(num==n){
                    big=dp[all-1];
                }
                else{
                    for(int i=0;i<all;i++){
                        if(one(i)==num){
                            big=max(big,dp[i]);
                        }
                    }
                }
                ans+=big;
                s+=n;
            }
            printf("Case #%d: %.5f
",cas++,ans);
        }
    }
    return 0;
}

Sawtooth http://acm.hdu.edu.cn/showproblem.php?pid=5047

根据另一个类似的题目可以猜想出，这种分割方法出来的平面个数应该是个一元二次函数，a*x^2+b*x+c，过3个点，0的时候是1个平面，题目还输入两组，三个进去解出方程的abc，然后输入x就可以输出值了。

#include<cstdio>
#include<cstring>
#define mt(a,b) memset(a,b,sizeof(a))
const int M=32;
class Hp { //高精度类
    int len,s[M];
public:
    void init(char ch[]) {
        len=1;
        mt(s,0);
        int i=0;
        while(ch[i]=='0'&&ch[i]!=0) i++;
        if(ch[i]!=0) {
            len=strlen(ch)-i;
            for(i=0; i<len; i++) {
                s[i]=ch[len-i-1]-48;
            }
        }
    }
    void print() { //输出
        int i=len-1;
        while(s[i]==0&&i>0) i--;
        for(; i>=0; i--) {
            printf("%d",s[i]);
        }
    }
    void add(int x) { //高精度加单精度
        int temp=0;
        s[0]+=x;
        while(s[temp]>9) {
            s[temp]-=10;
            temp++;
            s[temp]++;
        }
        if(s[len]!=0) len++;
    }
    void subtract(Hp a) { //高精度减高精度
        for(int i=0; i<len; i++) {
            s[i]-=a.s[i];
            if(s[i]<0) {
                s[i]+=10;
                s[i+1]--;
            }
        }
        while(len>1&&s[len-1]==0) len--;
    }
    void multiply(int b) { //高精度乘单精度
        int temp=0;
        for(int i=0; i<len; i++) {
            temp+=s[i]*b;
            s[i]=temp%10;
            temp/=10;
        }
        s[len++]=temp;
        while(s[len-1]>10) {
            s[len]+=s[len-1]/10;
            s[len-1]%=10;
            len++;
        }
        while(len>1&&s[len-1]==0) len--;
    }
    void multiply(Hp b) { //高精度乘高精度
        Hp c;
        mt(c.s,0);
        for(int i=0; i<len; i++) {
            for(int j=0; j<b.len; j++) {
                c.s[i+j]+=s[i]*b.s[j];
                c.s[i+j+1]+=c.s[i+j]/10;
                c.s[i+j]%=10;
            }
        }
        len=len+b.len;
        for(int i=0; i<len; i++) {
            s[i]=c.s[i];
        }
        while(len>1&&s[len-1]==0) len--;
    }
} A,B;
char a[M];
int main(){
    int t;
    while(~scanf("%d",&t)){
        int cas=1;
        while(t--){
            scanf("%s",a);
            A.init(a);
            B.init(a);
            A.multiply(B);
            A.multiply(8);
            B.multiply(7);
            A.subtract(B);
            A.add(1);
            printf("Case #%d: ",cas++);
            A.print();
            puts("");
        }
    }
    return 0;
}

 那个类似的题目被我找到了

折线分割平面 http://acm.hdu.edu.cn/showproblem.php?pid=2050

一样的解一元二次

#include<cstdio>
int main(){
    int t,n;
    while(~scanf("%d",&t)){
        while(t--){
            scanf("%d",&n);
            printf("%d
",2*n*n-n+1);
        }
    }
    return 0;
}

end