The 9th Zhejiang Provincial Collegiate Programming Contest

终于决定做套题了，这个是个开始！

首先给出两个解题报告网址，http://blog.watashi.ws/2274/zjpcpc2012/watashi大神的，另一个http://blog.sina.com.cn/s/blog_5123df350100zp0s.html

题目在zoj 3600 - 3610

zoj 3600 水题，不解释

zoj 3601 Unrequited Love 这个题还没做，不会

zoj 3602 Count the Trees 

           首先是树的区别，用的是map< pair <int,int> ,int >这个题最开始的时候不会，然后看到这个的时候就明白了，怎么区别不同的子树，区别开之后就是dfs的过程了，但是遇到了第二个难题，这个dfs居然会爆栈，然后就纠结怎么写模拟栈，考虑复杂了，写了好长时间，哎，看来水平还是不行啊！

#include <iostream>
#include <stdio.h>
#include <map>
#include <utility>
#include <string.h>
using namespace std;
#define N 110000
long long sum;
struct note
{
    int left,right;
}tree[2][N];

map < pair< int ,int > ,int > mymap;int cs;

pair <int ,int> p;
struct sta
{
    int root,type;
}mystack[N];
int indexx[2][N];
long long count[2][N*2];
int top;

void dfs(int kind,int root)
{
    mystack[0].root = root;
    mystack[0].type = 0;
    top = 0;
    while(top >= 0)
    {
        if(mystack[top].type == 0)
        {
            mystack[top].type++;
            if(tree[kind][mystack[top].root].left != -1)
            {
                top++;
                mystack[top].type = 0;
                mystack[top].root = tree[kind][mystack[top-1].root].left;
            }
        }else if(mystack[top].type == 1)
        {
            mystack[top].type++;
            if(tree[kind][mystack[top].root].right != -1)
            {
                top++;
                mystack[top].type = 0;
                mystack[top].root = tree[kind][mystack[top-1].root].right;
            }
        }else if(mystack[top].type == 2)
        {
            int l,r;
            if(tree[kind][mystack[top].root].left == -1) l = 0;
            else l = indexx[kind][tree[kind][mystack[top].root].left];

            if(tree[kind][mystack[top].root].right == -1) r = 0;
            else r = indexx[kind][tree[kind][mystack[top].root].right];

            p = make_pair(l,r);
            if(mymap [ p ] == 0) mymap[ p ] = cs++;
            count[kind][ mymap[ p ] ]++;

            indexx[kind][mystack[top].root] = mymap[p];
            top--;
        }
    }
}

int main()
{
    int t;
    cin >> t;
    while(t--)
    {
        int n,m;
        cin >> n >> m;
        for(int i = 1;i <= n;i++) scanf("%d%d",&tree[0][i].left,&tree[0][i].right);
        for(int i = 1;i <= m;i++) scanf("%d%d",&tree[1][i].left,&tree[1][i].right);

        cs = 1;
        sum = 0;
        mymap.clear();
        memset(count,0,sizeof(count));
        dfs(0,1);
        dfs(1,1);
        for(int i = 1;i < cs;i++)
        {
            //cout << count[0][i] << " " << count[1][i] << "\n";
            sum += count[0][i]*count[1][i];
        }
        cout << sum << "\n";
    }
}

zoj 3603 Draw Something Cheath 这个题也是简单题，不解释，

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
using namespace std;

int count[100];
int count2[100];

char str[120];

int main()
{
    int n;
    cin >> n;
    while(n--)
    {
        int m;
        cin >> m;
        for(int i = 0;i < 100;i++) count[i] = 100;
        for(int k = 0;k < m;k++)
        {
            scanf("%s",str);
            for(int i = 0;i < 100;i++) count2[i] = 0;

            for(int i = 0;i < strlen(str);i++) count2[str[i] - 'A'] ++;

            for(int i = 0;i < 100;i++) count[i] = min(count[i],count2[i]);
        }
        for(int i = 0;i < 100;i++)
        {
            if(count[i] != 0) while(count[i]--) cout << (char)('A'+i);
        }
        cout << "\n";
    }
    return 0;
}

zoj 3604 Tunnel Network 是一个难题吧，可以看最开始给出的第二个网址，我要说的是，这个和生成树问题很像，重要的是理解 Prüfer 编码与 Cayley 公式（

http://www.matrix67.com/blog/archives/682

），在这个基础上再做，题解中给出的方法其实就是Prufer编码，ok，其实代码还是很简单的，

#include <iostream>
#include <stdio.h>
#include <string.h>
using namespace std;

#define P 1000000007

long long Fpow(int a,int n)
{
    long long ret=1;
    long long A=a;
    while(n)
    {
        if (n & 1)
            ret=(ret*A)%P;
        A=(A*A)%P;
        n>>=1;
    }
    return ret;
}

int main()
{
    int t;
    cin >> t;
    while(t--)
    {
        int n,s;
        cin >> n >> s;
        long long re;
        if(n > s)
            re = s*Fpow(n,n-s-1)%P;
        else re = 1;
        cout << re << "\n";
    }
    return 0;
}

zoj 3605 Find the Marble 这个是个dp?? 和同学讨论的时候，突然提到了一个词，路径数！然后整个题都是围绕路径数展开的，当然最后好像与路径数有一点的差别。

#include <iostream>
#include <stdio.h>
#include <string.h>
using namespace std;

long long dp[55][55][55];//前i次选择j次，s到k的路径数

long long s,n,k,m;

int main()
{
    int t;
    cin >> t;
    while(t--)
    {
        memset(dp,0,sizeof(dp));
        int a,b;
        cin >> n >> m >> k >> s;
        dp[0][0][s-1] = 1;
        for(int i = 1;i <= m;i++)
        {
            cin >> a >> b;
            a--;b--;
            for(int j = 0;j <= k;j++)
            {
                for(int kk = 0;kk < n;kk++)
                {
                    dp[i][j][kk] = dp[i-1][j][kk];
                    if(j > 0 && kk != a && kk != b) dp[i][j][kk] += dp[i-1][j-1][kk];
                }
            }

            for(int j = 1;j <= k;j++)
            {
                dp[i][j][a] += dp[i-1][j-1][b];
                dp[i][j][b] += dp[i-1][j-1][a];
            }
        }
        int re = s-1;
        long long ma = 0;
        for(int i = 0;i< n;i++)
        {
            //cout << dp[m][k][i] << "\n";
            if(dp[m][k][i] > ma)
            {
                ma = dp[m][k][i];
                re = i;
            }
        }
        cout << re+1 << "\n";
    }
    return 0;
}

zoj 3606  Lazy Salesgirl 好一个线段树，想了好久，没想出来，然后看解题报告，数遍，直接贴代码，不懂得看文章最开始给出的网址

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
#define N 100010

#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
struct dat
{
    int pi,ti;
    int id,dis;
}pt[N];
struct note
{
    double sum[4];
    int size;
}data[N*4];

void PushUp(int rt)
{
    data[rt].size = data[rt<<1].size + data[rt<<1|1].size;
    for(int i = 0;i < 3;i++)
    {
        int kk = i+data[rt<<1].size%3;
        data[rt].sum[i] = data[rt<<1].sum[i] + data[rt<<1|1].sum[kk%3];
    }
}

void update(int k,int pi,int l,int r,int rt)
{
    if(l == r && r == k)
    {
        data[rt].size = 1;
        for(int i = 0;i < 3;i++)
            data[rt].sum[i] = pi*(i+1);
        return ;
    }
    int m = (l+r)>>1;
    if(k <= m) update(k,pi,lson);
    else update(k,pi,rson);
    PushUp(rt);
}

void init(int n)
{
    for(int i = 1;i <= n*4+10;i++)
    {
        for(int j = 0;j < 3;j++) data[i].sum[j] = 0;
        data[i].size = 0;
    }
}

bool cmp1(const dat a,const dat b)
{
    return a.ti < b.ti;
}
bool cmp2(const dat a,const dat b)
{
    return a.dis < b.dis || (a.dis == b.dis && a.id < b.id);
}
int main()
{
    int t;
    cin >> t;
    while(t--)
    {
        int n;
        cin >> n;
        for(int i = 1;i <= n;i++) scanf("%d",&pt[i].pi);
        for(int i = 1;i <= n;i++) scanf("%d",&pt[i].ti);

        sort(pt+1,pt+n+1,cmp1);
        for(int i = 1;i <= n;i++) pt[i].id = i;
        pt[0].ti = 0;for(int i = 1;i <= n;i++) pt[i].dis = pt[i].ti - pt[i-1].ti;
        sort(pt+1,pt+n+1,cmp2);

        double ans = 0;double re = 0;

        init(n);

        int now = 1;
        while(now <= n)
        {
            update(pt[now].id,pt[now].pi,1,n,1) ,now++;
            while(now <= n && pt[now].dis == pt[now-1].dis) update(pt[now].id,pt[now].pi,1,n,1),now++;

            //cout << data[1].size << "\n";
            //double tt = data[1].sum[0]/(num);
            double tt = data[1].sum[0]/(data[1].size);
            //µÃµ½tt;
            if(tt > ans)
            {
                re = pt[now-1].dis;
                ans = tt;
            }
        }
        printf("%.6lf %.6lf\n",re,ans);
        //cout << re << " " << ans << "\n";
    }
    return 0;
}

zoj 3607 Lazier Salesgirl 这个是zoj 3606的简化版，这个直接枚举就是了，不算太难吧

#include <iostream>
#include <stdio.h>
#include <algorithm>
using namespace std;

int data[2000],p[2000];

int dis[2000];
int dis2[2000];

int main()
{
    int t;
    cin >> t;
    while(t--)
    {
        int n;
        cin >> n;
        for(int i = 1;i <= n;i++) cin >> p[i];
        for(int i = 1;i <= n;i++) cin >> data[i];

        data[0] = 0;

        for(int i = 1;i <= n;i++) dis[i] = data[i] - data[i-1];

        for(int i = 1;i <= n;i++) dis2[i] = dis[i];

        sort(dis2+1,dis2+n+1);

        double re = 0;double ma = -1;

        for(int i = 1;i <= n;i++)
        {
            int con = 0;double sum = 0;
            for(int j = 1;j <= n;j++)
            {
                if(dis[j] > dis2[i]) break;
                con++;
                sum += p[j];
            }
            if(con != 0 && sum / con > ma) {ma = sum/con;re = dis2[i];}
        }
        if(ma == -1) {re = 0;ma = 0;}
        printf("%.6lf %.6lf\n",re,ma);
       // cout << re << " " << ma << "\n";
    }
    return 0;
}

zoj 3608` Signal Detection 计算几何问题，完全不懂，跳过

zoj 3609 Modular Inverse 坑die题，“ZOJ3609 Modular Inverse属于坑死你简单题。其实这些题就是拿来做ranklist的区分度用的，也让实力较弱的选手不会觉得整场太无聊。”出题者是这么说的，这个最大的坑是ax≡1 (mod m).  这儿是同余的意思，不要枚举x时用ax%m = 1;

而是ax%m = 1%m;

#include <iostream>
#include <stdio.h>
using namespace std;

int main()
{
    int a,m;
    int n;
    cin >> n;
    while(n--)
    {
        cin >> a >> m;
        int re = -1;
        for(int i = 1;i <= m;i++)
        {
            if(i*a%m == 1%m) {re = i;break;}
        }
        if(re == -1) cout << "Not Exist";
        else cout << re;
        cout << "\n";
    }
    return 0;
}

zoj 3610 Yet Another Story of Rock-paper-scissors 逗你玩？？哎，最近比赛“逗你玩”的题目越来越多了！很囧啊

#include <iostream>
#include <stdio.h>
using namespace std;
char b[1000],g[1000],x[100];
int main()
{
    int n;
    cin >> n;
    while(n--)
    {
        scanf("%s%s%s",b,g,x);
        cout << g << " will survive\n";
    }
    return 0;
}

，