ZROI week2

$$ZROI week2$$

除草机

首先考虑最少的拐点肯定是那种螺旋形状的，然后手玩几个数据发现和列数（行数）有关，且每增加1就是上一个状态加2，直接\(O(1)\)公式即可

吐槽：为啥\(n,m\)不给成\(10^{18}\)次方？？

收集隔膜

一眼就能看出双向搜索，首先打了个暴力企图对拍，先交了上去，\(WA?\)，发现和本地输出不一样，搞这个搞了半个小时还是没搞出为什么，然后我用其他的\(IDE\)测发现还是这个毛病，我的双向搜索就没敢交，也不知道改了哪里一个暴力就交过了。

最终得了60分。

• 巧妙的思路：由于贡献是一定的，所以考虑预处理出贡献，一遍\(dfs\)即可。
• 正解的思路：双向搜索，二分查找统计贡献。

翻转序列

想了一会儿想了个线段树的办法，然后疯狂写，一个半小时过去，死了

考虑一个对称中心\(T\)，统计出至少多少个字母对称就会产生一个贡献，显然信息量是\(O(N)\)的，这个统计完了之后可以用线段树区间加和区间\(max\)去求解。

然后就死了

事实证明我的方法还是对的

正解就是把线段树统计换成一些奇奇怪怪的，例如前缀和一类的东西。

树

当我写到这个题已经没时间了（10min），想了一会儿想出了正解，对于每条边建立两个点，链内同向，链间反向，在\(LCA\)处统计贡献即可。

懒得写了

还是敲了一手暴力交了上去，本想再拿点分发现没时间了。

正解其实和我思路差不多，考虑把没有贡献的地方缩成一个点，然后快速幂计算即可。

总结与心得

会的总是拿不上分有点难受，希望下次更加努力。

作业

bzoj 4800

题目大意：有\(n\)个物品，\(m\)块钱，给定每个物品的价格，求买物品的方案数

范围：\(n \leq 40,m \leq 10^{18}\)

题解：

如果直接凑是\(2^{40}\)的复杂度，发现贡献是可以合并的，考虑双向搜索，每一半用\(2^{20}\)去搞，最后合并起来即可。

复杂度：\(O(2^{20} + 2^{20})\)

代码：

#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int MAXN = 1e6 + 10;
ll a[50];
int lim;
int cnt1;
int cnt2;
ll A[MAXN];
ll B[MAXN];
int n;	
ll m;
ll ans;
int mid;

void dfs1(int now,ll sum) {
	if(now > mid) {
		A[++cnt1] = sum;
		return;
	}
	dfs1(now + 1,sum);
	if(sum + a[now] <= m) {
		dfs1(now + 1,sum + a[now]);
	}
}

void dfs2(int now,ll sum) {
	if(now > lim) {
		B[++cnt2] = sum;
		return;
	}
	dfs2(now + 1,sum);
	if(sum + a[now] <= m) {
		dfs2(now + 1,sum + a[now]);
	}
}

signed main () {
	scanf("%d %lld",&n,&m);
	for(int i = 1;i <= n; ++i) {
		scanf("%lld",&a[++lim]);
		if(a[lim] > m) --lim;
	}
	mid = lim >> 1;
	dfs1(1,0);
	dfs2(mid + 1,0);
	sort(A + 1,A + cnt1 + 1);
	sort(B + 1,B + cnt2 + 1);
	for(int i = 1,j = cnt2;i <= cnt1; ++i) {
		while(A[i] + B[j] > m and j) --j;
		if(j == 0) break;
		ans += j;
	}
	cout<<ans<<endl;
	return 0;
}

bzoj 1202

第一眼看到觉得是一个差分约束，第二眼觉得是个并查集。。。

并查集：考虑维护一个到根的路径权值，合并时判断即可。

#include <bits/stdc++.h>
using namespace std;
const int MAXN = 2e5 + 10;
int read () {
	int q=0,f=1;char ch=getchar();
	while(!isdigit(ch)){
		if(ch=='-')f=-1;ch=getchar();
	}
	while(isdigit(ch)){
		q=q*10+ch-'0';ch=getchar();
	}
	return q*f;
}
int f[MAXN];
int a[MAXN];
int find(int x) {
	if(x == f[x]) {
		a[x] = 0;
		return x;
	}
	int res = find(f[x]);
	a[x] = a[x] + a[f[x]];
	return f[x] = res;
}
int n,m;
int tag;
int main () {
//	freopen("1.in","r",stdin);
//	freopen("ans.out","w",stdout);
	int T = read();
	while(T--) {
		tag = 0;
		n = read(),m = read();
		for(int i = 1;i <= n + 1; ++i) {
			f[i] = i;
		}
		memset(a,0,sizeof a);
		while(m--) {
			int s = read(),t = read(),val = read();
			int l = find(s);
			int r = find(t + 1);
			if(l != r) {
				f[l] = r;
				a[l] = a[t + 1] - a[s] - val;
			}
			else if(a[t + 1] - a[s] != val) tag = 1;
		}
		cout<<(tag == 0 ? "true" : "false")<<endl;
	}
	return 0;
}

差分约束：按照输入的方法建一张图，跑可行的\(spfa\)即可。

#include <bits/stdc++.h>
using namespace std;
const int MAXN = 2e5 + 10;
int read () {
	int q=0,f=1;char ch=getchar();
	while(!isdigit(ch)){
		if(ch=='-')f=-1;ch=getchar();
	}
	while(isdigit(ch)){
		q=q*10+ch-'0';ch=getchar();
	}
	return q*f;
}

int n,m;

struct edge {
	int to;
	int nxt;
	int w;
}e[MAXN << 1];
int cnt;
int vis[MAXN << 1];
int inq[MAXN << 1];
int dis[MAXN << 1];
int head[MAXN << 1];
void add(int u,int v,int w) {
	e[++cnt].to = v;
	e[cnt].nxt = head[u];
	head[u] = cnt;
	e[cnt].w = w;
}
int spfa(int now) {
	vis[now] = 1;
	inq[now] = 1;
	for(int i = head[now];i;i=e[i].nxt) {
		int y = e[i].to;
		if(dis[y] > dis[now] + e[i].w) {
			dis[y] = dis[now] + e[i].w;
			if(inq[y] or spfa(y) == 0) return 0;
		}
	}
	inq[now] = 0;
	return 1;
}
int tag;
int T;
int main () {
	//freopen("1.in","r",stdin);
	//freopen("ans.out","w",stdout);
	T = read();
	while(T--) {
		n = read(),m = read();
		cnt = 0;
		memset(inq,0,sizeof inq);
		memset(vis,0,sizeof vis);
		memset(head,0,sizeof head);
		memset(dis,0x3f,sizeof dis);
		for(int i = 1;i <= m; ++i) {
			int x = read(),y = read(),z = read();
			add(x ,y - 1,z);
			add(y - 1, x , z);
		}
		tag = 0;
		for(int i = 0;i <= n; ++i) {
			if(!vis[i]) {
				dis[i] = 0;
				cout<<spfa(i)<<endl;
				if(spfa(i) == 0) {
					puts("false");
					tag = 1;
					break;
				}
			}
		}
		if(!tag) puts("true");
	}
	return 0;
}

bzoj 2679

吐槽：\(map\)定义的\(vis\)数组难道不能标记？

本意就是相当于将一些数划分成两个相等的集合，这里集合相等表示集合的数之和相等，数可以不划分完。

直接处理是\(3^n\)的，为什么？

对于每个数有三种情况：选到左边，选到右边，不选。

直接暴力的话就是\(3^{20}\)，显然是不行的。

考虑折半搜索，先处理一半的数据，再在另一半中找和，判断状态是否相同加和即可。

注意不要算0.

upd : 注意判重

#include <bits/stdc++.h>
using namespace std;
const int MAXN = 1 << 22;
int n;
int a[41];
int vis[MAXN];
int ans;
vector<int>v[MAXN];
map<int,int> mp;
map <int ,int > jud;
int cnt;
#define mid n / 2
void dfs1(int now,int sum,int sta) {
    if(now > mid) {
        if(!mp.count(sum)) {
            mp[sum] = ++cnt;
        }
        v[mp[sum]].push_back(sta);
        return;
    }
    dfs1(now + 1,sum + a[now],sta | (1 << (now - 1)));  
    dfs1(now + 1,sum - a[now],sta | (1 << (now - 1)));
    dfs1(now + 1,sum,sta);
}
 
void dfs2(int now,int sum,int sta) {
    if(now > n) {
        if(!mp.count(sum)) {
            return;
        }
        int res = mp[sum];
        if(jud[res] == 0) {
            sort(v[res].begin(),v[res].end());
            v[res].resize(unique(v[res].begin(),v[res].end()) - v[res].begin());
            jud[res] = 1;
        }
        for(int i = 0;i < v[res].size(); ++i) {
            if(!vis[v[res][i] | sta]) {
                ans ++;
            }
            vis[v[res][i] | sta] = 1;
        }
        return;
    }
    dfs2(now + 1,sum + a[now],sta | (1 << (now - 1)));
    dfs2(now + 1,sum - a[now],sta | (1 << (now - 1)));
    dfs2(now + 1,sum,sta);
}
 
int main () {
    cin >> n;
    for(int i = 1;i <= n; ++i) {
        cin >> a[i];
    }
    dfs1(1,0,0);
    dfs2(mid+1,0,0);
    cout<<ans - 1<<endl;
    return 0;
}

bzoj 3211

区间修改和查询

树状数组维护即可，对于开根操作暴力即可，数据比较水...

或者写个线段树维护单点开根的\(tag\)。

#include <bits/stdc++.h>
#pragma GCC optimize(2)
using namespace std;
#define lowbit(x) x & -x
#define int long long
const int MAXN = 3e6 + 10;
int c[MAXN];
int a[MAXN];
int num[MAXN];
int n;
int m;
int read () {
	int q=0,f=1;char ch=getchar();
	while(!isdigit(ch)){
		if(ch=='-')f=-1;ch=getchar();
	}
	while(isdigit(ch)) {
		q=q*10+ch-'0';ch=getchar();
	}
	return q*f;
}
void add(int x,int k) {
	for(int i = x;i <= n;i += lowbit(i)) {
		c[i] += k;
	}return;
}

int ck(int x) {
	int ans = 0;
	for(int i = x;i; i -= lowbit(i)) {
		ans += c[i];
	}
	return ans;
}

int get1(int l,int r) {
	return r - l + 1 - (num[r] - num[l - 1]);
}

void change(int l,int r) {
	if(ck(r) - ck(l - 1) <= get1(l,r)) return;
	for(int i = l;i <= r; ++i) {
		if(a[i] > 1) {
			int tmp = sqrt(a[i]);
			add(i,tmp - a[i]);
			a[i] = tmp;
		}
	}
}
//#include <ctime>

//clock_t st,ed;
signed main () {
	////freopen("11.in","r",stdin);
	//freopen("ans.out","w",stdout);
	//st = clock();
	n = read();
	for(int i = 1;i <= n; ++i) {
		a[i] = read();
		if(a[i] == 0) {
			num[i] = num[i - 1] + 1;
		}
		else num[i] = num[i - 1];
		add(i,a[i]);
	}
	m = read();
	for(int i = 1;i <= m; ++i) {
		int opt = read(),l = read(),r = read();
		if(opt == 1){
			cout<<ck(r) - ck(l - 1)<<endl;
		}
		if(opt == 2) {
			change(l,r);
		}
	}
//	ed = clock();
//	cout<<(ed - st) / 1000<<endl;
	return 0;
}