Visible Lattice Points SPOJ - VLATTICE

题意:

有一个n∗n∗n的三维直角坐标空间,问从(0,0,0)看能看到几个点。

题解:

本题是二维的一个升级版,升级成三维
用莫比乌斯反演来做

代码:

#include <bits/stdc++.h>
#include <unordered_map>
#define debug(a, b) printf("%s = %d\n", a, b);
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> PII;
clock_t startTime, endTime;
//Fe~Jozky
const ll INF_ll= 1e18;
const int INF_int= 0x3f3f3f3f;
void read(){};
template <typename _Tp, typename... _Tps> void read(_Tp& x, _Tps&... Ar)
{x= 0;char c= getchar();bool flag= 0;while (c < '0' || c > '9')flag|= (c == '-'), c= getchar();while (c >= '0' && c <= '9')x= (x << 3) + (x << 1) + (c ^ 48), c= getchar();if (flag)x= -x;read(Ar...);
}
template <typename T> inline void write(T x)
{if (x < 0) {x= ~(x - 1);putchar('-');}if (x > 9)write(x / 10);putchar(x % 10 + '0');
}
void rd_test()
{#ifdef ONLINE_JUDGE
#elsestartTime = clock ();freopen("data.in", "r", stdin);
#endif
}
void Time_test()
{#ifdef ONLINE_JUDGE
#elseendTime= clock();printf("\nRun Time:%lfs\n", (double)(endTime - startTime) / CLOCKS_PER_SEC);
#endif
}
const int maxn=1000020;
int mu[maxn];
int prime[maxn];
int vis[maxn];
int cnt=0;
void get_mu(int N){mu[1]=1;for(int i=2;i<=N;i++){if(!vis[i]){prime[++cnt]=i;mu[i]=-1;}for(int j=1;j<=cnt&&i*prime[j]<=N;j++){vis[i*prime[j]]=1;if(i%prime[j]==0)break;mu[i*prime[j]]=-mu[i];}}
}
int main()
{//rd_test();get_mu(1000000);int t;read(t);while(t--){ll n;cin>>n;ll sum=3ll;for(int i=1;i<=n;i++){sum+=mu[i]*(1ll*n/i)*(n/i)*(n/i+3ll);}cout<<sum<<endl;}//Time_test();
}

Visible Lattice Points SPOJ - VLATTICE相关推荐

  1. poj 3090 Visible Lattice Points(离线打表)

    这是好久之前做过的题,算是在考察欧拉函数的定义吧. 先把欧拉函数讲好:其实欧拉函数还是有很多解读的.emmm,最基础同时最重要的算是,¢(n)表示范围(1, n-1)中与n互质的数的个数 好了,我把规 ...

  2. POJ 3090 Visible Lattice Points 【欧拉函数】

    <题目链接> 题目大意: 给出范围为(0, 0)到(n, n)的整点,你站在(0,0)处,问能够看见几个点. 解题分析: 很明显,因为 N (1 ≤ N ≤ 1000) ,所以无论 N 为 ...

  3. POJ 3090 (欧拉函数) Visible Lattice Points

    题意: UVa 10820 这两个题是同一道题目,只是公式有点区别. 给出范围为(0, 0)到(n, n)的整点,你站在原点处,问有多少个整点可见. 对于点(x, y), 若g = gcd(x, y) ...

  4. POJ_3090.Visible Lattice Points

     1 /*  2  *  3  *因为是关于对角线对称,所以只需要求解矩形的下半部分在乘以2就可以了  4  *满足题意的点对的要求是x和y互质,可以使用欧拉函数  5  *需要注意的是特殊的点对(1 ...

  5. SPOJ - VLATTICE

    链接 题意:三维平面,找从(0,0,0)看(n,n,n)能看到的点 题解:很明显就是求gcd(i,j,k)==1的(i,j,k)对数,改一下公式即可,记得要算平行坐标轴的三个平面,还有含0的三个坐标 ...

  6. SPOJ-VLATTICE Visible Lattice Points-莫比乌斯反演

    需要将问题分解一下. 我们需要求的是这个立方体从(0,0,0)能看到的点的个数.可是在三个含有(0,0,0)的面上我们没有办法和其他的一起进行分析(因为含有坐标是0,而我们的数论工具都是从1开始的), ...

  7. 莫比乌斯反演 做题记录

    来自Peterwuyihong 的题单. 前置知识 前置芝士1 数论分块 UVA11526 H(n) P2261 [CQOI2007]余数求和 P2260 [清华集训2012]模积和 其中有一个式子需 ...

  8. python乌龟画花_python turtle画花

    [LeetCode] Reconstruct Original Digits from English 从英文中重建数字 Given a non-empty string containing an ...

  9. ACM 数学类题目推荐

    转:http://blog.sina.com.cn/s/blog_6635898a0100magq.html 1.burnside定理,polya计数法     这个大家可以看brudildi的< ...

最新文章

  1. android webview mailto,Webview email link (mailto)
  2. 在UltraWebGrid单元格中绑定dropdownlist
  3. ASP.NET(C#) 四舍五入、进一法、舍位(取整,舍去小数,向负无穷舍入)函数
  4. ros构建机器人运动学模型_古月私房课 | ROS机械臂开发中的主角MoveIt!
  5. php访问服务器文件路径,PHP与服务器文件系统的简单交互
  6. linux操作系统进程间通信IPC之共享存储映射
  7. 学校为什么要单位接收函_签了三方,想毁约怎么办?这几点你必须要知道!
  8. python读取配置文件使用_如何使用Python3读取配置文件
  9. shell命令的退出状态码(exit status)
  10. Keil破解过程中注册机被识别为病毒文件
  11. 2021-03-31【Caltech行人数据集的 .vbb 转 .xml 】(学习记录)
  12. 红帽linux系统安装指导
  13. HCSE设计知识点50个
  14. 《图解HTTP》笔记
  15. 中国护照可以免签的10个旅游天堂国家(地区)
  16. Interspeech2022论文解读 | LODR:一种更好、更轻量的语言模型融合新方式
  17. 【PS】ps基础绘画球体
  18. 2022年linux面经必备(持续更新中)
  19. mysql根据字母和数字进行排序
  20. 【串口助手】Python从零开始制作温湿度串口上位机

热门文章

  1. 一个富二代仅凭“1+1”就压制了全世界的数学家两个世纪......
  2. 为什么那么好的女孩子还单身?
  3. 为啥饮料瓶大都是圆的,牛奶盒却是方的?
  4. 159个故事串起三千年大历史!这套“儿童版史记”太无敌了!
  5. 为什么说人工智能是一个大谎言
  6. 陶哲轩的10岁与30岁
  7. 搞笑动图:这些痛,只有程序员懂…
  8. java中有stdin_在java中听stdin的后台进程
  9. saiku 连接 MySQL_Saiku连接mysql数据库(二)
  10. 捷达vs7测试_捷达VS7——品质硬核!