Flipping Game

PTA | 程序设计类实验辅助教学平台 (pintia.cn)icon-default.png?t=N7T8https://pintia.cn/problem-sets/91827364500/exam/problems/91827370509?type=7&page=0题意 

        给定两个长度为n的字符串,代表初始灯泡状态和目标灯泡状态,经过k个回合,按下m个不同的开关,问有多少种方案能使灯泡从初始状态到目标状态

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <string>
#include <cstring>
#include <vector>
#include <queue>
#include <deque>
#include <cstdlib>
#include <stack>
#include <list>
#include <set>
#include <map>
#include <unordered_map>

using namespace std;

const int mod = 998244353;

int t;
int n, k, m, ans;
char start[105], target[105];
long long C[105][105], dp[105][105];

int main()
{
    C[0][0] = C[1][0] = C[1][1] = 1;//初始化组合数数组
    for(int i = 2;i < 105;i ++ )
    {
        C[i][0] = C[i][i] = 1;

        for(int j = 1;j < i;j ++ )
        {
            C[i][j] = (C[i - 1][j] + C[i - 1][j - 1]) % mod;
        }

    }

    scanf("%d", &t);

    while(t -- )
    {
        ans = 0;

        scanf("%d %d %d", &n, &k, &m);
        scanf("%s %s", start, target);

        for(int i = 0;i < k;i ++ )
        {
            for(int j = 0;j <= 105;j ++ )
            {
                dp[i][j] = 0;
            }
        }

        int num = 0;//记录初始灯状态与目标灯状态相差个数

        for(int i = 0;i < n;i ++ )
        {
            if(start[i] != target[i])
            {
                num ++;
            }
        }

        //dp[i][j]代表在进行第i次操作后,仍有j个灯泡状态与目标状态不符
        dp[0][num] = 1;

        for(int i = 0;i < k;i ++ )
        {
            for(int j = 0;j <= n;j ++ )
            {
                if(dp[i][j])
                {
                    for(int k = 0;k <= j;k ++ )
                    {
                        if(n - j >= m - k)
                        {
                            int to = j - k + m - k;
                            //C[j][k]表示在j个不符灯泡中挑选k个灯泡的情况数 C[n - j][m - k]表示在符合的灯泡中挑选m - k的情况数
                            dp[i + 1][to] = (dp[i + 1][to] + C[j][k] * C[n - j][m - k] % mod * dp[i][j] % mod) % mod;
                        }
                    }
                }
            }
        }
        printf("%lld\n", dp[k][0]);
    }

    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
以下是一个基于C语言的LDPC Modified Weighted Bit Flipping算法的代码示例: ```c #include<stdio.h> #include<stdlib.h> #include<math.h> #include<time.h> #define N 2048 //码长 #define K 1024 //信息位数 #define M 1024 //校验位数 #define p 0.01 //BSC信道误码率 #define MAX_ITER 50 //最大迭代次数 #define tau 2.0 //重要性参数 double H[M][N]; //LDPC矩阵 int d[N]; //接收端接收的码字 double r[N]; //接收端的对数似然比 double q[N]; //对数似然比的绝对值 double x[N]; //解码后的信息位 int perm[N]; //随机置换序列 double c[M]; //校验和 int iter_cnt = 0; //迭代次数计数器 //生成LDPC矩阵 void gen_LDPC_matrix() { int i,j,k; int col[N]; int row[M]; srand((unsigned)time(NULL)); for(i=0;i<M;i++) for(j=0;j<N;j++) H[i][j]=0; for(i=0;i<N;i++) perm[i]=i; for(i=0;i<M;i++) row[i]=i; for(i=0;i<N;i++) col[i]=i; for(i=0;i<N;i++) { k = rand()%(N-i) + i; j = perm[i]; perm[i] = perm[k]; perm[k] = j; } for(i=0;i<K;i++) { for(j=0;j<N/K;j++) { H[i][col[i*N/K+j]] = 1; } } for(i=K;i<M;i++) { for(j=0;j<N/K;j++) { H[i][col[(rand()%(N/K))+j*N/K]] = 1; } } } //BSC信道模拟 void BSC_channel() { int i; for(i=0;i<N;i++) { if((double)rand()/RAND_MAX < p) d[i] = 1-d[i]; r[i] = log((1-p)/p); if(d[i]==1) r[i] = -r[i]; q[i] = fabs(r[i]); } } //判断校验方程是否满足 int check_equations() { int i,j; int flag = 0; for(i=0;i<M;i++) { c[i] = 0; for(j=0;j<N;j++) { if(H[i][j]==1) c[i] += x[j]; } c[i] = fmod(c[i],2); if(c[i]!=0) flag = 1; } return flag; } //Modified Weighted Bit Flipping算法 void MWBF() { int i,j,k; double w[N]; double sum_w = 0; double f[N]; double delta[N]; double sum_delta = 0; double p0,p1; while(iter_cnt < MAX_ITER && check_equations()) { iter_cnt++; for(i=0;i<N;i++) { w[i] = 1/(1+exp(q[i])); if(x[i]==1) w[i] = 1-w[i]; sum_w += w[i]; } for(i=0;i<N;i++) { f[i] = 1; for(j=0;j<M;j++) { if(H[j][i]==1) f[i] *= (2*c[j]-1); } delta[i] = log((1-w[i])/w[i]) + f[i]/sum_w/tau; sum_delta += fabs(delta[i]); } while(sum_delta>1e-6) { k = 0; for(i=1;i<N;i++) { if(delta[i]>delta[k]) k = i; } x[k] = 1-x[k]; p0 = p/(1-p); p1 = (1-p)/p; for(i=0;i<M;i++) { if(H[i][k]==1) { if(c[i]==0) p1 *= (1-2*w[k]); else p0 *= (1-2*w[k]); } } r[k] = log(p0/p1); if(x[k]==1) r[k] = -r[k]; q[k] = fabs(r[k]); sum_delta = 0; for(i=0;i<N;i++) { f[i] = 1; for(j=0;j<M;j++) { if(H[j][i]==1) f[i] *= (2*c[j]-1); } delta[i] = log((1-w[i])/w[i]) + f[i]/sum_w/tau; sum_delta += fabs(delta[i]); } } sum_w = 0; sum_delta = 0; } } int main() { int i; double BER; gen_LDPC_matrix(); for(i=0;i<N;i++) d[i] = rand()%2; BSC_channel(); MWBF(); BER = 0; for(i=0;i<K;i++) BER += x[i] ^ d[i]; BER /= K; printf("BER = %f\n",BER); return 0; } ``` 该代码实现了一个N=2048,K=1024的LDPC码的解码,使用的是BSC信道。算法采用了Modified Weighted Bit Flipping算法,重要性参数tau默认为2。代码中使用了对数似然比,避免了浮点数下溢问题。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Miracleresss

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值