2019 ICPC Asia Nanchang Regional (L, C, E, G)

2019 ICPC Asia Nanchang Regional (计算客)

ICPC-Nanchang 终榜

---

L
签到。
C
签到题。时隔久远已经不记得了到底是个啥了，尽管当时过了还是当作练习题再做一下。(并且应用Vim编辑一下代码）
如果给定$i$要求$j$，共有 $2^{i最高位的1之后的所有0的个数}$ 种可能的组合情况。
如果给定$n$要求$i$，可以枚举$i$ 在$n$的哪一位$1$上出现第一个$1$，之后的1通过组合数计算出现不同次数的0的个数对于$j$的影响即可。
需要预处理$n$的某位数之后出现的0的个数、2的幂次和3的幂次。

const int maxn = 1e5 + 10;
const int M = 1e9 + 7;
char s[maxn];
int b[maxn];
ll pow2[maxn], pow3[maxn];
int quick_pow(int a, int p){
	ll res = a, ans = 1;
	while(p){
		if(p & 1) ans = (ans * res) % M;
		res = (res * res) % M;
		p >>= 1;
	}
	return ans % M;
}
	
void ini(){
	pow3[0] = 1;
	for(int i = 1; i < 100000; i++){
		pow3[i] = pow3[i - 1] * 3 % M;
	}
	pow2[0] = 1;
	for(int i = 1; i < 100000; i++){
		pow2[i] = pow2[i - 1] * 2 % M;
	}
}
int main(){
//    Fast;
	ini();
	int t; scanf("%d", &t);
	while(t--){
		scanf("%s", s); int top =  (int)strlen(s);
		int k = 0;
		for(int i = top - 1; i >= 0; i--){
            b[i] = b[i + 1] + (s[i] == '0');
            if(s[i] == '1') k++;
		}
		ll ans = 0;
		for(int i = 0; i < top; i++){
			if(s[i] == '1') ans += (pow2[b[i]] * pow3[--k]) % M;
			ans %= M;
		}
		printf("%lld\n", ans + 1);
	}
    return 0;
}

E
签到题。因为只对白边的选取数量加以限制，所以可以将黑边全部选上，然后再用白边按照权值从大到小构造一次MST，如果构造完成后仍然可以再选取若干条边就完全按照权值选取白边。
这里有一个小小的存图问题，如果用前向星存两条有向边作为一条无向边，在排序后不能保证那些权值相同的不同边之间的顺序关系，因而不能较好地记录某条边是否被访问过的问题。为此可以直接用 {起点、终点、权值}的结构体来存边，这样子也是一种可行的图存储方式。
WA了好几发的原因： vector用爽了结果对应数组的大小开小了；本地重定向输入流忘关了
计蒜客的这个数据有点弱…我的假数据范围都能AC…感觉不是场上真正的测试数据…

const int maxn = 5e4 + 5;
const int maxe = 5e5 + 5;
struct star{
    int u, v, w;
};
vector<star> vec;
int n, m, k;
int parent[maxn];
int find(int x){
    return x == parent[x]? x: parent[x] = find(parent[x]);
}
void merge(int x, int y){
    int r1 = find(x), r2 = find(y);
    if(r1 != r2) parent[r2] = r1;
}
int vis[maxe];
int isok(){
    for(int i = 2; i <= n; i++)
        if(find(i) != find(i - 1)) return 0;
    return 1;
}
bool cmp(star a, star b){
    return a.w > b.w;
}
void ini(){
    for(int i = 0; i <= n; i++) parent[i] = i;
    vec.clear();
    memset(vis, 0, sizeof vis);
}
int main(){
//    Fast;
    int t; scanf("%d", &t);
    while(t--){
        scanf("%d %d %d", &n, &m, &k);
        ini();
        ll ans = 0;
        for(int i = 0; i < m; i++){
            int u, v, w, op;
            scanf("%d %d %d %d", &u, &v, &w, &op);
            if(!op){ans += w; merge(u, v);}
            else vec.push_back({u, v, w});
        }
        sort(vec.begin(), vec.end(), cmp);
        int top = (int)vec.size();
        for(int i = 0; i < top && k; i++){
            int u = vec[i].u, v = vec[i].v;
            if(find(u) != find(v)){
                vis[i] = 1; k--;
                merge(u, v);
                ans += vec[i].w;
            }
        }
        for(int i = 0; i < top && k; i++){
            if(vis[i]) continue;
            ans += vec[i].w;
            k--;
        }
        if(isok()) printf("%lld\n", ans);
        else puts("-1");
    }
    return 0;
}

---

补题
G
在场上看这道题的时候，第一反应就是先对那个奇怪的大模数分解下因子，然后分出来了后却没什么进一步的想法…此外队友对阶乘打表后并没有发现2803之后均为0…最后我本着不让机子空的原则又去瞎捣鼓打表，最终发现了2803这个数字…后来2个小时硬是没写出来…
菜是原罪…

不过其实再看这道题我还是感觉雾里看花花不语…
正确做法是：根据2803将数进行划分，枚举每条两个端点均不为0的线段，记录其吃的食物数量和获得的食物重量，用前缀和优化时间复杂度。之后离线所有的询问，按照询问的$k$大小进行排序，双指针跑一边$dis$数组和$query$数组即可。

坑点
在求$dis$数组的时候，因为是在取模下的前缀和减法，所以一定要记得对减法进行正数取模！

代码

const int M = 998857459;
const int maxn = 1e5 + 10;

int top;
ll f[3000], a[maxn], sum[maxn];
ll ans[maxn];

vector<int> vec;
ll dis[maxn];

struct node{int id, q;};
bool cmp(node a, node b){return a.q < b.q;}
vector<node> v;

void ini(){
    f[0] = 1;
    for(int i = 1; i < 3000; i++) f[i] = (f[i - 1] * i) % M;
    memset(ans, -1, sizeof ans);
}

int main(){
//    Fast;
    ini();
    int n, query; scanf("%d %d", &n, &query);
    for(int i = 1, temp; i <= n; i++){
        scanf("%lld", a + i);
        if(a[i] < 2803){
            vec.push_back(i);
            temp = (int)f[a[i]];
        }
        else temp = 0;
        sum[i] = (sum[i - 1] + temp) % M;
    }
    top = (int)vec.size();
    for(int i = 0; i < top; i++) for(int j = i; j < top; j++){
        dis[vec[j] - vec[i] + 1] = max(dis[vec[j] - vec[i] + 1], (sum[vec[j]] - sum[vec[i] - 1] % M + M) % M);
    }
    for(int i = 0, temp; i < query; i++){
        scanf("%d", &temp);
        v.push_back({i, temp});
    }
    sort(v.begin(), v.end(), cmp);
    int p = 0;
    for(itn i = 1; i <= n; i++){
        while(p < query && dis[i] >= v[p].q){
            ans[v[p].id] = i;
            p++;
        }
    }
    for(int i = 0; i < query; i++) printf("%lld\n", ans[i]);
    return 0;
}

总而言之，这道题的综合算法自己掌握的还是很不扎实…其中求$dis$数组的想法和后面离线查询等算法思想值得学习！

dis数组好像在尚未发现2803即完全考虑1e5时讨论到，因为无法接受的1e5所以拒绝了该想法？

---

参考博客：
2019 ICPC Asia Nanchang Regional（C,E,G,K,L）