螺旋方阵

本文探讨如何通过给定的奇数n和整数m,在由1到n*n构成的螺旋矩阵中找到对应的行号和列号。算法通过计算对角线数值、确定所在层以及计算与目标值的距离来定位目标值的位置。

Problem Description

Given an odd number n, we can arrange integers from 1 to n*n in the shape of a spiral. The figure below illustrates the spiral made by integers from 1 to 25.

1 2 3 4 5

------------------------

1 | 21 22 23 24 25

2 | 20 7 8 9 10

3 | 19 6 1 2 11

4 | 18 5 4 3 12

5 | 17 16 15 14 13

As we see above, each position in the spiral corresponds to a unique integer. For example, the number in row 1, column 1 is 21, and integer 16 is in row 5, column 2.
Now, given the odd number n (1<=n<=32768), and an integer m (1<=m<=n*n), you should write a program to find out the position of m.

Input

The first line of the input is a positive integer T(T<=20). T is the number of the test cases followed. Each case consists of two integer n and m as described above.

Output

For each case, output the row number and column number that the given integer is in, separated by a single whitespace. Please note that the row and column number are both starting from 1.

Sample Input

3
3 9
5 21
5 16

Sample Output

1 3
1 1
5 2

Author

HYNU 


代码:

/*题目数值大,要推理出来,找下标规律
先把对角线一半的数值推出来,并保存
再确定在哪一层之后,算出距离,
*/
#include<cstdio>
#include<iostream>
using namespace std;
int a[32769];
int main()
{
    int i, j, t, m, n;
    int count = 0, x, y;
    while(scanf("%d", &t) != EOF)
    {
        while(t--)
        {

            scanf("%d%d", &n, &m);
            count = 0;
            for(i = 0; i < n; i++) //得出对角线的数值
                if(i == 0) a[i] = 1;
                else a[i] = 8 * i - 2 + a[i - 1];

            while(a[count] < m)  count++; //确定在哪一层

            i = j = n / 2 + 1 - count; //确定行列坐标的大概位置

            x = a[count] - m; //距离 左下半部分(竖行) 的距离(从行列到具体数值的距离)
            if(x <= count * 4)
            {
                if(x <= count * 2) //列
                    i = i + x;
                else//行
                {
                    i = i + count * 2;
                    j = j + x - count * 2;
                }
            }
            else
            {
                y = x - count * 4; //距离 右上半部分(竖行) 的距离(从行列到具体数值的距离)
                if(y <= count * 2 - 1) //列
                {
                    i = i + count * 2 - y;
                    j = j + count * 2;
                }
                else//行(比本行的对角数还要大)
                {
                    y = y - (count * 2 - 1);
                    i = i + 1;
                    j = j + count * 2 - y;
                }
            }
            printf("%d %d\n", i, j);
        }
    }
    return 0;
}


内容概要:本文介绍了一种基于蒙特卡洛模拟和拉格朗日优化方法的电动汽车充电站有序充电调度策略,重点针对分时电价机制下的分散式优化问题。通过Matlab代码实现,构建了考虑用户充电需求、电网负荷平衡及电价波动的数学模【电动汽车充电站有序充电调度的分散式优化】基于蒙特卡诺和拉格朗日的电动汽车优化调度(分时电价调度)(Matlab代码实现)型,采用拉格朗日乘子法处理约束条件,结合蒙特卡洛方法模拟大量电动汽车的随机充电行为,实现对充电功率和时间的优化分配,旨在降低用户充电成本、平抑电网峰谷差并提升充电站运营效率。该方法体现了智能优化算法在电力系统调度中的实际应用价值。; 适合人群:具备一定电力系统基础知识和Matlab编程能力的研究生、科研人员及从事新能源汽车、智能电网相关领域的工程技术人员。; 使用场景及目标:①研究电动汽车有序充电调度策略的设计与仿真;②学习蒙特卡洛模拟与拉格朗日优化在能源系统中的联合应用;③掌握基于分时电价的需求响应优化建模方法;④为微电网、充电站运营管理提供技术支持和决策参考。; 阅读建议:建议读者结合Matlab代码深入理解算法实现细节,重点关注目标函数构建、约束条件处理及优化求解过程,可尝试调整参数设置以观察不同场景下的调度效果,进一步拓展至多目标优化或多类型负荷协调调度的研究。
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值