2012 MUTC 7 总结

题解链接：http://page.renren.com/601081183/note/865867924

hdu题号：4360~4369

　　今天状态不太好。中午，我带着疲惫的身躯去到机房，但是开门的人迟迟没来，对于我们的比赛开始时间也就延迟了几分钟。今天也是很久都没有队伍带榜，不过我一打开题目就看懂了两个，一个是1006，另外的是1001。

　　比赛刚开始的时候，我1001还没更新题意，所以没敢试这题，不过初步可以知道是最短路。然后，我观察了一下1006的输入输出，我很快就发现这个对称图形只需统计其中一部份的可操作点的个数，然后快速幂一下，答案就出来了。我第一遍打出来的程序出现了一个问题，就是点的数目居然数多了。于是我又想了一下，画一下图，很快我又发现了，这是要统计1/8的正方形，而不是1/4,。于是我就赶紧改了，测试数据通过以后就立马扔了上去，不过返回了一个TLE。然后，我就很郁闷的想了一下，不能偷懒直接用bool数组来储存点的状态，而是应该用一个点的结构体来储存，在统计前排一下序，然后挑出非重复点，最后用1/8块中的点数来减去点的个数，这样就得到可操作点的个数了。

　　代码最快速度只需15ms，差距实在是大啊！

1006（hdu 4365）的代码：

#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>

using namespace std;
typedef __int64 ll;
const int maxn = 5005;
const ll mod = 100000007;

void scan(int &k)
{
    char c;
    while ((c = getchar()) && (c > '9' || c < '0'));
        k = c - '0';
    while((c = getchar()) && (c <= '9' && c >= '0'))
        k = k * 10 + c - '0';
}

struct point
{
    int x;
    int y;

    bool operator<(const point & a) const
    {
        if (a.x != x) return x < a.x;
        return y < a.y;
    }
} p[maxn];

int cal(int n, int m)
{
    int cnt = 0;
    int t;

    t = n / 2;
    if (n & 1) t++;
    if (!m) return t * (t + 1) / 2;
    sort(p + 1, p + m + 1);
    cnt++;
    for (int i = 2; i <= m; i++)
    {
        if (p[i - 1].x == p[i].x && p[i - 1].y == p[i].y) continue;
        cnt++;
    }

    /*
        if (n & 1){
            t = n / 2;
            for (int i = 1; i <= t + 1; i++){
                for (int j = i; j <= t + 1; j++){
                    if (!mp[i][j]) cnt++;
                }
            }
        }
        else{
            t = n / 2;
            for (int i = 1; i <= t; i++){
                for (int j = i; j <= t; j++){
                    if (!mp[i][j]) cnt++;
                }
            }
        }
     */
#ifndef ONLINE_JUDGE
    printf("t  %d\n", t);
#endif
    return t * (t + 1) / 2 - cnt;
}

void con(int n, int &x, int &y)
{
    int t = n >> 1;

    if (n & 1)
    {
        if (x > t) x = n - x + 1;
        if (y > t + 1) y = n - y + 1;
    } else
    {
        if (x > t) x = n - x + 1;
        if (y > t) y = n - y + 1;
    }
    if (x > y)
    {
        t = x;
        x = y;
        y = t;
    }
}

ll work(ll k, int m)
{
    ll ret = 1;

    while (m){
        if (m & 1) ret *= k, ret %= mod;
        k *= k;
        k %= mod;
        m >>= 1;
    }

    return ret;
}

bool deal()
{
    int n, m, k;
    int x, y;

    if (scanf("%d%d%d", &n, &m, &k) == EOF) return false;
    /*
    for (int i = 1, end = (n >> 1) + 1; i <= end; i++){
        for (int j = 1; j <= end; j++){
            mp[i][j] = false;
        }
    }
     */
    for (int i = 1; i <= m; i++)
    {
        scan(x);
        scan(y);
        x++;
        y++;
        con(n, x, y);
#ifndef ONLINE_JUDGE
        printf("con  x %d  y %d\n", x, y);
#endif
        //mp[x][y] = true;
        p[i].x = x;
        p[i].y = y;
    }

    m = cal((ll)n, m);
    printf("%I64d\n", work(k, m));

    return true;
}

int main()
{
#ifndef ONLINE_JUDGE
    freopen("in", "r", stdin);
#endif
    while (deal());

    return 0;
}

　　我出完这题以后，我的队友已经将1003的转移方程写好了，于是我就立马让位置给他，让他赶紧打上去。dp的题代码 长度通常都不大，所以他很快就打好了。可惜的是，难得让他完成一题，结果这题是要用单调队列优化的，所以结果是TLE。然后我们花了好长时间想优化的方法，可惜还是没想到。赛后，师兄告诉我们用单调队列两边各扫一次，得到两个单调队列，然后，对于递增和递减的区间分别用两种不同的队列更新。据说还可以用二分搜索来加速，思路我已经有了，代码明天再打吧！

　　然后我就尝试这做1001，我把输入输出都打好了，就连拆点的操作也想好了，不过当时状态真实相当不好，忘记了spfa的写法，搞到搞半天yy了个dij+heap，最后还要写烂了，最终只好放弃。回来宿舍后，打了一遍，sample轻松通过，不过就是没想到这里可以有自环，还可以有等多个自环原地转啊转的！所以在我查看数据之前都没有意识到这个问题。当然，改过来以后就ac了。

1001（hdu 4360）的代码：

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <queue>
#include <algorithm>
#include <map>
#define ONLINE_JUDGE 1
using namespace std;

typedef __int64 ll;

const int maxv = 1500;
const int maxe = 15000;
const ll inf = 1000000000000;

map<char, int> lt;
struct edge{
    int s, t;
    ll l;
    bool operator < (const edge &x) const{
        return s < x.s;
    }
}E[maxe << 1];

ll dis[maxv << 2];
int cnt[maxv << 2], me[maxe << 1];
int q[maxv << 2], qh, qt;
bool inq[maxv << 2];

void init(){
    lt['L'] = 0;
    lt['O'] = 1;
    lt['V'] = 2;
    lt['E'] = 3;
}

void spfa(int v, int e, int s){
    qh = qt = 0;
    sort(E, E + e);

    int p = 0;

    for (int i = 0; i < v; i++){
        inq[i] = false;
        dis[i] = inf;
        cnt[i] = 0;
    }
    for (int i = 0; i < e; i++){
        while (p <= E[i].s) me[p++] = i;
    }
    while (p <= v) me[p++] = e;
#ifndef ONLINE_JUDGE
    for (int i = 0; i < 20; i++){
        printf("me %d : %d    E : %d %d %I64d    dis %I64d\n", i, me[i], E[i].s, E[i].t, E[i].l, dis[i]);
    }
#endif

    for (int i = me[s]; i < me[s + 1]; i++){
        int et = E[i].t;

        if (dis[et] > E[i].l) dis[et] = E[i].l;
        cnt[et] = 1;
        if (!inq[et]){
            q[qt++] = et;
            inq[et] = true;
        }
    }
    while (qh < qt){
        int cur = q[qh++];
#ifndef ONLINE_JUDGE
        printf("pass %d\n", cur);
#endif

        inq[cur] = false;
         for (int i = me[cur]; i < me[cur + 1]; i++){
            int et = E[i].t;

            if (dis[et] > dis[cur] + E[i].l){
                dis[et] = dis[cur] + E[i].l;
                cnt[et] = cnt[cur] + 1;
                if (!inq[et]){
                    q[qt++] = et;
                    inq[et] = true;
                }
            }
            else if (dis[et] == dis[cur] + E[i].l && cnt[et] < cnt[cur] + 1){
                cnt[et] = cnt[cur] + 1;
                if (!inq[et]){
                    q[qt++] = et;
                    inq[et] = true;
                }
            }
        }
    }
}


void deal(int cc){
    int n, m, N;

    scanf("%d%d", &n, &m);
    N = n << 2;
    for (int i = 0; i < m; i++){
        int s, t, len, kd;
        char buf[3];

        scanf("%d%d%d%s", &s, &t, &len, buf);
        s--; t--;
        kd = lt[buf[0]];
        E[i << 1].s = kd * n + s;
        E[i << 1].t = ((kd + 1) % 4) * n + t;
        E[i << 1].l = len;
        E[i << 1 | 1].s = kd * n + t;
        E[i << 1 | 1].t = ((kd + 1) % 4) * n + s;
        E[i << 1 | 1].l = len;
    }

    m <<= 1;
    spfa(N, m, 0);

    if (dis[n - 1] < inf) printf("Case %d: Cute Sangsang, Binbin will come with a donkey after travelling %I64d meters and finding %d LOVE strings at last.\n", cc, dis[n - 1], cnt[n - 1] / 4);
    else printf("Case %d: Binbin you disappoint Sangsang again, damn it!\n", cc);
}

int main(){
    int T;

#ifndef ONLINE_JUDGE
    freopen("in", "r", stdin);
    freopen("out", "w", stdout);
#endif
    init();
    scanf("%d", &T);
    for (int i = 1; i <= T; i++){
        deal(i);
    }

    return 0;
}

　　最后半个钟，我的队友才辛苦看懂1005这原本很水的模拟题。然后就手忙脚乱的打上代码，提交了上去。思路是对的，sample也是对的，不过我一开始的时候没有处理行末的空格和最后的空行，所以直到比赛结束我们都没有通过这题。赛后改了这个问题以后就立马过了.......T^T    这下玩大了，那时还以为没改会先得到PE，结果oj实在是狠，直接返回wa，搞到我们先是觉得是想法错误了，所以更改的方向也错了......

1005（hdu 4364）的代码：

#include <cstdio>

int mt[4][4];
int op[4][4] = {2, 3, 1, 1, 1, 2, 3, 1, 1, 1, 2, 3, 3, 1, 1, 2};

int deal(int OP, int in){
    switch(OP){
        case 1: return in;
        case 2: {
            in <<= 1;
            if (in > 0xff){
                in ^= 0x1b;
            }
            return in;
        }
        case 3:{
            int tmp = in << 1;
            in = tmp ^ in;
            if (in > 0xff){
                in ^= 0x1b;
            }
            return in;
        }
    }
    return 0;
}

int main(){
    int T;
#ifndef ONLINE_JUDGE
    freopen("in", "r", stdin);
#endif
    scanf("%d", &T);
    while (T--){
        for (int i = 0; i < 4; i++){
            for (int j = 0; j < 4; j++){
                scanf("%X", &mt[i][j]);
#ifndef ONLINE_JUDGE
                printf("%d ", mt[i][j]);
#endif
            }
#ifndef ONLINE_JUDGE
            puts("");
#endif
        }
        int out;

        for (int i = 0; i < 4; i++){
            for (int j = 0; j < 4; j++){
                out = 0;
                for (int k = 0; k < 4; k++){
                    out ^= deal(op[i][k], mt[k][j]);
                }
                if (j < 3) printf("%02X ", out & 0xff);
                else printf("%02X\n", out & 0xff);
            }
        }

        if (T) puts("");
    }

    return 0;
}

　　吸取教训了，以后不能靠oj来判断是否格式错误，而是要根据题意来判断！

　　继续努力！

——written by Lyon

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