ed div2 D. Bicolored RBS*

给出n个匹配好的括号  显然n为偶数

定义:()深度1

(())深度2  以此类推

 

要求 对括号进行01染色   (当然染完色匹配是正确的)

使得  最大深度最小化

 

贪心

用链式结构   以相邻的()进行轮流01染色直到染完即可

#include<bits/stdc++.h>
using namespace std;
//input by bxd
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define repp(i,a,b) for(int i=(a);i>=(b);--i)
#define RI(n) scanf("%d",&(n))
#define RII(n,m) scanf("%d%d",&n,&m)
#define RIII(n,m,k) scanf("%d%d%d",&n,&m,&k)
#define RS(s) scanf("%s",s);
#define ll long long
#define pb push_back
#define CLR(A,v)  memset(A,v,sizeof A)
//
#define inf 0x3f3f3f3f
const int N=(int)5e5+5;
char ans[N],s[N];
int le[N],ri[N];

void del(int x)
{
    int L=le[x],R=ri[x];
    le[R]=L;
    ri[L]=R;
}

int main()
{
    int n;RI(n);RS(s+1);
    
    rep(i,1,n)
    {
        le[i]=i-1;
        ri[i]=i+1;
    }
    ri[0]=1;
    le[n+1]=n;
    
    int cnt=0;
    char t='0';
    while(cnt<n)
    {
        int head=ri[0];
        while(head<=n)
        {
            
            if(s[head]=='('&&s[ri[head]]==')' )
            {   
                int q=head;int w=ri[head];
                ans[head]=ans[ri[head]]=t;cnt+=2;
                head=ri[ri[head]];
                del(q);del(w);
            }
            else head=ri[head];
        }
        if(t=='0')t='1';
        else t='0';
    }

    rep(i,1,n)
    cout<<ans[i];

    return 0;
}
View Code

 

今早重跑居然T了...

直接递归比较快  

on

#include<bits/stdc++.h>
using namespace std;
//input by bxd
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define repp(i,a,b) for(int i=(a);i>=(b);--i)
#define RI(n) scanf("%d",&(n))
#define RII(n,m) scanf("%d%d",&n,&m)
#define RIII(n,m,k) scanf("%d%d%d",&n,&m,&k)
#define RS(s) scanf("%s",s);
#define ll long long
#define pb push_back
#define CLR(A,v)  memset(A,v,sizeof A)
//
#define inf 0x3f3f3f3f
const int N=2e5+5;
int ans[N],n,t[N],sta[N],pos;
char s[N];

void color(int L,int R,int k)
{
    rep(i,L,R)
        ans[i]=ans[t[i]]=k,color(i+1,t[i]-1,k^1),i=t[i];
}

int main()
{
    pos=0;
    RI(n);RS(s+1);
    rep(i,1,n)
        if(s[i]=='(')sta[++pos]=i;
        else
            t[i]=sta[pos],t[sta[pos]]=i,pos--;

    color(1,n,0);
    rep(i,1,n)
        printf("%d",ans[i]);

}
View Code

 

转载于:https://www.cnblogs.com/bxd123/p/10873183.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
for a2=1:zx ZKX=a2*gj; for b2=1:zy ZKY=b2*gj; r=sqrt((ZKX-JSX)^2+(ZKY-JSY)^2); Xt=abs(JSX-ZKX); %计算点到钻孔的x距离 if(a2==2&&b2==1) continue; end rbs=((ZKX-a1*gj)^2+(ZKY-b1*gj)^2)/gj^2+1; for j=1:nj if (j==1) z=1; elseif(j==nj) z=H-1; else z=(j-1)*dz; end for k=1:nj if(k==1) a=0; b=dz/2; elseif(k==nj) a=H-dz/2; b=H; else a=(2*k-3)*dz*0.5; b=(2*k-1)*dz*0.5; end rydis=(a+b)/2; jsdis=z; [v,Rap,Iap,Rlamd,Ilamd] = untitled55(rydis,jsdis); rr=r; parfor i=1:LL t=i*dt; aa=integral(@(x)0.25*exp(v*Xt*0.5 / Rap)*exp(-v * sqrt(rr*rr + (z - x).*(z - x))*0.5 / Rap).*erfc((sqrt(rr*rr + (z - x).*(z - x)) - v * t)*0.5 / sqrt(Rap*t))./sqrt(rr*rr + (z - x).*(z - x))/(2 * 3.1415926*Rlamd),a,b); ab=integral(@(x)0.25*exp(v*Xt*0.5 / Rap)*exp(v*sqrt(rr*rr + (z - x).*(z - x))*0.5 / Rap).*erfc((sqrt(rr*rr + (z - x).*(z - x)) + v * t)*0.5 / sqrt(Rap*t))./sqrt(rr*rr + (z - x).*(z - x))/(2 * 3.1415926*Rlamd),a,b); ac=integral(@(x)0.25*exp(v*Xt*0.5 / Iap)*exp(-v * sqrt(rr*rr + (z + x).*(z + x))*0.5 / Iap).*erfc((sqrt(rr*rr + (z + x).*(z + x)) - v * t)*0.5 / sqrt(Iap*t))./sqrt(rr*rr + (z+ x).*(z + x))/(2 * 3.1415926*Ilamd),a,b); ad=integral(@(x)0.25*exp(v*Xt*0.5 / Iap)*exp(v*sqrt(rr*rr + (z + x).*(z + x))*0.5 / Iap).*erfc((sqrt(rr*rr + (z + x).*(z + x)) + v * t)*0.5 / sqrt(Iap*t))./sqrt(rr*rr + (z + x).*(z + x))/(2 * 3.1415926*Ilamd),a,b); aa(isnan(aa)) = 0;ab(isnan(ab)) = 0;ac(isnan(ac)) = 0; ad(isnan(ad)) = 0; Tj(i,j,k,rbs)=(aa+ab-ac-ad); %Tj(i,j,k,rbs)=(aa+ab); end end end end end优化 代码
05-29
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值