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 3604Tunnel 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;
}

，