UVALive 6838 Flipping Parentheses(线段树、单点更新、区间查询)

题目链接:
UVALive 6838 Flipping Parentheses
题意:
给出一个长度为n个串,每个字符只能是’(‘或’)’,而且左括号和右括号个数相等,在操作的过程中,要保证这个串的任意前缀串的左括号个数都要大于等于右括号的个数。
有Q个操作。每次操作输入一个坐标t将(坐标从1~n)这个坐标下的括号取反,即左括号变右括号,右括号变左括号。
对每次输入输出最靠前的一个坐标,使得改变这个坐标下的字符(取反),仍然满足任意前缀子串的左括号个数大于等于右括号个数?
分析:
将左括号个数减去右括号个数设为平衡度。显然整个串的平衡度为0,而且任意前缀串的平衡度大于等于0.
假设将id位置的括号从左变成右,那么从这个位置一直到串尾的每个位置上的前缀串的平衡度都增加2,相反,如果由右括号变成左括号,那么从这个位置一直到串尾的每个位置上的前缀串的平衡度都减少2.
首先考虑将第t位置的括号由左变右的情况,。显然需要找到最靠前的右括号将其变为左括号,假设坐标为id,如何找到id呢?可以在线段树中保存当前区间的右括号的个数,那么就相当于单点查询了。找到id后还需注意,因为我们需要的是操作完成后最靠前的右括号,如果id>t,那么因为需要将t位置的左括号改为右括号,那么这是id应该就是t,先变成右括号,再变成左括号。然后就是区间更新,更新从id位置到串尾的线段树区间最小前缀平衡度。
为什么要保存区间最小前缀平衡度呢?是因为接下来将第t位置的括号由右变左需要用到。我们需要找到的是一个左括号,同样假设下标id,将其变为右括号,而且保持任意前缀
子串的平衡度非负,这是我们应该先更新在查找。如果从前往后找的话,那么从id以后的子串的前缀平衡度就难以保证了。所以我们需要从后往前找。
找到一个这样的最靠前的位置:从这个位置往后(包括这个位置)的前缀平衡度恒大于等于2,因为左括号变右相当于平衡度减2.也就是从后往前找到第一个前缀平衡度小于2的后一个位置就是我们要随着t改变的位置。

更新区间前缀平衡度最小值要用lazy标记。
还有一个特例需要注意。就是要改变的位置是第一个括号和最后一个括号时,只能是在把它变回来。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <climits>
#include <cmath>
#define lson(x) (x<<1)
#define rson(x) ((x<<1)|1)
using namespace std;
const int MAX_N=300010;

int n,Q;
char s[MAX_N];
int sum[MAX_N];

struct Segtree{
    int left,right;
    int num1,num2,balance_min,add;
}segtree[MAX_N<<2];

inline void push_up(int cur)
{
    segtree[cur].num1=segtree[lson(cur)].num1+segtree[rson(cur)].num1;
    segtree[cur].num2=segtree[lson(cur)].num2+segtree[rson(cur)].num2;
    segtree[cur].balance_min=min(segtree[lson(cur)].balance_min,segtree[rson(cur)].balance_min);
}

inline void push_down(int cur)
{
    int add=segtree[cur].add;
    if(add!=0){
        segtree[lson(cur)].balance_min+=add;
        segtree[rson(cur)].balance_min+=add;

        segtree[lson(cur)].add+=add;
        segtree[rson(cur)].add+=add;

        segtree[cur].add=0;
        return ;
    }
}

inline void build(int left,int right,int cur)
{
    segtree[cur].left=left;
    segtree[cur].right=right;
    segtree[cur].add=0;
    if(left==right){
        segtree[cur].balance_min=sum[left];
        int val=s[left]=='('?1:0;
        segtree[cur].num1=val;
        segtree[cur].num2=1-val;
        return ;
    }
    int mid=(left+right)>>1;
    build(left,mid,lson(cur));
    build(mid+1,right,rson(cur));
    push_up(cur);
}

inline int query1(int a,int b,int cur)
{//查询第一个右括号
    int left=segtree[cur].left;
    int right=segtree[cur].right;
    if(left==right){
        //printf("query1=%d\n",left);
        return left;
    }
    push_down(cur);
    int tmp=segtree[lson(cur)].num2;
    int mid=(left+right)>>1;
    int res;
    if(tmp) res=query1(left,mid,lson(cur)); //左儿子有右括号
    else res=query1(mid+1,right,rson(cur));
    push_up(cur);
    return res;
}

inline void update_dandian(int cur,int id,int flag)
{//flag=0单点左改右,flag=1单点右改左
    int left=segtree[cur].left;
    int right=segtree[cur].right;
    if(left==right){
        //printf("update1_dandian:left==%d\n",left);
        segtree[cur].num1=flag;
        segtree[cur].num2=1-flag;
        return;
    }
    push_down(cur);
    int mid=(left+right)>>1;
    if(id<=mid) update_dandian(lson(cur),id,flag);
    else update_dandian(rson(cur),id,flag);
    push_up(cur);
}

inline void update_qujian(int a,int b,int cur,int flag)
{
    int left=segtree[cur].left;
    int right=segtree[cur].right;
    if(left==a&&right==b){
        segtree[cur].add+=flag;
        segtree[cur].balance_min+=flag;
        return ;
    }
    push_down(cur);
    int mid=(left+right)>>1;
    if(a>mid) update_qujian(a,b,rson(cur),flag);
    else if(b<=mid) update_qujian(a,b,lson(cur),flag);
    else {
        update_qujian(a,mid,lson(cur),flag);
        update_qujian(mid+1,b,rson(cur),flag);
    }
    push_up(cur);
}

inline int query2(int cur)
{
    int left=segtree[cur].left;
    int right=segtree[cur].right;
    if(left==right){
        return left;
    }
    push_down(cur);
    int mid=(left+right)>>1;
    int res;
    int tmp1=segtree[lson(cur)].balance_min;
    int tmp2=segtree[rson(cur)].balance_min;
    //printf("tmp1=%d tmp2=%d left=%d right=%d\n",tmp1,tmp2,left,right);
    if(tmp2<2) res=query2(rson(cur)); //从右往前找到第一个前缀平衡度小于2的位置
    else res=query2(lson(cur));
    push_up(cur);
    return res;
}

int main()
{
    freopen("Gin.txt","r",stdin);
    while(~scanf("%d%d",&n,&Q)){
        scanf("%s",s+1);
        memset(sum,0,sizeof(sum));
        for(int i=1;i<=n;i++){
            sum[i]=sum[i-1]+(s[i]=='('?1:(-1));
        }
        build(1,n,1);
        for(int i=0;i<Q;i++){
            int t;
            scanf("%d",&t);
            if(s[t]=='('){ //左改右
                if(t==1){
                    printf("1\n");
                    continue;
                }
                int id=query1(1,n,1); //id是现有的第一个右括号位置
                if(id>=t){ 
                    printf("%d\n",t);
                    continue;
                }
                s[id]='(';
                update_dandian(1,id,1);//单点更新,右改左
                update_qujian(id,n,1,2);//区间更新,右改左
                s[t]=')';
                update_dandian(1,t,0);//单点更新,左改右
                update_qujian(t,n,1,-2);//区间更新,左改右
                printf("%d\n",id);
            }else { //右改左
                if(t==n){
                    printf("%d\n",n);
                    continue;
                }
                update_dandian(1,t,1);//单点更新,右改左
                update_qujian(t,n,1,2);//区间更新,右改左
                s[t]='(';
                int id=query2(1)+1;
                update_dandian(1,id,0);//单点更新,左改右
                update_qujian(id,n,1,-2);//区间更新,左改右
                if(id>t) id=t; 
                s[id]=')';
                printf("%d\n",id);
            }
        }

    }
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是LDPC Modified Weighted Bit Flipping算法的C语言实现: ```c #include <stdio.h> #include <stdlib.h> #include <time.h> #include <math.h> // 定义LDPC码参数 #define N 256 // 码字长度 #define K 128 // 信息长度 #define M 128 // 码字检验矩阵的行数 #define Q 6 // 权重因子,取值范围[1, 8] // 定义模拟参数 #define MAX_ITER 100 // 最大迭代次数 #define ERR_THRESHOLD 1e-6 // 误码率阈值 // 定义随机数生成器 #define RAND() ((double)rand() / RAND_MAX) // 生成LDPC码矩阵,采用随机方式生成 void generateLdpcMatrix(int H[M][N]) { int i, j, k, rowWeight, colWeight; // 初始化H矩阵 for (i = 0; i < M; i++) { for (j = 0; j < N; j++) { H[i][j] = 0; } } // 随机生成每行的权重 for (i = 0; i < M; i++) { rowWeight = rand() % (Q + 1) + 1; // 随机生成每行的列下标 for (j = 0; j < rowWeight; j++) { do { k = rand() % N; } while (H[i][k] != 0); H[i][k] = 1; } } // 计算每列的权重 for (j = 0; j < N; j++) { colWeight = 0; for (i = 0; i < M; i++) { colWeight += H[i][j]; } // 如果权重为0,则重新生成一行 if (colWeight == 0) { do { i = rand() % M; } while (H[i][j] != 0); H[i][j] = 1; } } } // 生成随机的信息位 void generateMessage(int message[K]) { int i; for (i = 0; i < K; i++) { message[i] = (int)(RAND() + 0.5); } } // 生成LDPC码字 void generateCodeword(int H[M][N], int message[K], int codeword[N]) { int i, j; // 初始化码字向量 for (i = 0; i < N; i++) { codeword[i] = 0; } // 计算码字向量 for (i = 0; i < M; i++) { for (j = 0; j < N; j++) { codeword[j] ^= H[i][j] & message[i]; } } } // 生成高斯噪声 void generateNoise(double sigma, double noise[N]) { int i; for (i = 0; i < N; i++) { noise[i] = sigma * sqrt(-2.0 * log(RAND())) * sin(2.0 * M_PI * RAND()); } } // AWGN信道 void awgnChannel(double sigma, int codeword[N], double received[N]) { int i; generateNoise(sigma, received); for (i = 0; i < N; i++) { received[i] += codeword[i]; } } // 计算似然比 void computeLLR(double sigma, double received[N], double llr[N]) { int i; for (i = 0; i < N; i++) { llr[i] = 2.0 * received[i] / (sigma * sigma); } } // 计算权重向量 void computeWeight(int H[M][N], double llr[N], double weight[N]) { int i, j; for (j = 0; j < N; j++) { weight[j] = 0.0; for (i = 0; i < M; i++) { weight[j] += H[i][j] * llr[j]; } } } // 计算误码率 double computeBer(int message[K], int decMessage[K]) { int i; double err = 0.0; for (i = 0; i < K; i++) { if (message[i] != decMessage[i]) { err += 1.0; } } return err / K; } // Modified Weighted Bit Flipping译码算法 void MWBFDecoder(int H[M][N], double llr[N], int decMessage[K], int *iterCnt) { int i, j, k, l, rowWeight, colWeight, maxIdx; double weight[N], prob[N], maxProb, tmp; int syndrome[M], decCodeword[N], dec[N]; // 初始化译码结果 for (i = 0; i < K; i++) { decMessage[i] = 0; } // 初始化译码次数 *iterCnt = 0; // 初始化译码向量 for (i = 0; i < N; i++) { decCodeword[i] = (llr[i] < 0.0) ? 1 : 0; } // 迭代译码 while (1) { // 判断是否达到最大迭代次数 if (*iterCnt >= MAX_ITER) { break; } // 计算校验和 for (i = 0; i < M; i++) { syndrome[i] = 0; for (j = 0; j < N; j++) { syndrome[i] ^= H[i][j] & decCodeword[j]; } } // 如果校验和为0,则译码成功 for (i = 0; i < M; i++) { if (syndrome[i] != 0) { break; } } if (i >= M) { break; } // 计算权重向量 computeWeight(H, llr, weight); // 根据权重向量进行修正 for (i = 0; i < N; i++) { if (decCodeword[i] == 0) { prob[i] = exp(-weight[i]); } else { prob[i] = exp(weight[i]); } } for (i = 0; i < K; i++) { dec[i] = 0; rowWeight = 0; for (j = 0; j < N; j++) { if (H[i][j] != 0) { rowWeight++; } } for (j = 0; j < rowWeight; j++) { maxProb = -1.0; for (k = 0; k < N; k++) { if (H[i][k] != 0 && prob[k] > maxProb) { maxProb = prob[k]; maxIdx = k; } } dec[i] ^= decCodeword[maxIdx]; prob[maxIdx] = -1.0; } decMessage[i] = dec[i]; } // 计算修正后的码字向量 for (i = 0; i < N; i++) { decCodeword[i] = 0; colWeight = 0; for (j = 0; j < M; j++) { if (H[j][i] != 0) { colWeight++; tmp = 2.0 * weight[i] / rowWeight; decCodeword[i] ^= (tmp >= 0.0) ? 0 : 1; tmp = fabs(tmp); if (tmp > 1.0) { tmp = 1.0; } prob[i] = tmp; } } } // 更新译码次数 (*iterCnt)++; } } int main() { int i, j, iterCnt, message[K], decMessage[K], codeword[N], decCodeword[N]; double sigma, received[N], llr[N], weight[N], ber; // 初始化随机数种子 srand((unsigned)time(NULL)); // 生成LDPC码矩阵 int H[M][N]; generateLdpcMatrix(H); // 生成随机信息位 generateMessage(message); // 生成LDPC码字 generateCodeword(H, message, codeword); // 设置信噪比 sigma = 0.2; // AWGN信道传输 awgnChannel(sigma, codeword, received); // 计算似然比 computeLLR(sigma, received, llr); // Modified Weighted Bit Flipping译码 MWBFDecoder(H, llr, decMessage, &iterCnt); // 误码率 ber = computeBer(message, decMessage); printf("BER = %e (iterCnt=%d)\n", ber, iterCnt); return 0; } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值