Prime Land

http://poj.org/problem?id=1365

题意:给定一个数字n的拆分形式,然后让你求解n-1的值;

解析:直接爆搞

// File Name: poj1365.cpp
// Author: bo_jwolf
// Created Time: 2013骞?0鏈?9鏃?鏄熸湡涓?21:29:25

#include<vector>
#include<list>
#include<map>
#include<set>
#include<deque>
#include<stack>
#include<bitset>
#include<algorithm>
#include<functional>
#include<numeric>
#include<utility>
#include<sstream>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<ctime>

using namespace std;
vector<int> prime, ans;
const int maxn = 40000;
int unprime[ maxn ];
int main(){
	unprime[ 0 ] = unprime[ 1 ] = true;
	for( int i = 2; i < maxn; ++i ){
		if( !unprime[ i ] ){
			prime.push_back( i );
			for( int j = i + i; j < maxn; j += i ){
				unprime[ j ] = true;
			}
		}
	}
	int n, p;
	string line;
	while( getline( cin, line ),line[ 0 ] != 48 ){
		ans.clear();
		istringstream stream( line );
			long long sum = 1;
		while( stream >> n >> p ){
			while( p-- ){
				sum *= n;
			}
		}		sum -= 1;
		for( int i = prime.size() - 1; i >= 0; --i ){
			if( sum % prime[ i ] == 0 ){
							ans.push_back( prime[ i ] );
				int temp = 0;
				while( sum % prime[ i ] == 0 ){
					sum /= prime[ i ];
					temp++;
				}
				ans.push_back( temp );
			}
		}
		for( int i = 0 ; i < ans.size(); ++i ){
			cout << ans[ i ] << ( i == ans.size() - 1 ?'\n':' ' );
		}	
	}
return 0;
}


 

% 参数设置 grid_size = 50; % 500m 10m land_size = 500; tree_area = 10; safety_radius = 2.5; heights = [5, 10, 15, 20, 25]; canopy_radius = [2.8, 5.5, 8.5, 11.9, 14.5]; % 定义最大树木数目 maximum_trees = grid_size^2; % 网格中最多能种植的树木数目 % 输入已经种植的树木数目 N_prime = input('已经种植的树木数目: '); % 初始化变量 x = zeros(grid_size, grid_size); h = ones(grid_size, grid_size) * 5; % 假设所有树的初始高度为5米 % 初始化总成本 total_cost = 0; % 遍历网格 for i = 1:grid_size for j = 1:grid_size % 检查安全距离 safe = true; for k = max(1, i-1):min(grid_size, i+1) for l = max(1, j-1):min(grid_size, j+1) if i ~= k || j ~= l if sqrt((i-k)^2 + (j-l)^2) * tree_area < 2 * safety_radius safe = false; break; end end end if ~safe break; end end % 如果满足安全距离条件,尝试种植树木 if safe && N_prime < maximum_trees x(i, j) = 1; % 计算最佳树高 min_cost = inf; best_height = 0; for height = heights canopy_r = interp1(heights, canopy_radius, height); if (i-1) * tree_area + canopy_r <= land_size && (j-1) * tree_area + canopy_r <= land_size cost = 10 * height + 10; if cost < min_cost min_cost = cost; best_height = height; end end end h(i, j) = best_height; N_prime = N_prime + 1; % 更新已种植的树木数目 end end end % 计算结果 remaining_trees = maximum_trees - N_prime; % 在已经种植的树木基础上还能种植的树木数目 total_cost = sum(sum((h * 10 + 10) .* x)); fprintf('在已经种植的树木基础上还能种植的树木数目: %d\n', remaining_trees); fprintf('总成本: %d\n', total_cost);请分析一下此代码的错误
05-25
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值