水仙花数是指一个N位正整数(N>=3)，它的每个位上的数字的N次幂之和等于它本身。例如:153 = 13+53+ 33。本题要求编写程序，计算所有N位水仙花数。

好的，我可以帮你编写这个程序。以下是 Python 代码实现：

n = int(input("请输入位数N："))
start = 10 ** (n-1)  # 起始数值
end = 10 ** n - 1  # 终止数值

for num in range(start, end+1):
    sum = 0
    temp = num
    while temp > 0:
        digit = temp % 10
        sum += digit ** n
        temp //= 10
    if num == sum:
        print(num)

代码解释：

1. 首先输入位数 `n`，确定水仙花数的位数范围。
2. `start` 和 `end` 分别表示水仙花数的起始数值和终止数值。
3. 通过 `for` 循环遍历所有位数为 `n` 的整数。
4. 在循环中，通过取余和整除运算依次计算每个数位上的数字并求和。
5. 如果计算结果等于原数，则说明是水仙花数，将其输出。

希望能够帮到你！