ACM-ICPC 2017 Asia Shenyang Solution

A: BBP Formula

https://www.cnblogs.com/LzyRapx/p/7802790.html

#include <bits/stdc++.h>
using namespace std;

#define ll long long

inline ll qpow(ll x, ll n, ll mod)
{
    ll base = x;
    ll ans = 1;
    while (n)
    {
        if (n & 1) ans = (ans * base) % mod;
        base = base * base % mod;
        n >>= 1;
    }
    return ans;
}

inline double BBP(int n, ll k, ll b)
{
    double res = 0.0;
    for (int i = 0; i <= n; ++i)
        res += qpow(16, n - i, (8 * i + b)) * 1.0 / (8 * i + b);
    for (int i = n + 1; i <= (n + 500); ++i)
        res += powf(16, n - i) * 1.0 / (8 * i + b); 
    return k * res; 
}

inline char print(double tmp)
{
    int x = int(tmp);
    if (x >= 0 && x <= 9)
        return x + '0';
    return x - 10 + 'A';
}

int t, n;

inline void Run()
{
    scanf("%d", &t);
    for (int kase = 1; kase <= t; ++kase)
    {
        scanf("%d", &n); --n;
        double ans = BBP(n, 4, 1) - BBP(n, 2, 4) - BBP(n, 1, 5) - BBP(n, 1, 6);
        ans = ans - (int)ans;
        if (ans < 0) ans += 1.0; 
        ans *= 16.0;
        printf("Case #%d: %d %c
", kase, n + 1, print(ans));
    }
}

int main()
{
    #ifdef LOCAL
        freopen("Test.in", "r", stdin);
    #endif

    Run();

    return 0;
}

B: Bridge

留坑。

C: Empty Convex Polygons

留坑。

D: Defense of the Ancients

留坑。

E: Five-roune Show Hand

留坑。

F:Heron and His Triangle

题意：给出一个n，找出一个最小的t满足 t >= n 并且 t-1 t t + 1 三条边组成的三角形的面积的整数

思路：根据海伦公式，然后打表，找出前几项是

4

14

52

194

724

然后发现 F(n) = 4F(n - 1) - F(n - 2)

可以发现，这个数增长的很快，60多项就会超过10^30了 用JAVA打个表，然后暴力找或者二分找都可以

import java.math.BigInteger;
import java.util.Scanner;

public class Main
{

    public static void main(String[] args)
    {
        Scanner in = new Scanner(System.in);
        BigInteger ans[] = new BigInteger[200 + 10];
        ans[0] = BigInteger.valueOf(4);
        ans[1] = BigInteger.valueOf(14);
        for(int i = 2; i <= 200; ++i)
        {
            ans[i] = ans[i - 1].multiply(BigInteger.valueOf(4));
            ans[i] = ans[i].subtract(ans[i - 2]);
        }
        int T = in.nextInt();
        for(int cas = 1; cas <= T; ++cas)
        {
            BigInteger n = in.nextBigInteger();
            for(int i = 0; i <= 200; ++i)
            {
                if(ans[i].compareTo(n) >= 0)
                {
                    System.out.println(ans[i]);
                    break;
                }
            }
        }
    }
}

G: Infinite Fraction Path

题意：给出长度为n的字符串，对于第i位的数，它和第(i^2 + 1) % n位的数有一条有向边，从一个数出发，走n - 1 步，得到一个长度为n的字符串，输出字典树最大的那个

思路：考虑搜索

贪心的想法肯定是第一步是字母序最大的那个字母开始

两条剪枝：

第一条：每次扩展一步，如果有一个点扩展出来到这一位的字母小于其它点过来的 剪掉

第二条：每次扩展一步，如果下标被相同步数到达这里的访问过，剪掉

#include<bits/stdc++.h>

using namespace std;

#define N 200010
typedef long long ll;

int n;
char str[N];
char ans[N];
int used[N];

struct node{
    ll idx;
    int step;
    inline node(){}
    inline node(ll idx, int step) :idx(idx), step(step){}
};

queue<node>q;

inline void Init()
{
    while(!q.empty()) q.pop();
    memset(ans, 0, sizeof ans);
    memset(used, 0, sizeof used);
}

inline void BFS()
{
    node st, now;
    while(!q.empty())
    {
        st = q.front();
        q.pop();
        if(st.step > n - 1) continue;
        if(str[st.idx] < ans[st.step]) continue;
        if(str[st.idx] > ans[st.step])
            ans[st.step] = str[st.idx];
        now.idx = (st.idx * st.idx + 1) % n;
        now.step = st.step + 1;
        if(str[now.idx] < ans[now.step]) continue;
        if(used[now.idx] == now.step) continue;
        used[now.idx] = now.step;
        q.push(now);
        if(str[now.idx] > ans[now.step])
            ans[now.step] = str[now.idx];
    }
}

int main()
{
    int t;
    scanf("%d",&t);
    for(int cas = 1; cas <= t; ++cas)
    {
        Init();
        scanf("%d", &n);
        scanf("%s", str);
        char Max = 0;
        for(int i = 0; i < n; ++i)
        {
            Max = max(Max, str[i]);
        }    
        for(int i = 0; i < n; ++i)
        {
            if(str[i] == Max)
                q.push(node(i, 0));
        }
        BFS();
        printf("Case #%d: ", cas);
        for(int i = 0; i < n; ++i)
        {
            printf("%c",ans[i]);
        }
        printf("
");
    }
    return 0;
}

H:Legends of the Three Kingdoms

留坑。

I:Little Boses

ull的范围是2 ^ 64 - 1

#include <bits/stdc++.h>

using namespace std;
#define ull unsigned long long

int t;
ull a[4];

inline bool check()
{
    ull D = 1ull << 62;
    for (int i = 0; i < 4; ++i)
        if (a[i] != D)
            return false;
    return true;
}

int main()
{
    scanf("%d", &t);
    while (t--)
    {
        for (int i = 0; i < 4; ++i) scanf("%llu", a + i);
        if (check())
        {
            puts("18446744073709551616");
            continue;
        }
        else
        {
            ull sum = 0;
            for (int i = 0; i < 4; ++i) sum += a[i];
            printf("%llu
", sum);
        }
    }
    return 0;
}

J:New Self-describing Sequence

留坑。

K:Rabbits

题意：有n只小兔子，如果两个小兔子中间有间隙，那么在这两个小兔子之外的兔子可以跳到某一个间隙中，求最多能跳多少步

思路：显然 答案是 最左边那个小兔子到最右边的左边的小兔子的间隙数和最右边小兔子到最左边的右边的小兔子的间隙取max

#include <bits/stdc++.h>

using namespace std;

#define N 510

int t;
int n;
int arr[N];

int main()
{
    scanf("%d", &t);
    while (t--)
    {
        scanf("%d", &n);
        for (int i = 1; i <= n; ++i)
            scanf("%d", arr + i);
        int ans = 0;
        ans = arr[n - 1] - arr[1] + - n + 2;
        ans = max(ans, arr[n] - arr[2] - n + 2);
        printf("%d
", ans);
    }
    return 0;
}

L:Tree

题意：给出一棵树，k种颜色，每种颜色的边集是选取最少的边使得这些边覆盖所有这种颜色的点，求k中颜色的边集并

思路：显然，如果一条边存在于所有颜色的边集当中，那么它连接的两边的端点个数一定>= k

#include <bits/stdc++.h>

using namespace std;

#define N 200010

struct node
{
    int to, nx;
    inline node() {}
    inline node(int to, int nx) : to(to), nx(nx) {}
}edge[N << 1];

int head[N], pos;
int used[N];
int sum[N];

inline void Init()
{
    memset(head, -1, sizeof head);
    memset(used, 0, sizeof used);
    memset(sum, 0, sizeof sum);
    pos = 0;
}

inline void addedge(int u, int v)
{
    edge[++pos] = node(v, head[u]); head[u] = pos;
    edge[++pos] = node(u, head[v]); head[v] = pos;
}

int ans;
int t, n, k;

inline void DFS(int u)
{
    used[u] = 1;
    sum[u] = 1;
    for (int it = head[u]; ~it; it = edge[it].nx)
    {
        int v = edge[it].to;
        if (used[v] == 1) continue;
        DFS(v);
        sum[u] += sum[v];
    }
    if (sum[u] >= k && (n - sum[u] >= k)) ++ans;
}

int main()
{
    scanf("%d", &t);
    while (t--)
    {
        scanf("%d%d", &n, &k);
        Init();
        for (int i = 1, u, v; i < n; ++i)
        {
            scanf("%d%d", &u, &v);
            addedge(u, v);
        }
        ans = 0;
        DFS(1); 
        printf("%d
", ans);
        
    }
    return 0;
}

M:Wandering Robots

题意：给出n * n 的矩形，有k个障碍物，机器人刚开始在(0, 0) 点，它会等概率的选择停留在原地，或者走向相邻的可走的格子（即没有障碍物的），时间过去了很久，求机器人在(x, y) (x + y >= n - 1) 的点的概率

思路：LTS大神推出答案就是 (满足条件的点的可能性) / (所有点的可能性)

先公式求出所有点的可能性

然后枚举障碍物，减去多加的可能性

首先障碍物本身要减去自身的可能性，如果它的相邻点不是障碍物，还要减去一

#include<bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef pair <int, int> pii;

inline ll gcd(ll a, ll b)
{
    while(b ^= a ^= b ^= a %= b);
    return a;
}

ll n, k;

map <pii, int> mp;

int Move[][2] = 
{
    0, 1,
    0,-1,
    1, 0,
   -1, 0,
};

inline bool ok(int x, int y)
{
    if (x < 0 || x >= n || y < 0 || y >= n) return false;
    return true;
}

int main()
{
    int t;
    scanf("%d",&t);
    for(int cas = 1; cas <= t; ++cas)
    {
        scanf("%lld %lld", &n, &k);
        mp.clear();
        for (int i = 1, x, y; i <= k; ++i)
        {
            scanf("%d%d", &x, &y);
            mp[pii(x, y)] = 1;
        }
        printf("Case #%d: ", cas);
        if (n == 1)
        {
            if (k == 0)
                puts("1/1");
            else 
                puts("0/0");
        }
        else
        {
            ll fenmu = (n - 2) * (n - 2) * 5 + 16 * (n - 2) + 12;
            ll fenzi = (n - 2) * 8 + 5 * (n - 1) * (n - 2) / 2 + 9;
            for (map<pii, int>::iterator it = mp.begin(); it != mp.end(); ++it)
            {
                int x = it->first.first, y = it->first.second;
                int tmp = 5;
                if (x == 0 || x == n - 1)
                    tmp--;
                if (y == 0 || y == n - 1)
                    tmp--;
                fenmu -= tmp;
                if (x + y >= n - 1)
                    fenzi -= tmp;
                for (int i = 0; i < 4; ++i)
                {
                    int dx = x + Move[i][0];
                    int dy = y + Move[i][1];
                    if (!ok(dx, dy)) continue;
                    if (mp.count(pii(dx, dy)) == 1) continue;
                    --fenmu;
                    if (dx + dy >= n - 1) 
                        --fenzi;
                }
            }
            ll G = gcd(fenzi, fenmu);
            fenzi /= G, fenmu /= G;
            printf("%lld/%lld
", fenzi, fenmu);
        }
    }
    return 0;
}