POJ2083 Fractal【分形回溯】

本文介绍了一种递归生成盒分形的算法,通过定义递归函数printBox(n,x,y),在指定坐标绘制n度盒分形。文章详细解释了递归过程,包括递归边界条件和递归调用的坐标计算。

n度的盒分形的规模为3(n-1),即n度的盒分形图为一个长宽为3(n-1)的正方形。
设置递归函数printBox(n,x,y)生成以坐标(x,y)为左上角的n度盒分形。
1)递归边界: 若n=1,则在(x,y)输出‘X’
2)若n>1,则计算n-1度的盒子的规模 m = 3^(n-2),分别在左上方, 右上方,中间,左下方和右下方画出5个n-1度的盒子。
对于左上方的n-1度的盒子,左上角的坐标为(x,y),递归printBox(n-1,x,y)生成;
对于右上方的n-1度的盒子,左上角的坐标为(x+2m,y),递归printBox(n-1,x+2m,y)生成;
对于中间的n-1度的盒子,左上角的坐标为(x+m,y+m),递归printBox(n-1,x+m,y+m)生成;
对于左下方的n-1度的盒子,左上角的坐标为(x,y+2m),递归printBox(n-1,x,y+2m)生成;
对于右上方n-1度的盒子,左上角的坐标为(x+2m,y+2m),递归printBox(n-1,x+2m,y+2m)生成

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<math.h>
using namespace std;

char a[735][735];

void dfs(int cur,int x,int y)
{
    if(cur==1)
    {
        a[x][y]='X';
        return ;
    }
    int s=pow(3,cur-2);
    dfs(cur-1,x,y);
    dfs(cur-1,x,y+2*s);
    dfs(cur-1,x+s,y+s);
    dfs(cur-1,x+2*s,y);
    dfs(cur-1,x+2*s,y+2*s);
}

int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        if(n==-1)
            break;
        memset(a,' ',sizeof(a));
        dfs(n,1,1);
        int s=pow(3,n-1);
        for(int i=1;i<=s;i++)
        {
            a[i][s+1]='\0';
        }
        for(int i=1;i<=s;i++)
        {
            printf("%s\n",a[i]+1);
        }
        printf("-\n");
    }
    return 0;
}

POJ 2083 是一个经典的编程竞赛题目,名为 "Get Many Persimmons"。该问题的核心在于通过模拟和贪心策略找到最优解法。以下是关于 POJ 2083 的解决方案及相关信息。 --- ### 题目概述 在一个二维平面上有若干个柿子树节点,每个节点有一个坐标 (x, y) 和对应的权值 w_i。你可以从原点出发,每次可以选择向右或者向上移动一格,并收集经过的柿子树上的果实。目标是在不超过最大容量 K 的情况下最大化收获的柿子数量。 --- ### 解决方案分析 #### 方法选择 这是一个典型的 **动态规划** 或者 **背包问题变体** 的应用案例[^1]。具体而言,我们需要考虑以下几点: - 路径的选择受到方向限制(仅能向右或向上)。 - 收集柿子的过程类似于多重背包问题中的物品选取逻辑。 #### 动态规划建模 设 `dp[x][y][k]` 表示当前位于 `(x, y)` 并已携带重量为 k 的柿子的最大价值,则状态转移方程如下: \[ dp[x][y][k] = \max(dp[x-1][y][k], dp[x][y-1][k]) \] 如果当前位置 `(x, y)` 存在柿子树且其权重为 `w`、价值为 `v`,则进一步更新: \[ dp[x][y][k+w] = \max(dp[x][y][k+w], dp[x][y][k]+v) \] 边界条件设定为当超出范围时返回负无穷大以排除非法路径。 #### 时间复杂度优化 由于朴素实现的时间复杂度过高(O(N*M*K)),实际操作中可采用滚动数组技术降低空间消耗至 O(K)[^2]。 ```cpp #include <bits/stdc++.h> using namespace std; const int MAXK = 1e4 + 5; int N, M, K; pair<int, pair<int,int>> tree[105]; // weight, value, position(x,y) bool cmp(const pair<int,pair<int,int>>& a,const pair<int,pair<int,int>>& b){ if(a.second.first != b.second.first) return a.second.first<b.second.first; return a.second.second<b.second.second; } int main(){ ios::sync_with_stdio(false); cin.tie(0); while(cin>>N>>M>>K && !(N==0&&M==0)){ memset(tree,0,sizeof(tree)); for(int i=1;i<=N;i++)cin>>tree[i].second.first>>tree[i].first>>tree[i].second.second; sort(tree+1,tree+N+1,cmp); vector<vector<int>> dp(M+1,vector<int>(K+1,-1)); //-1 means not reachable yet. dp[0][0]=0; // start at origin with no persimmon collected. for(int idx=1;idx<=N;idx++){ auto &[wi,val,pos]=tree[idx]; int x=pos.first, y=pos.second; for(int prev_y=M;y-prev_y>=0;prev_y--){ for(int k=K;k-wi>=0;k--){ if(dp[prev_y][k]!=-1){ dp[y][min(k+wi,K)] = max(dp[y][min(k+wi,K)], dp[prev_y][k]+val); } } } } int res=0; for(auto &row : dp) for(auto val: row)res=max(res,val); cout<<res<<"\n"; } } ``` --- ### 关键点解析 1. **路径约束**: 只允许向右或向上移动意味着我们只需关注前驱状态即可完成状态转移[^3]。 2. **容量管理**: 类似于经典背包问题,需确保不超载的同时追求收益最大化。 3. **效率提升**: 使用滚动数组减少内存占用并提高运行速度。 ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值