个人如果做网站赚钱,要给公司做一个网站怎么做的吗,wordpress变慢,wordpress xml怎么开Visible Lattice Points SPOJ - VLATTICE
题意#xff1a;
有一个n∗n∗n的三维直角坐标空间#xff0c;问从(0,0,0)看能看到几个点。
题解#xff1a;
本题是二维的一个升级版#xff0c;升级成三维 用莫比乌斯反演来做
代码#xff1a;
#include bits/stdc.h…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 pairint, 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 maxn1000020;
int mu[maxn];
int prime[maxn];
int vis[maxn];
int cnt0;
void get_mu(int N){mu[1]1;for(int i2;iN;i){if(!vis[i]){prime[cnt]i;mu[i]-1;}for(int j1;jcnti*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;cinn;ll sum3ll;for(int i1;in;i){summu[i]*(1ll*n/i)*(n/i)*(n/i3ll);}coutsumendl;}//Time_test();
}