51nod算法马拉松32

比赛链接

马拉松是真的难…
应该是我太菜了TnT

A 特殊表示法

同bzoj1534
题解戳这里

B 小树

题意：

n个有标号的点，其中m个是叶子节点。
问有多少数的形态。

做法：

引入prufer序列。（不知道的百度）

发现叶子节点不出现在序列里。于是我们钦定m个叶子节点，剩下的点放到序列里，每个点至少出现一次。
然后这个相当于n个球放到m个盒子里，球不同盒子不同，盒子非空。公式是第二类斯特灵数S(n,m)*m!。
关于球盒问题，给一个博客链接，这个比较全。
然后S(n,m)不能递推，要用这个公式

（其实很多东西自己百度百科一下就可以QAQ）

代码：

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;

const int N = 1000010;
const int MOD = 1e9 + 7;
int n, m;
LL fac[N], inv[N];

inline LL ksm(LL x, int p)
{
    LL ret = 1;
    while(p) {
        if(p&1) ret = ret*x%MOD;
        x = x*x%MOD; p >>= 1;
    } return ret;
}
inline void prepare()
{
    fac[0] = 1;
    for(int i = 1; i <= 1000000; i ++) fac[i] = fac[i-1]*i%MOD;
    inv[1000000] = ksm(fac[1000000], MOD-2);//inv[i]表示i阶乘的逆元
    for(int i = 999999; i >= 1; i --) inv[i] = inv[i+1]*(i+1)%MOD;
}
inline LL C(int n, int m){ if(!m || n==m) return 1; return fac[n]*inv[m]%MOD*inv[n-m]%MOD; }
inline LL S(int n, int m)
{
    LL ret = 0;
    for(int k = 0; k <= m; k ++) {
        if(k&1) ret = (ret-C(m, k)*ksm(m-k, n)%MOD+MOD)%MOD;
        else ret = (ret+C(m, k)*ksm(m-k, n)%MOD)%MOD;
    }
    return ret;
} //s(n, m) = 1/m!*sigma(k=0 to m){(-1)^k*C(m,k)*(m-k)^n};
int main()
{
    prepare();
    scanf("%d%d", &n, &m);
    printf("%lld\n", C(n, m)*S(n-2, n-m)%MOD);
    return 0;
}

C 如何愉快地与STL玩耍

题意：

有一个set的序列，有两种操作：
1. 在序列[l,r]内每个位置插入一个数
2. 询问[l,r]内的数去重以后的第k小

做法：

这道题和zjoi2013 k大数查询很像。
那道题是树套树，题解戳这里。
但是这道题要去重，所以不能这么做。

标算是bitset优化线段树，线段树每个节点套一个bitset。
具体操作见代码。
（这还是我第一次写bitset…QwQ）

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cctype>
#include<bitset>
#define mid (l+r>>1)
#define lc o<<1
#define rc o<<1|1
using namespace std;
typedef long long LL;

inline LL read()
{
    char ch = getchar(); LL x = 0; int op = 1;
    for (; !isdigit(ch); ch = getchar()) if (ch == '-') op = -1;
    for (; isdigit(ch); ch = getchar()) x = x*10+ch-'0';
    return op*x;
}
inline void writeAbs(LL a){ if (a >= 10) writeAbs(a/10); putchar('0'+a%10); }
inline void write(LL a){ if (a < 0) putchar('-'); writeAbs(abs(a)); }

const int N = 10001, M = 65540;
int n, q, tot;
bool vis[M<<1];//vis[i]=1表示i号节点被完全覆盖过
bitset<N> T[M<<1], tag[M<<1], tmp;

inline void ask(int o, int l, int r, int x, int y)
{
    if(l == x && r == y) { tmp |= T[o]; return; }
    if(vis[o]) tmp |= tag[o];
    if(y <= mid) ask(lc, l, mid, x, y);
    else if(x > mid) ask(rc, mid+1, r, x, y);
    else ask(lc, l, mid, x, mid), ask(rc, mid+1, r, mid+1, y);
}
inline void insert(int o, int l, int r, int x, int y, int p)
{
    T[o][p] = 1;
    if(l == x && r == y) { vis[o] = 1; tag[o][p] = 1; return; }
    if(y <= mid) insert(lc, l, mid, x, y, p);
    else if(x > mid) insert(rc, mid+1, r, x, y, p);
    else insert(lc, l, mid, x, mid, p), insert(rc, mid+1, r, mid+1, y, p);
}
int main()
{
    n = read(), q = read();
    while(q --) {
        int opt = read(), x = read(), y = read(), w = read();
        if(opt == 1) insert(1, 1, n, x, y, w);
        else {
            tmp.reset(); ask(1, 1, n, x, y);
            if(w == 0 || w > tmp.count()) { puts("-1"); continue; }
            int l = 0, r = 10000;
            while(l < r) {
                int m = l+r>>1, t = ((tmp>>l)<<(l+10000-m)).count();
                if(t >= w) r = m; else { l = m+1; w -= t; }
            }
            write(l), puts("");
        }
    }
    return 0;
}

DEF都还没看过= =先留坑，以后哪天突发奇想做了再写。