Bootstrap

bitset优化

bitset优化背包
d p i , k ∣ = d p i − 1. k − j ∗ j dp_{i,k}|=dp_{i-1.k-j*j} dpi,k=dpi1.kjj 其中 i i i代表到第 i i i个区间,选完数得到的平方和为 k k k 这种情况是否存在
j j j表示第 i i i次选数为 j j j

枚举 i , j , k i,j,k i,j,k时间复杂度为 O ( n l 3 ) O(nl^3) O(nl3)

用bitset表示dp中的第二维,将枚举 k k k这一部分 O ( l 2 ) O(l^2) O(l2)优化为 O ( l 2 / 64 ) O(l^2/64) O(l2/64)

std::bitset<1000005> dp[105];
void solve()
{
	int n;
	cin>>n;
	std::vector<int> l(n+1),r(n+1);
	for (int i=1;i<=n;i++)
		std::cin>>l[i]>>r[i];
	dp[0][0]=1;
	for (int i=1;i<=n;i++)
		for (int j=l[i];j<=r[i];j++)
			dp[i]|=dp[i-1]<<(j*j);
	std::cout<<dp[n].count();
}

简单字符串问题
暴力枚举中心点 O ( n 2 ) O(n^2) O(n2),用两个bitset顺逆标记每个字母所在的位置上,之后直接枚举中心点,利用位运算求重合部分数量 O ( n 2 / 64 ) O(n^2/64) O(n2/64)

char a[100005];
void solve(){
    int n;
    cin>>n;
    vector<bitset<100000>> vis1(26),vis2(26);
    for (int i=0;i<n;i++){
        cin>>a[i];
        vis1[a[i]-'a'][i]=1;
        vis2[a[i]-'a'][n-i-1]=1;
    }
    ll ans=0;
    for (int i=0;i<n;i++){
        ans+=((vis2[a[i]-'a']>>(n-i))&(vis1[a[i]-'a']>>(i+1))).count();
    }
    cout<<ans<<"\n";
}

CF333E Summer Earnings
暴力做法直接枚举三个点,看最小的一条边,不断更新最大值 O ( n 3 ) O(n^3) O(n3)

换种角度:考虑枚举某个三角形中最小的边(两点),判断是否存在另一点到两点的距离大于此边

转换一下,先求出所有的边,然后进行从小到大排序,不断加边。如果当前加的这条边(枚举的最小的边)的两端点 x , y x,y x,y,存在合法的另外一点 z z z(此前存在 x x x z z z, y y y z z z的两条边),则该边即为答案

如果用数组打标记,枚举 z z z的方法,复杂度仍为 O ( n 3 ) O(n^3) O(n3)
所以用bitset代替数组,直接用位运算和count判断是否存在某个 z z z 复杂度降为 O ( n 3 / 64 ) O(n^3/64) O(n3/64)

struct segment{
    int x,y;
    double len;
};
double road(int x1,int y1,int x2,int y2){
    return (x1-x2)*(x1-x2)+(y1-y2)*(y1-y2);
}
int main(){
    int n;
    cin>>n;
    vector<int> ax(n+1),ay(n+1);
    for (int i=1;i<=n;i++){
        cin>>ax[i]>>ay[i];
    }

    vector<segment> b;
    for (int i=1;i<=n;i++){
        for (int j=i+1;j<=n;j++){
            b.push_back({i,j,road(ax[i],ay[i],ax[j],ay[j])});
        }
    }

    sort(b.begin(),b.end(),[&](segment i,segment j){
        return i.len>j.len;
    });

    vector<bitset<3005>> vis(n+1);
    for (auto i:b){
        if ((vis[i.x]&vis[i.y]).count()){
            cout<<fixed<<setprecision(9)<<(double)(sqrt(i.len)/2);
            return 0;
        }
        vis[i.x][i.y]=1;
        vis[i.y][i.x]=1;
    }
	return 0;
}
;